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SYNOPSIS 


This study is concerned with the problem of estimating 
the parameters and testing the equality of location parameters 
6j[ (i = l/2/»»»,K) of K(> ?,) two~parameter exponential distribu~ 
tions E(0^,a) based on K independent type II censored samples# 
'Here the Common scale parsimeter o is assumed to be unknown# 

Type II censored sample is an ordered sample in which a known 
number of smallest (left) observations and/or largest (right) 
observations are missing* An ordered sample is obtained by 
rearranging the variates in an ascending order of magnitude# 


Most of the work in this field is done when complete or 
right censored samples are available# However, there are 
situations when some smallest observations are also not available 
In the present work, main attention has been paid to type II 
doubly censored samples* 


Let *..,X^^2o = 1^2, ...,K) be K 

ri+l ^i ®i 

independent ordered samples from EiQ ) v^ith r^ > 0, Sj^ ^ 0 
and rj_+l < nj^-Sj^* Based on these observations, the Least Square 
(LS) and Maxiinum Likelihood (ML) estimators of 6^ (i = 1,2 ,#*#,k) 


derived under the assumption ^ ^2 


and a are obtained* For K = 2, the LS and ML estimators are 

= 0^ =0* Some distribution 
theory results regarding these estimators are also obtained- 
A brief comparison of the LS and ML estimators is made by using 
the mean square error criterion* 
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For K = 2/ the test statistics based on LS and ML 
estimators are proposed for testing the null hypothesis 0^ = 0^ 
against one-sided alternatives- These statistics are equivalen 
to 

where 


cr’'^ = P]_/d is the pooled estimator of Oj, 


K ^r®i 
E { S 
i=:l j=rjL+l 

K 






d = S ( nj^~rjL'^Sj|_-l ) and K = 2» 
i=l 

Against the alternative 0^ > © 2 / large values of T lead to the 
rejection of the null hypothesis- The null and non-null 
distributions of T are derived- Some exact and approximate 
critical points of T are tabulated- Power values along v/ith 
a normal approximation are also tabulated- From thiS/ it is 
concluded that the test T is more sensitive for the left 
censoring than for the right censoring- 


The tost statistics ^ based on LS estimators^ 
ML estimators and LR test procedures respective! y^ are propose< 
for testing 0^ = ©2 against 0^ © 2 - These are given by 

Vj_ = !T~qj_l 


(i = 1,2), 
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and for Y > 0, 


X = 




^ ^ Hh JC 'q a -p /\ 

Y 1 ^ c7^+fY-,p)^ exp(~P/a^) 


{ d'*'" a^+(r^+f)Y~P}^^'^ {P-.d* + r2Y} ^ 


wh 


ere = b^-b^, q 2 = 


r.+l 


bj_ = E ■(nj_~j+l)'''^ (i = 1,2), 
j=l 


m 


^ = log { nj/(nj_~rj_) } (i = 1,2), 

a = P^/d'-' , d* = d+2, 
f = n^+n^ - r^-^r2, 

R = { 


fl n ^'^ Cn - r ^) 

1=1 


■( I ' m ) 


) exp (a*). Y = 


x(^)+s. x('() -(n.-r.) x()-) } 

• ^ ^ J n.~s. ^^1 r-+l^ 

i=l j=ri+l 11 1 


and is -the solution of the equation 


V^o 

e ° = {1 + 


^ 2 ^ 


d* d 


^ i/U +--; r 7 


Hi "'" 


d* + fY~P 


The expression for X given above is for r^ >0, r 2 > 0-- 
Considerably simpler expressions are obtained if r^ = 0 an^or 
r 2 =0# For Y < 0, X is obtained by replacing t.^,n 2 /r^/r 2 by 

n„,n. ,r<,, r. respectively in the above expression* 

JL ^ « 
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The exact critical points and power values of the tests 
and are evaluated* Due to the complex nature of X , only 
simulated critical points and power values of K are tabulated# 
On the basis of power calculations, it is concluded that the tesi 
V 2 is somewhat more biased than and there is very little 
difference between and X • Since the statistic X is far moj 
complicated than , use of is recommended# 

For testing K(> 3 ) populations, two test statistics are 
proposed. We first consider the case of right censoring, so the 
the smallest observation of each sample is available* For 
simplicity of notations let Xj_ = (i = 1 , 2 , and 

X(^) = min( X^,X 2 / # » *, X-j^) » The statistic given by 

Tj_ = {X^~min(X2,X3, ##-,Xj^)}/a'=^ 

is proposed for testing ; 0 ^ = ®2 ~ *** ~ ® against 

Hi : 81 > max( 02 / 03 / • » statistic 

T„ = { max (x. ) ~ min (X. )}/a* 

is proposed for testing against H 2 : at least one is 
different from 0# The exact critical points and power values 
of the tests Ti and T 2 are evaluated# For K = 3 case, the 
performance of statistic Ti is studied for different combinatiP 
of ni,n 2 /n 2 and d# It is observed that, the test Ti is more 
sensitive for changes in ni than changes in n 2 /n 2 and d* For 
K = 3 and equal sample size case, the performance of Ti and T 2 


XV 


along with that of and is studiod* Tho statistics and 
arc given by 


and 


U. = [ max {n,.(x.-x.)^ n, ( X.~X . ) } 1/da * 
2<i<K ^ 1 i 111 


K 

^ 4=2 

i=l 


n^( X^-X^ ^ j )/{ ( K~1 )a*} . 


The test is actually the LR test and has been proposed by 
other #jathors • On the basis of power calculations carried out, 
T^ is recommended for testing against and the LR test 

statistic is recommended for testing against "^ 2 * 

If the smallest observation is missing in each sample of 
size n, then the following test statistics are studied ; 


V = CX^^^ - min (X^^^)}/a* 

2<i<K ^ 

for testing against 

= { max (X^^^) ~ min (X^^^)}/a* , 
l<i<K ^ l<i<K 

and 

Vg = [ max { (X^^^~X^^ (X^^^~X^^^)}]/a* 

2<i<K - z z 

for testing against For K = 3, exact critical points 

and simulated power values of these tests are tabulated* On the 
basis of these calculations is recommended for testing 

against but nothing can be said regarding the preference of 
V 2 over as in some cases V 2 performs better and in some cases 
performs better* 



CHAPTER I 


INTRODUCTION AND SUMMARY 


1 *1 Scope * 

Exponential distribution is often proposed for modelling 
the lifetime distributions of items like electronic components, 
mechanical breakdowns, light bulbs etc* [ Davis (1952); Epstein 
(1958)/ Proschan(l963)/ Nelson (1975)]* In a two-parameter 
exponential distribution, the location parameter is interpreted 
as the minimum (or guarantee) time, before which no failures 
occur, and the scale parameter, as the mean life measured from 
the location parameter as the starting point* 

There are several situations, where the complete 
sample is neither available nor desirable* Since life— testing 
experiments are usually destructive, this limits the number of 
items to be tested (Sinha and Kale 1980, p* 18)* Moreover, in 
the ordered samples frequently found in biological data either 
some smallest and/or some largest observations are not available 
[ipsen (1949)]* An ordered sample is obtained by rearranging 
the variates in an ascending order of magnitude. In an ordered 
sample, if a known number of smallest (left) values or largest 
(right) values or both are missing, then such a sample constitutes: 
a type II censored sample# 

Several problems dealing with the estimation and testing 
of hypothesis of a two— parameter exponential distribution are 
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discussed in several books^ for example see Sarhan and 
Greenberg (1962)^ Mann^ Schafer and Singpuirwalla (1974)/ Bain 
(1978)/ Sinha and Kale (1930)* Conpletc and right censored 
samples have been considered by many authors (Walsh (19&0), 
Halperin (1952)/ Epstein and Sobel (1953)/ Hogg and Tanis (1963)/ 
Grubbs (1971)/ Kumar and Patel (1971 ), Dubey (1973)/ Weiaman ^ ^ 
(1973)/ Khatri (1974), Pemg (1978), Mathai (1979), Regal (i960), 
Bhattacharyya and Mohrotra (1981), Hsiah (1981), Goria (1982), 
Mohrotra and Bhattacharya (1982), Singh (1983), Singh and 
Narayan (198 3) etc*) * 

Although the left censored samples have not been 
considered so thoroughly as the right censored samples, there 
are situations in vAiich some smallest observations are not 
available* For example, in e3g>erimental biology, n animals 
are tested for antibodies after a certain period of time* Only 
(n~r) of these samples contain measurable aniounts while r of 
the animals develop the antigen at a level too low for measurement 
by the prevailing technique (Ipsen, 1949) • This gives rise 
to a censored sample from left* Another example where the 
smallest order statistic is difficult to observe is the failure 
time of human kidney* It is not easy to tell the failure time 
of one kidney since noticeable syraptoms occur only when both 
kidneys fail* Similarly, in a transistor set four battery 
cells may be used- Failure of only one coll may not affect 
the performance, and hence the first failure time may go 
unreported* But when two cells fail, the transistor may not 
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xvork properly and the second smallest order statistic becomes 
the first available observation • 

Based on doubly censored samples/ the Least Square (LS) 
estimators and Maximum Likelihood (ML) estimators of the 
parameters of a two-parameter exponential distribution were 
derived by Sarhan (1955) and Kambo (1978) respectively# Recently/ 
Tiku (l98l) and Khatri (l98l) have considered the problem of 
testing equality of location parameters of two-paramoter 
exponential distributions based on type II doubly censored 
samples* 

One of the problems arisi-ng in life-testing experiments 
is tho estimation and the comparison of minimum (or guarantee) 
time of the items* Related to this problem, tho present v/ork 
is concerned with the problem of estimating the parameters and 
testing the equality of location parameters of K(> 2) two- 
parameter exponential distributions# Here tho scale parameters 
are assumed to bo equal taut unknown# In this thesis# main 
attention has been paid to type II doubly censored samples* 

The following topics are studied in this thesis : 

1# Comparative study of the I^IL and LS estimators of tlie 
parameters # 

2,#. Testing equality of location parameters against one-sided 
alternatives • 

Testing equality of location parameters against two-sided 
alternative* 


3 
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Generalized statistics for K right censored samples* 

5* Generalized statistics in equal sample case^ when one 
observation is missing on the left* 

The results obtained are compared with the existing 
results* Suitable graphs and tables are provided to support 
the thvjory, v/herover necessary* The notations vjhich are used 
consistently in the text arc givt^n in the next section and th<i 
subsequent sections describe briefly the above mentioned topics* 


1 • 2 • Notations and abbreviations « 


As far as possible random variables will be designated 
hy upper case letters, and their realizations (observations) by 
the corresponding lower case letters* 


(i) 

1 


< X, 


(i) 





Xl( i) = (X 


( i) i) 


X 


(i) 

n..~s 


) 


ith ordered sample of size 

nj^ (i 5= 1/2/***,K) with superscript 

i dropped if there is only one sample* 

L. 

ith type II censored sample in which 
first r^ and last s^j^ observations 
are missing, where rj_ >0, > 0, 

^and r^L+l < 


P^(x),P(x) = P [x < x] 
. f^(x), p(x) 

E(X) 

Var(x) 

mse(t) 

se(t) 

LR 


cumulative distribution function of X 

probability density ftinction of X 

mean of X 

variance of X 

mean square error of T 

standard error of T 

likelihood ratio 
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KP test 
LS 

pdf 

jpdf 

cdf 

•w»r*to 

MVU 

DP 


E(e,a) 
N(u, n^) 


Kumar and Patel (1971) test 
least square 
maximum likelihood 

independent identically distributed 
probability density function 
joint pdf 

cumulative distribution function 

with respect to 

minimum variance unbiased 

degrees of freedom 
two-parameter exponential 
distribution with location 
parameter 0 and scale parameter c 

"normal distribution with mean ju, 

, 2 
and variance 17 

"central chi— square distribution 
with y DP 


a 


f 


8 * 


a 


o' 




A 



d 



1 

BCa^b) = X t^”^(l-t)^"'^ 

o 

1 (a) = X t^~^ dt, 

o 


~LS estimators of a and 0 
respectively 

LS and ML estimators under the 
hypothesis 0^^ ~ 

is distributed according to 

is approximately equal to 

is identically equal to 

dt,a >0, b > 0 

a > 0 



6 


Q (xis) =; dt/(p-l)i / X > 0, p = for s > -1 

^ X 

= e3<p £~x(1+s)Hx(1+s)} V{ jl (1+s)^} 
j=0 

X 

L (x|s) =/ ^t/Cp-l)! , X > 0, p = 1,2^... 

^ o 






f Xp/pl / s = 1 

l/il~s)^ ~ 2 ejtp£-J-x( 1-s)) x^( l~s)^'“^/jl , s 1 

j=0 
K 

= E ( n •*-rj“'S . ) 

x=l 

K 

= E ) = d* - K 

i=l 


2 


(i) 


U “"i 
.(i) 


2 { S X''. ' + ^ 

i=i j=r,+i J ^ 




1 *3 » Comparative study of the ML and LS es-timators of -the 
parameters » 

Let x:^4.i/X:^+ 2' ’ *'^n~s ^ ‘^yp® censored sample 

from an exponential distribution E(0,a) with pdf 


f(x?8/Cr) =: exp {-( x— 0 )/ct} , 0<a<oo, 0<x<«»# 

SuXhatme (19 36) obtained the best unbiased ostirtators of 0 and 
0 based on a complete sample (r = 0, s =0 case) of size n* 
Lloyd (1952) discussed the technique of estimating 0 and a by 
applying general LS theory to an ordered sample* By applying 



7 


Lloyd's method Sarhan (1954/1955) and Greenberg and Sarhan (1962) 
have obtained the best linear unbiased estimators of mean and 
standard deviation for daable and middle censoring case as well* 

Let li = 1,2,. ..,K) be K 

independent type II doubly censored samples from E(e^/a)# The 
LS estimators 6'^ and a* of 6jj^(i = 1/2/**«/K) and a are obtained 
in Section 2*2* These are given by 


r^tl . 

whore b, = I (nj^-j+l)“ 

j=l 


1,2/.*.,k) and a* = 

K 

(i = 1/2, -»*,K)/ d = E (n .'-rj~s .—1) 

1=1 ^ ^ ^ 


K 




and P- = E { y X^. + s. 

1=1 jlri+l J ^ "i-®l 




In Section 2 »3/ it is shown that (i = 1,2, •♦»/K) and 

IL ^ 

are independent, and 2da*/a has ^ ^2d distribution » 


ML estimators of the parameters for a complete and 
right censored sample were dcirivod by Sukhatmo (19 36) and 
Epstein and Sobel (1954) respectively* For a doubly censored 
sample, ML estimators wore discussed by Tilcu (1967) • He has 
obtained the ML equations and has mentioned that, those equations 
do not have oqDlicit solutions* Hence, ho obtained a modified 
ML estimator of the paramo bers* Using those equations, 

Kambo (1978) has obtained 03<plicit expressions for the ML 

A A 

estimators* In Section 2*4, the ML estimators and O of 
(i = 1,2, -••/K) and a are obtained. These are given by 
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®i ~ ^r^+1 ■*' ^ “ 1/2/**»,K) and a = 

■¥t 

where d = d+K- 

Kambo (1978) compared the minimum variance unbiased 
(MVU) estimators with the ML estimators/ •v\flien a doubly censored 
sample is available » He has shown that for a single sample 
MSE(a^) > MSE(a) and verified numerically that MSE(6^) can be 

A 

greater or less than MSE(6;, By this, he concluded that/ some 
times I/IL estimators are better than MVU estimators* Similar 
type of comparison is done in Section 2 *5 and it is concluded 
that in general for K > 3/ MSE(a^) is less than MSE(a)* However/ 
for K < 2, MSE(a'^^) is greater than MSECo), 

Epstein and Tsao (1953) derived the ML estimators 
under tho hypothesis 0^ = ©2 = ® ^or right censored samples* 

In Section 2*6/ the LS estimators and ML estimators for doxobly 
censored samples are discussed under the same hypothesis * A 
brief comparison has been done between LS estimators and MIL 
estimators in this case also* 

1 «»4 * Testing equality of location parameters against one-sided 
alternatives • 

For the right censored samples from two populations 
KuiCiar and Patel (1971) have proposed a test statistic for testing 
Hq : 6^ “ H 2 : A ^ 2 “ Weinman et ^* (1973) 

extended it for testing H^ against one-sided alternative 
* ®1 ^ ®2* statistic is W = (X^^ )/a*/ where 
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Cf* « pooled estimator of a • They obtained the 

critical point c as 

r d[{ nj^/(nj_+n2 ^'^^-l] /ng if < an2/(l-*cc) 

c ! 

|_ d [l-{n 2 /(n|j^+n 2 ) (l—0!;)}^'^'^]/n^ otherwise, 

where a is the chosen level of significance. Note that c > 0 
for > an^/ (l-oc.) . They obtained the power function P(<P) of 
W for c > 0 and <f> = (ej-e^Vcr > 0 as 

P(<p) =s 1~{ ( n. tn,, ) } [(n-+n^) S (d<p/c)Vil 

^ i=0 

^ d-1 

-{n^/Cl+njc/d) i £ {<p(d+n2c)/c} Vi! ] 

i=0 

~n-,<p <3.-1 . 

-[rue ^ 2 {<P ( d-n- c)/e} Vi ! ] 

i=0 

. l/{(nj^+n2)(l-nj^c/d)V, 

provided that n^c d* For the case n^c = d, the last term in 
P ( <P ) becomes 

-n- cp ^ 

-n 2 e (d<P/c) /(n^+n 2 )dl» 

For c < 0, P(<^) ) is given by 

P(<P) = l~n 2 expC-n^^cp) (1— n^c/d) /(n^+n 2 )» 

In Chapter III, a test statistic based on ML and LS 
estimators is proposed for testing against : 0^ > 02 * 
This test is equivalent to 

T = i4lli - 
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The distributions of the statistic T under and are 
derived- Approximations for null distribution in terms of 
Student's t and normal distributions are studied- Approximate 
critical points obtained from above approximate null distributions 
are also tabulated along with the exact critical points- It is 
observed that the normal approximatiori is better if ^2' 

otherwise Student's t approximation is better- Some exact and 
normally approximated power values arc tabulated- Also, the 
variation in power diie to different combination of r^^ and ^2 
is plotted* With power function as a base, it is concluded 
that the test T is unbiased and it is more sensitive for the 
variations in r^ compared to the variations in r^* Furthcsr, 
the performance of the test is not seriously affected for 


variations in right truncations for fixed values of 



Epstein and Tsao (195 3) discussed the LR test 
procedure for testing various t 3 fpes of hypotheses based on right 
censored samples- Kumar and Patel (l97l) proposed a test based 
on 1 for testing : 0^ =02 against 

Dubey (1973) and Weinman ^ ad- (1973) derived the power ftinction 
of KP test- Weinman et al- compared the KP test with LR test 
and they concluded that in general LR test performs better than 
KP test- The power function P(<iP) of KP test with critical 
point c was obtained by Weinman et ^* For <p = (62~6j_)/d > 0 
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it is given by 

d-1 

P((p) = l~exp(”d<p/c) 2 {6ip/c)^/±l 

i=l 

+ exp(-nj^<p){ (l+n^^c/d) ^~( l-n^c/d)~^} 

^3. X 

+ { n^/(n2_+n2 )} (l+njc/d)*"^ ^-dfP/c ^ { q3( n^+cVc)} 

i=0 

+ {n2/(n^+n2)} (l-n^c/d)''^ 2 { <P(~nj^+(yc)} Vi* 

for l-n^c/d /= O, However, for 1-n^c/d = O the third and last 
terms of P(<p) become respectively 

2 n 2 /(nj^+n 2 )} exp(-n^ <p)and ~n 2 (d«P/c)'^ e3^(~n^(p)/£ d I (n^+n 2 )}'» 

The power P(9) for <P < 0 is obtained by interchanging n 2 and n^ 
and evaluating the above expression for 1 q)l • 

Recently, Tiku (1981 ) considered the problem of 
testing : 0^ = ®2 ^2 * ®1 ®2 doubly 

censored samples- He proposed the test statistic 

u = 

rj^+1 

for testing H^ against H 2 and obtained its null distribution as 
^2 . r 

fy(u) =h[ 2 (~1)^(.^') B(f-i-j,r 3 _+l){l+h 2 (j)u}“^~^ 
j=0 ^ 

^1 . r 

+ 2 (~ 1 )^( .^)B(f+j,r 2 +l){l-i-h- ( j)u}~'^"'^ ] , u > 0, 

j=0 ^ 

2 

whore H = 7T { B( n^-r^/ r^+1} hj_(j) = ( n£~rjj^+j )/d (i = 1,2,), 
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f = nj^+n 2 '-rj^~r 2 and d = f- 2 -«S 2 ~S 2 • 

Khatri (l98l) has derived the non-null distribution of U as 
follows : 

g^Cu) = H S (~l)^'(^^)B(f+j,r 2 +l) S (d(p/u)^{l+h. (j)u} 
j =0 i=0 

^2 . r* 

. exp(-d(p/u)/il + H E (-1 )^' ( ? )b( f+j,r- +l){ l+h„( j )u}~*^''^ 

j _0 J 

# ej<p {•-h 2 (j)d(p} - H (-l)^( j^')B(f+j,rj^+l)£ l-h 2 ( j)u}“^~^ 

d . . 

* exp{-dfp/u} [ E (d(j)/u)^{ 1-h^ ( j )u} / i| I , u > 0, 

i=l 

where <p = (0^-02)/h >0# This final expression is slightly 
wrong due to some integration errors • Further^ he has not .taken 
into account the singularities at u =s l/h 2 (j) for j = 0,l/2/*-»/ 
(for (P > O), as has been done by Weinman et (197 3) for 

right censored samples# We give correct form of this expression 
in Section 4 #2# 

In Chapter IV, two statistics defined by 
= |T-q^l and V 2 = 1 T-q 2 1 based on LS and ML estimators 
respectively are proposed for testing against H 2 / where 

- 4212^'^*- 11 = bl-^2- 12 = 

r^^+l 

b^ = E (nj:~j+l)~^ (i = 1 , 2 ) and m^ = log{ ny'(n^-rj_) } ( i=l, 2 ) . 
i=l 

The null and the non-null distributions of these statistics are 
(derived# As a special case, the non— null distribution of Tiku's 
statistic(U) is also obtained# 



13 


The LR test statistic X is derived for testing 
against For r^^ > 0, r 2 ^ 0^ Y > 0, it is given by 

XT *4*X* 

’ S ^ *A =^1+* ^ f. ^2 ' 

{draQ+(rj_+f)Y~P} ^ £P-d* 


n, n. 


“1 “? -(n.~r. ) -(n^-r,) ^ 

where const = ^^2*" ^2^ * e:xp(d ),a = P^/o, 


x5.^^+SiX^^^ -(ni-rj.)x^^|. ], 

i=l-j=r.+lJ ^ ^i“®i ^ 


P- = E [ S 


Y = X^ -X^ , P = P^ +( n«— r,,)Y and a_ is the solution of 
r 2 tl 12 2 o 

the equation 


Y/a^ 


1 + 


r^Y 


^ A 

P-d a 


'/ 


1 t 




d* Sq + fY-P i 


For Y < 0, X is obtained by replacing n^,n 2 gr^fr 2 by Uj/n^/r^/r^ 
respectively in the above expression# 


The caritical points of the tests ^ 

tabulated# The corrparativo performance of all the four test 
statistics is studied- Since the distribution of X is not 
easy to derive, only simulated critical points and power values 
are used- Using the power function as a base, it is concluded 
that the test U is more biased than and V^/ is more biased 
than ,and there is very little difference in the power values 
of X and • Since the statistic X is very complicated, while 
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•the statistic is considerably simple^ we strongly recommend 
the use of the test statistic in such situations* 

Generalized statistics for K right censored samples* 


Khatri (1974) derived -the LR test and two test 

statistics and U^/ using union intersection principle for 

testing : 6^=B 2 =* - 6 against II 2 i at least one 6^ 

is different from 9, based on K independent right censored 

samples “ l/2,*»*/K) from E(0j^,a)<» For 

^^i i 

simplicity of notations^ let (i = 1/2,***/K) and 

= min(X^,X 2 / *• »»^Xj^) - Then Khatri^ s statistics are given 
by 


U. 


K 

S n 
i=l 


^(X-X(i) 


)/da^. 


U. = r max {n.(X.-X,. O}-)/d0* 
1<1<K ^ ^ 


and = [ max {n^ (X.-X. ),n-, (X^-X. )} 1/da** 

^ 2 <i<K 1 ^ T 1x1 

He has obtained the null distributions of and Further, 

he has dj.scussed thoir non-null distributions without carrying 
out any power calculations* Singh (193 3) also discussed the 
LR pixjceduro for testing against H 2 * He obtained the LR test 
as 


U 


4 


K 

s n,(Xi- X( ^ ) )/{ (K-1 )a*} = dU^/ ( K-1 ) 


and has shown that, the null distribution of is ^2(K-l),2d ' 
However, he has not studied the power function of U^* 
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It is reasonable to study a test for testing against 

the alternative It : 0^ > max (0j)» For this purpose, in 

2<i<K 

Chapter V a test statistic given by 

T ^ “•min ^ ^ 2 ^ ^3 f * * " ^ ^ ^ /d* 

is proposed* Although there are several tests (as mentioned 
above) for testing against yet we propose another test 

based on 

= { max (X. ) - min (X. )}/a** 
l<i<K 

Chapter V is mainly devoted for studying the performance of 
for K = 3 and different combinations of n^/n^/n^ and d* 

For comparison purposes the performance of statistics 
and for K = 3 and n^^ = n 2 = n^ = n is also studied f Note 
that for equal sample sizes U 2 = nT^/d* The necessary critical 
points of those tests are tabulated* Some exact power values 
of the tests and are tabulated for points in the 
parametric space satisfying 0^^ > '> ©g* Simulated power 
values of T^,T 2 ,U 2 and are also tabulated, since the 
expressions for povrer functions of and provided by 

Khatri (1974) aro extremely complicated* Prom theso calculations 
of power values, it is concluded that 

( i) the test T^ is more sensitive to changes in n^ 
compared to 

(ii) the power of the test is considerably higher than 
that of other throe statistics. 
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(iii) the test T2 performs slightly better than if 
( 9 2"*® 2^ small, otherwise reverse is the case, 

( iv) the test T2 performs better than if (6^-02) is 
largo, othor\^^ise reverse is the case* 

Finally, the LR test (or equivalently is 

recommended for general alternative hypothesis since its 

critical points are easy to evaluate from the P-distribution . For 
the specific altomativas like the statistic Tj is recommended, 
since its critical points are available in a compact form and 
its power is considerably higher than that of ether tests* 

1 *7 * Generalized statistics in equal sample case* when one 
observation is missing on the left * 

In Chapter VI, tests for the equality of location 
parameters of K populations are considered, when the smallest 
observation is missing and atleast second smallest observation 
is available in each sample of equal size n* Here the tost 
defined by 

~ min 
2 <i<K 

is proposed for testing : 0^ t= ©2 = • » » = 6 ^. against 

^1 * ®1 ^ max (62/63/ •• • Simil ir to the previous section, 

the statistics V2 -.md given by 

Vo = { max (xl ~ min {yL^\/a* 

3 .<i<K 

V- = [ max { (4^^~X^^^),(X^^^-X^^^)}]/0* 

2 <i<K ^ ^ z ii 


and 
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are proposed for testing against : atleast one 0^ is 
different from 6# Compared to two-sample case (K = 2 )^ the 
LR test is much more complicated even for K = 3 » Consequently, 
this has not beesn studied at all. 

For K = 3 / the exact critical points of all the three 
tests are tabulated for some selected values of n and d- Since 
it does not appear simple to evaluate the non-null distributions 
of those statistics, the power of those tests are calculated by 
Monte-Carlo techniques for 0^ > 0 ^ > 63- From this study, we 
have recommended for testing against a specified alternative 
like Hj^,and for testing against the alternative H2/ the test 
is recommended if ©^^-02 is ejoected to be very small, otherwise 
is recommended* 



CHAPTER II 


ESTIMATION OF THE PARAMETERS AND BASIC DISTRIBUTION THEORY 

2«1* Introduction « 

In this chapter, derivation and the comparison of 
the Least Square (LS) and Maximum Likelihood (ML) estimators 
of tho location and the scale parameters of two-parameter 
exponential distributions are discussed* Some of the results 
are established, which are \ised in later chapters* 

Let X ,-/X X be a type II censored sample 

from an exponential distribution E(0,a) with pdf 

-1 

(2-1-1) f(x;e,a) = a ■ exp{-(x-e )/cj}, 0< a <«>, e<x<'»« 

Many researchers have investigated the problem of 
estimation of tho scale parameter a and the location parameter 
C* For a complete sample (case r = s = O) of size n, 

Sukhatme (19 36) obtained the best unbiased estimators of 0 
and a as 

n 

(2-1.2) 0* = X--aVn/ 0^= ( E X.-nX^)/d, 

j=l ^ 

where d = n-r~s-l = n-1 for r = s = 0# 

Lloyd (1952) discussed the technique of estimating 0 
and 0 by applying general LS theory to an ordered sample. 

Type II censoring on the right (case r = O) was considered by 
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Sarhan (1954)/ and he obtained the LS estimators of 0 and a 
as 


(2*1*3) 


n— s 

X^-a*/n/ = ( 2 Xj+sX^^^-nX^V^i* 

j —1 


The LS estimators of 0 and a based on a type II doubly censored 
sample was discussed by Sarhan (1955)/ and is given by 


r+1 . n-s 

(2.1*4) 0^ = ^/ =( _E^_^^X^.tsX^_3-(n-.r)x^^^}/d. 


Generalization for K > 2 type II doubly censored samples from 
E(0^/a) are discussed in Section 2.2* The necessary distribution 
theory results are obtained in Section 2*3* 

The ML estimators of 0 and a for a complete sample 
was obtained by Sukhatme (1936) as 


/V yy 

(2*1.5) 0 r. X^/ a = ( E Xj-nX^)/d* / 

j =1 

where d'^ = n-je-s ?= n for r = s = 0» For type II censoring on 
the right/ the ML estimators was discussed by Epstein and Sobel 
(1954)* This is given by 


A A 

(2.1.6) 0 = X- / a = ( E X.+sX ~nX.)/d*. 

j=l 

Tiku (1967) has derived the ML estimators of the parameters 
based on a type II doubly censored sample. He has obtained the 
following ML equations : 
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(2.1.7) ^ [l-f 

and 


( 2 , 1 . 8 ) 


£ -J— i±iy] « o 

"^r+l ‘ ■' 


e, 1 

S N ,1 _ 


+ S 

n n-s n P( z 


r ^r+1 ’ 

j 2 


r+l‘ 


r+l 


j = 0, 


where f(z) = exp(-2) and F(z) = l~exp(~z) are the pdf and cdf of 
Z = (X-0)/a respectively# According to Tilcu (1967)/ ML equations 
( 2 . 1 . 7 ) and (2-1.8) do not have explicit so3.utions (for r > O) 
due to the presence of the term f( 2 )/P( 2 ). Ho therefore obtained 
modified ML estimators and of 0 and a, which are 

g iven by 


XT 4*1 

(2-1 .9) 0 ^=X ^ - a , £ — ^r-r 

mod r+1 mod n~ 3 +l 

1 =r+l 

In equations ( 2 . 1 *7 ) and ( 2 -1 *8 ) / Kambo ( 19 78 ) 
eliminated ^ ^ V^( ) and obtained explicit expressions for 

the ML estiniators. These are given by 


( 2 ,1 .10) 


A 

0 = 


X 


r+1 


+ a log (l-r/n), 


A 

a 


a 


mod* 


Note that/ for r = 0/ equation (2»1*10) reduces to equation 
( 2 ,1 .6 ) which in turn reduces to equation ( 2 ,1 .5 ) for s = 0 • 


Prom equations (2.1.4) and (2.1.10)/ it is easy to see 


that 


r+1 


( 2.1 .11) 0^ = X_.^ ~ 0* £ 


r+1 


j=i 


/ o 


* d 0 
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In Section 2 *4, ML estimator's for K(> 2) samples are 
derived. A brief study of MSE of ML estimators and MVU estimators 
is made in Section 2.5. 


Epstein and Tsao (1953) considered the problem of 
testing equality of two exponential distributions based on 
right censored san^les. They derived the ML estimators under 
the hypothesis = ®2 = These arc given by 


( 2 . 1 . 12 ) ©^ = min 


and 


A 






a xsl 


3=1 


i 1 


whore d* = S (nj^-rj_~s j, ) with = r^ =0# For doubly censored 
i=l 

samples, tho LS estimators and the ML estimators under the 
hypothesis ©^=© 2=6 are discussed in Section 2.6. 

2.2. LS estimators of the parameters for K(^ 2) samples case . 

In this section, tho LS estimators for K samples are 


d er ived 


Dd. Lot ^rt+l'^r^4-2' ' ^ ^ 1,2, ...,K) be K 


independent type II censored samples from E(©jj_,a). Denote 

2 

E(x) = a 7 and Var(x) = o "D, 


where 


E(Xj^^) = 

Var(X^'^^) 


0i + 
= 

3 


c:.“' = s l/(n.-g+l). 


/_i) 

g=l 
j 


= S l/(ni-g+l)^'. 


g=l 
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x(i) ' 

ro 




e 


1 

CD 

1 

X == 

X( 2 ) 

(TO 

, X( i) = 

y.( i) 

/ r = 

ro 

j e = 

®2 

ro 

f 

m 

* 

» 

x(k) 

ro 

^ ro 

1 

i 

IT^ +2 

» X 

^(i) 

' to 

a 

^ ro 



A 


1 0 ^ # 
ro 


0 c(i)i 

rv) ro 

0 1 0 C( 2 ) 


ro ro 


O 0 * 

ro ro 


1 c(k) 


/ C(i) 


IT » “f*!, 

,(i) 

' t ^+2 


,( i) 
n.-s^ 


Cl*X(K+l) 

and rank(A) = K+i/ D is the dispersion matrix of X/a, and is 

ro ro 

given by 


D = 

ro 


5i 


i 2 

00 s 

0 

ro 


> £2 

0 9 0 

0 

CO 

i 

0 

ro 

0 0 0 

£k 

d*xd* 

(i) 

r^+1 

A±) 

^r^+1 

0h i 

’ ' ^r^+1 

( i) 
r^+l 

ri+2 

0 i 

A±) 

^ri+2 

’(i) 

^r^+1 

^(i) 

ri+2 

0 

.(i) 

00 d 

TX£-S 


dixdi 


K 


d* = 2 < 3 .*/ * 3 ^ = ^i“^i'”®i' 3 “ n^-s^/ i=l/ 2 ^--*,K; 

i=l 

and 1 and 0 are the matrices of suitable orders with all entries 

ro ro 

as 1 and 0 respectively* 



Following Lloyd (1952), the LS estimator of 'i in 
this setup is given by 


( 2 - 2 . 1 ) 


r 

ro 


pi- 


s' 

ro 


(A' Q A)‘'^(A' Q x) , 

iNi rsi fNj ro 


where Q = D 

ro CO 


-1 


Qi 0 

rsU. ro 

0 Q, 

CO tO'* 


0 0 . 

ro ro 


That is. 


col 


a. 

1 


# 

0 


o 

ro 

0 

CO 


gK 


and Qi = g-L 




0 


0 


(njL-rj^-l)^’+(nj_~rj^~2)^ -"(n^-r£~2) ♦ 0 

-( s 


0 


0 


» <• . i. s 


Hence, A' q A 

to ro 4ro 


1 

\ 

0 

•t # • 0 

;V“i 

0 

41 

l/a2. 

♦ *1 # 0 

iV®2 

1 

j 

# 

0 

0 

... 1/aj^ 

1 

1 ^c/ ^ 



_ 


f 

’"l 

^2 


1 

' Q 

^1 

^2 


i ^ 

) 


(K+1)X(K+1) 


r^+l 


bj 

1 

= 

ri+l 

1 

= E 

j=l 

(Hi- 




r • +-1 



3 . • 

1 

- ^r^+l 

1 

= E 
j=l 

(n^- 

'j+l) 


K 



K ■ 

Q = 

. d+ 2 (bt/a.), 

J A X o- 

d = 

E dj and 

■J* 


i=l 


i=l 


where 
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Now, for a part ion ed . matrix 


R S 

(TO ro 


g/ ^ 
L.<ro 


(where R is non~singalar) , 


the inverse is (see, for example Rao 1973, p« 33) 

,-11 


R~^ + F E“^ P' 

4N) CO ro CO 

ro CO 

IL« 

where E = T-S' R~^ S and P = R~^ S. 


~F E 

CO ro 


E 

CO 




CO CO CO fO CO 


CO ro ro 


Applying this we have 


(2.2 .2) (A'Q A)”^ = i 


CO CO CO 


aj^a+t^ 

^1^2 

... b3_!:^ 

bj^bj 

0 

A 

a2d-Fb2 

0 

0 

0 0 0 


0 

0 0 0 3.j^ClHhl3^ 

-bi 

-'"2 

... -b^ 


Now 


(A' Q x) 

CO CO CO 




/a 

^rj^-l-r K 
K 


^ + 2 b.X^ ^?-,/a . 

1 i^l 1 i'i+1 ^ 


“■^2 

C 

0 

0 

1 


ni-Si 


whore 


K ri ( * \ 

E [ 2 -(n^-'r.)xl„^^J .Substituting 

1=1 j =r j|_+l •• •’ - 


1 1 

for (A' Q A)*"^ and (A'Q x) in equation (2*2*1) we obtain for 

ro CO no CO ro ^ 

i =: 1/2/ 

el = V * <’* = "’l/’^ ' 


( 2*2*3) 



where 
( 2 
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j=l 

Note that for K = 1, ecjaation (2i2#'^) reduces to equation (2^3 - 4 ) 
obtained by Sarhan (1955)* 

The variance covariance matrix of r* is (A'q 

ro ro *v» ro 

From equation (2.2*2) wo have 


(2.2.5) VaHe"^) [a^+b^/d] (i = 1,2*..*/K), 

( 2 *2 -6 ) Var(a^^) = a 
and 

( 2 .2 . 7 ) Covar(6t,a''^) = -b^d V<a/Covar(e*,e^)=b^bjaVd( i/£j=l, 2/ . ^ 

Hence / 

( 2 .2 .8 ) Var(e*-eJ) = [a^+a 2 +(b^--b 2 ) • 

2*3. Distribution theory . 


In this section, the distribution of the statistics are 
discussed, which are main tool for deriving the distribution of 
the ijS estimators • 


Theorem 2.?.1 . 
distributed. 


XI) ^(2) 

r^ +1 ' rp +1 ' 


„(k) 

' r^tl 


and a* 


are independently 


Proof. The ipdf of X(l),X(2) , #,X(k) is 


f ( x( 1 ) / x( 2 ) , 



i=:l 


. x(K) 2jr • ^ ^ 

n^l [ l~exp{-(x^^|j-e j^)/a}] ^ 

X 1 


r^L* Sj^t a 


^i'"^i‘’®i 



where 

( 2 . 2 - 4 ) 
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r^+1 

b. =* £ . 

j=l ^ 

Note that for K = 1, equation (2»2»^) rec3ucGs to eq^aation (2-3 -4) 
obtained by Sarhan (1955). 

The variance covariance matrix of 7* is (A'q 

1-0 ro rs) 

From equation {2 •I* 2) we have 


(2-2.5) Var(e'^) = [ a^+b^/d] (i = 1,2,.»*,K), 

(2-2-6) Var(a‘^^) = a^'/d 

and 

( 2 -2 .7 ) Covar(e ) = -bj_a ^d, Covar(e J,ep=b^bja V<a( i/^j=i, 2/ - -,Y. 

Hence / 

(2-2-8) Var(e*-e|) = [a^4a2+(bj_~b2) V^. ] • 

2 - 3 • Distribation theory • 


In this section, the distribution of the statistics are 
discussed, -which are main tool for deriving the distribution of 
the tS estimators * 


Theorem 2-3-1 - 
distributed- 


Xl) ^(2) 
r^ +1 * r2 +i ‘ 


^(K) 

' rjr+1 


and a'^ 


are independently 


Proof- The 1pdf of X(1),X(2),- .-,X(K) is 

CO ro 


f ( x( 1 ) , x( 2 ) , 

ro ro 


K 

n 

i=l 


..,x(K);0 

n^i [ l~e3q){-(x^j]_3^-e j^)/a}] 


Til s^t 




^i 

— X 
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-SiVa 

i X 


n .-s. 

X X /■ 4 

E (x!.^-e.)/a} 

j =r^+l ^ 


For the ith sample (i 
t r an s f o rma t io ns , 


for (i = 1,2, 

1,2,.,.,K), considering the 


( 2.3*1) 


^r^+1 ~ ^^i"'^i^^^r^+l~®i^/^'' 


- (n^-j+l)(x^.^^-Xj^^)/a (j = rji^+2/ • - 


and noting that ^ 0, the corresponding inverse transformation 


xs 


S i) 

■i' 


’'rT+l = “ + ®1' 




.( 1 ) 


-( i) 


(i) 


^nf-Si " ^ ^rjL+l^^’^i’'^i^'^ "^r^+ 2 ^^ ^i'"^!"’^ ^ 




Hence, 


n .~s. 

^ V V Mi) Q \ J ^ o'! 

0 S y = £ (x ~e ) + sA:x -G . ) 

j=rj^+l ^ j=rj^+i J X X 


Following standard methods^, the Jacobian of the transformation 
is 


,.K s J ! o 

J = n 


XXX 


isrl (n^-rj_)! 
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The jpdf of (j = Tj^+l/ • • i ~l,2,»»'p,K) is then 

obtained as 

f(;g(l) / y( 2) .. y(K)) 


K 


^ f - ' .W./ //_ „ \t T i r ^ -S^K 


T\^ r j_ ! ’ ( r j_ | ^ "'^r ^ +1^ ^ ’^i~’^i ^ ^ 1 .<,^j ^ 


J=r ^+1 


( 2-3 -2) 


K 

= n 

i=l 


n.l 


.(i) //„ _ r -,(i) 


pnip^pi ®^f-yr,+i» 


i^-^ J 


^i-Si 


n [exp{-y^>^] for > 0 


(i) 


j =^^^+2 

This shows that aro independently distributed with the 

following density functions : 


( 2-’=' -^) 

and 

p(i) 'y’ = r-unr 

( 2 -3 -4 ) 

j 

for j = 


Hote thatz 

(2-3-5) 

K ’^i'^^i / . -v 

E s y ' ^ 

i=l j =rj_+2 ^ 

v/hich is 

free from yi^l- ( 

ri+i 


til r , 

Kn.ir ' ) i [l*-®5cp{-y/(njj_-~rjL)}] exp(-y),y>p 


,K, 


da'^/a. 


From equations (2*3-2) to (2-3-5) and transformation (2*’^-l)^ 
it is easy to see that 

independently distributed* This completes the proof of the 


theorem 



2v3-l» 2d0*/a has a chi~ square distribution with 2d 

degrees of freedom (DF) i 


££2££* making the transformation = 2y^. j=r.-:-2. ....n.-s- 

J J 1 xx 

X = in equation we see that Z^. are 

iid chi-square variates with 2 DP, Consequently, 

,, n.-s. 

K 11 , V 

S 2 z\^ = 2da*/c^ 

i-l j=r^+2 ^ 

has a chi-squaro distribution with 2dDF, 

In particular, the pdf of Vi = da*/a is given by 

(2 *3 *6) f(w) = c"*^/ ITd) , w > 0, 

Corollary 2«3«2« The pdf of (i = 1, k) is given by 


(2.^.7) f 


X 


( i) 


[l-expC-Cx-© . )/a} ] 

(x) = exp£-(nj^~r^) (x-9 




B ( r j^+1 , nj_~rjL ) c 

for X > 6^, where B(rj^+1, nj-rj^) = rj^l (nj^-r^-l ) i/rij^U 

Proof » The proof follows directly from the marginal pdf of 

( i) 

^‘r-+l fTOm equation (2 -3 **3) and the transformation given by 
equation (2*3*l)» 

The distribution of the LS estimator 0* of 0 isl, 2, . * •,K) 
is given in the following corollary *. 

C orollary 2«3*3 « The pdf of = (0*~e^)/a (i = l,2,.-o,K) 

, , , Pi (y)/ y < 0 

(2.3^) f„ (y): ' 


xs 



29 


where 


1 ri d-1 {-(rij^-r.+j+d/h. )y}^ 

P^Cy) = G (~1)^( ^)exp{dy/b.} s i— i 

j-G ^ h=0 h 1 ( n^~r^+j+d/^;i^^ 


P 2 (y) = G j^)exp{~(njL-*r£+j )yl/( n^-rj_+j+(Vbj_)^/ 




G = (d/bj_) /B( rj^+l^nj^-rj; ), = 2 (n.~j+l)* 

j =1 


and is given in equation (2*2*’^)« 

Proof » The proof follows immediately on making suitable 
transformation and using Theorem 2*3 *1 along with equations 
( 2 .3-6) and (2*3*7) * 

2-4_* ML estimators of the parameters for K(> 2) samples case* 
Let X(i) (i = 1/2/ •••/K) be (K > 2) independent type II 

ro 

doubly Censored samples from E(0^/a)» Then the likelihood 
function is 


L ( 0 ^ / 0 * ♦ * / O t ^( 1 ) / !x^( 2 ) / * • */ ^( k ) ) 


K n . ! 
77 


n.j a ^ ''^^r.+l“®i^'^‘^ r. 

-i: [l~e ^ 1 


X 1 


XI i ) 

+ S il - 0i(i=l/ 2/ > 0, 

j=rj_+l ^ i 

v;here d* = nj_~rx'-Sx* For i = make the substitution 


3 


(x^.^-ex)/c7 (j = rj_+l/ ^«.,n^~Sj_), 
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then we have 


^(©1,0 2^ »• 1 ^(1 ) , 2( 2) / O » •, z(K) ) 

„(i) 


K n . 1 a 

X * 




i=l 


[1- 


■'r.+l r. 


] exp [~s.z. 


n .~s . 
--S 


for >0 (i == 1,2,..-^K). 


Differentiating log L(©, ,0 | 2 ( l) , zC 2 ), • » z(K) ) w-r-.to 
0^ (i = 1/ 2/ <» •»/!<) and a we get the following (K+1) likelihood 
ccruations ; 


(2-4-1) r^/[v 


i) 

^i+l 


~1 J “• ( ^2_”'^i ) ~ ® (^ ~ 1/2^»»<'/K) 


and 


^r^+1 

(2.4.2) S [ S 2^/^+Sizi^i, -riZ^^l,/£e -.l}]--d*eO, 


K ""r^i 


i=l j er j_+l ^ 


^ ^r®i ^ 


Equation (2-4-1) simplifies to 


®i = ^rt+3 ■*■ ^ (l-r^/n^) 

(2-4.3) = - a 

whore mj^ = log C rij_/( n^~rj|_) } (i = 1/2/«..*/K)« 

Substitution of (2-4-1) in equation (2-4-2) gives 

( 2 * 4 . 4 ) a = da“/d* = P^/d"^, 

where a* is given in equation (2-2*3)- 

A A 

The distributions of a and 0j_can bo easily obtained by 
applying ThoOrom 2*3*1 and its comllaries* 
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* Comparison of ML and LS estljna'tors * 

Epstein and Sobel (1954) showed that/ the LS estimators 
given in equation (2»1*3) are the minimum variance unbiased 
(MVU) estimators of 0 and a respectively# It follov7S that 
0* (i = 1/ 2/ »■ •»,!<) and a* given in equation (2»2«3) are the 
MVU estimators of 0^^ (i = 1/2/«»»/K) and cr respectively 
[see/Sarhan and Greenberg 196 2, p» 368 ]* 

Then from equations (2»2»^) and (2 #4 *4) we have 

(2.5-1) a* = d'''‘‘a/d , e* = -b 

Using the equations (2-2.5), (2-2.6) and (2-S-l), the 
moans, variances and mean square errors (MSE) of ML estimator 

A 

a are respectively 
( 2-5.2) E(a) = da/d* 

(2-5-3) Var(a) = daV<a*^ 

(2-5-4) MSE(a) = (d+K^)aV(d+K)^' 

From relations (2-2-5) and (2-5-3), it can be seen that 
Var(a) < Var(a*)- 

The relative efficiency E of a w-r-to a* is given by 
E = MSE(a*)/MSE(a) = (d+K) 

Note that, E is less than 1 whenever d > K/(K-2). Since in 
general d is large, E is less than 1 for most values of K > 3- 
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This shows that in general for K > 3^ o* is a better estimator 

A 

than a • Hov/ever, for K < 2, the value of E is greater than 1 

A it 

and the ML estimator cr even though biased/ is bettor than cr # 
This agrees with the conclusions drawn by Kambo (1973) for 
K = 1 » Further/ for fixed K, E tends to 1 as d tends to 
infinity. Table 2.5.1 gives E for K = 1(1)6(2)10 and 
d = 2 ( 2 ) 10 ( 5 ) 30 . This table shows that for K > 3 and moderate 
values of d/ the estimator is considerably better than O'. 

The equation (2.4*^) also gives 

E( 0 j^) = 0 i+C' ( bj|^~A and Var(e^) = { a^^+A^/d} o^'. 

Hence # 

(2.5.5) MSE(e^) = {(aj_4A|/d) + (b^-A^)^! 
where = dmj^/d’^ (i = 1,2/».*/K). 

A 

In the general case direct comparison of the MSB of 0^ 
and 0* is difficult. But from equations (2.5*5) and (2*2.5)/ 
it is clear that, for rj_ = 0 (i = I/2/.../K) MSE(0*) < MSE(0^) . 
For the two samples case, n^ =5/ ®1 “ ^2 “ 

r 2 = 12 and S 2 = 1/ the MSE of estimators are calculated from 
equations (2-2.5) and (2*5.5)/ and are given by 

MSE(e^) = 0.05 330^ < MSE(0^) = 0*8000^ 

MSE(e|) = 1 .43240^ > MSECe^) = 1 .36810^'. 


and 
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This shows that nothing can be concluded about the relative 
efficiencies of those estimator's* The proper estimator out 
of those two is the one with the smaller MSE* Similar 
conclusions are drawn by Kambo (19 78) for the case of one 
population* 


For later usO/ wo nov/ evaluate the variance of (8^-0 2 )* 


From equations (2*4*'^) and (2*4*4)/ we have 

S1-O2 = 

where m = ir^— m 2 « From Theorem 2*3*1 and equation (2*2*6) we 


get 

(2-5-6) Var(0j^-e2) = 


2*6* Estimators of the parameters under the hypothesis e^r=e2=0* 


The LS estimators and the ML estimators of the 
parameters 6 and a under the hypothesis 0^ = based on 

two independent typo II censored sanples arc discussed in this 
section • 

2 *6 *1 • LS cstliaators of the paramotors * In this subsection/ 
all unspecified noteitions are as given in Section 2*2 with K = 2, 
Denote E(x) = B T and VaHx) = d^D, 

ro oJ CNO ro 



1 C(l) 


e 

B =: 


and T = 


<X> 

1 c( 2) 

ro 

a 


fsjj ro 


— - 


where 
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Then the LS estimator of r is given by 


( 2 »6 »1 ) 


21^ 

<vO 


O 

o 


{ B * Q B )“^ ( B ' Q X ), 


CS 3 CS> CU 


ro ro 4 ro 


where ( B' Q B) 


ro (po rs> 


1/^^ + 1/^2 ^ 2/^2 

bf/ai + ^2^3 - 2 Q 


(2*6.2) = Jl . 

ro <vj CO _ ^ 


Q 


-( b ^/ a ^+ b 2/ a 2) l / a ^^+ l / a ^ 


where Q* = |B" Q B| = {(b^-b^)^ + d( a^+a2)}/a^a2. 
Further^ 


( B ' Qx ) = 


ro <0 <0 


( 1 ) ^ C 2) ' 

1 ■'' ^ r 2+]/ 2 


^ ( i) 


Simplification of equation ( 2 * 6 - 1 ) gives 


(1) (2) 
+1 ^r„+l 


( i) 


^ I. -I TJL i. qtj. ^9 ^ ^i^r .+1 

So +-af-> - <4 + srf - 


° - Q - 


i] 


and 


(1) ^(2) 

b. b- ^r- +1 ^r«+l 

(2.6.4) J* = A- [-(— + 4 - 1 

o Q,, a^ a2 a^ a2 


■} 


2 . 

Substituting for (P- + S 4.1/^i^ ^ equation ( 2 * 6 . 3 ) 

i=l ^i"*""^ 

by (2*6.4)/ and simplifying, v/e get the unbiased LS estimators of 


e and a as 
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(2.6.5) e*= 


anc3. 


( 2 «6 *6 ) 0^ — ^ ^ } /{ <3.( ajl^+a2 ^ 

!L 1L 

^ »-.i 9 

The variance covariance matrix of ^ is (B' Q b) ‘‘‘CT » 

«n> rs» <sj 

The equation (2.6.2) now gives 

(2.6.7) Var(eQ) = Qa^/Q% Var( 0 *) = ( a^+a 2 )a 
and 


Covar(0j/aQ) = ’~(^i^2'^2^1 


2.6 . 2 . ML estimators of the parameters . 

The likelihood function of X(l) and X( 2) is given by 

ro to 

2 n ! 

L(G,a |x(l) ,x( 2) ) = , g- ' /y [l-exp{ x^^lj^-ei/d] ^ 


r-i-Si 

. exp [- hsA^^K ~e)+ E )}]/X^^l- >e(i==l,2),a > 0 . 

^ ^ ^ I j=rj_+l ^ ^i"^^ 

( 1 ) ( 2 ) 

Without loss of generality, let otherwise, we 

( 2 ) ( 1 ) 

can simply relabel the samples- Now y = x^ ;+i-a, +1 i ° 

2 x> 

"5^ \ XT 

(2.6.8) L(0,a ix(l) ,x( 2) ) = Const [l-exp{~(x^^ .~0)/0l] 


€o ro 


-r^+l 


. [l-oxp{- 2 _ (34 )|j^- 9 Vo!] ^ exp 


( 1 ) 

for yiiri-, > e, 0 > 0, 
r^+1 - 
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(i) 


(l) 


where P = S [ .£ x. -s < +1 ] 


1=1 3 =r^+l 


1 1 


According to zero or non- zero values of and r^/ we have four 
different cases# 

Case (l): r- / r^j^O . The likelihood function ( 2 #6 #8 ) becoTnes,r 


L (0 , a 1 x( 1) , x( 2 ) ) = Const- exp [— P+f (x^^ 0 )}]/o’ 


3 * 


for 




(1) 


> 0, 0 > 0i 


A ( 1 ) A « 

It is clear that L is maximum for and 0^ = P/d - This 

has been also obtained by Epstein and Tsao (1953). 

Case ( II) : r^^ > 0/ r 2 = 0. Now, the likelihood function (2.6w8) 
simplifies to 


( 2 -6 .9 ) 


Hi-s^ 


LC0,ai z(l),z( 2 )) 

.o ,0 

= £ 2 ^ [1-0 S 

^^r^il > 0/ d > 0, 


z 


( i). 


i *"1 j=r.+l 
1 


■where z 


( i) fi±) 


(x. ~0)/0/ j = r^-i-l/ . . ./ n^-s . and i = 1,2. 


-j •— j - , - . ^ -X'-"' '1 1 

The maximizing equations for L(0,0 1 z(l ) , z( 2) ) arc 


(2-6-10) r^/[exp {z^^|^}-l] - f = O 
and 

(2.6.11) [exp{4l>^}-l] H- d* = 0. 

Solving the equations (2-6-10) and (2-6-11)/ we get 



37 


(2.6.12) e. = log 

No-to that . 

Case (III) ; =0, ^ 2 ^ 


function (2*6.8) roc3ucos to 


(1+r^/f), = P/d*. 

In this case the likelihood 


( 2 . 6 . 13 ) L(e,ff !x(l) # x( 2 ) ) = Const, { l-qw} 

ro CS5 

for 0<w<l, 0<q<l/ 

( 1 ) 

where q = oxp(-Y/Cf) and w = exp C*"(x^ -6)/G} * We maximize 

it w*r. to 0 first/ which is equivalent to maximising it w.r» to 
w* Unlike the Case (II)/ the maximum is not necessarily at 0 

{ 1 ) ( 1 ) 

point which is less than x^ , but could be at x^ also. 

Towards this end/ consider the function 

^2 f 

g(w) = (l~qw) w / 0 < w < 1/q# 

which is plotted in Figure 2*6.1. 



FIGURE 2.6-1. A function related to likelihood function- 
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Note that 1/q >1/ g(o) =0, g(l) = (l-q) and g(l/q) =0* 
Differentiating log g(w) w*r» to w and equating it to zero, 
we see that the ma:cimum of g(w) in the range 0 to 1/q occurs 
at where 

(2*6*14) Wq = f/{q(r2+f)} • 

Note that Wq may be greater than 1* 

Consequently, the differentiation of likelihood 
function w*r* to 0 gives the ML estimator, only if ^ 1 
w >1, then the likelihood function will attain the maximum 
w»r* to 0 at a point where Wq = 1/ that is at exp{ — ~6 )/cr} =» 1* 
Thus, the maximum in this case is at 0^ = thus have 

two possibilities. 

Case (i) : w^ < 1* With same notations as used in equation 
( 2 *6 -9 ) , we have 


L(0,a I z( l) , z( 2) ) 

* CO • CO 


= Const* 






+ E 

i j=rj^+l 



for ^ 0 > O* 

Diff erv-ntiating log L(©,a i z( 1 ), z( 2) ) wr* to 0 and o, we will 

CO 40 

get ML equations as 

- f = 0 . 

£* ^ 


( 2 .6 *15) 
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Solving these equations we obtain 

(2.6-16) 0^ = log (-^), = (P~fY)/d*. 

Case (ii) ^ Wq > 1 • In this case^ as shown above 0^ = ^ . 

Substihuting this in equation (2.6*1?), vre have to maximize 

(2.6.17) L(a,0 1 x(l ), x( 2) ) = Const, a ^ 

Q ro <%o 

(for y > 0, CT > 0, ) 

w.r. to O . For this, the maximizing equation is 

(2.6.18) d*a -!- r2y/{ e^'^*^~l} - P = 0* 

A 

The solution of this equation gives the ML estimator of 0 • 
This can be solved by using Newton~Raphson method with P/d* 
as an initial value. 

In Case (III), the procedure of choosing proper ML 
estimator is as follov/s : 

Calculate the quantities a = y/log (itry'f) S-nd 
b = (P~fy)/d*. Depending on the values of 'a' and 'b' we have 
throe possibilities. 

Case (l) ; a < b» In this case, the relevant JiL estimators arc 
given in Case (i). Since, the liTcclihood function attains its 
maximum v;ithin the pertinent range, < 1 . Equivalently 

> y/log (l+r^^f), 

whore = b* 
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C ase (2) s a > b* In this case, the ML estimators are given by 
6 q == and which is the solution of the equation ( 2*6 •is) ^ 

Note that, for 


f(a) = d'b + r 2 y/(e^^^~l) ~ P# we have 
f(0) =: < 0, 

f(a) = d^*'(a-'b) > 0 

and f'(a) = > 0, 

hence f?quation (2*6»18) has a unique solution in (0,a) • 

Case ( •^) t a = b i> In this boundary case, the ML estimcitors 

given in Case (i) and Case (ii) turn out to be same with 
* (1 ) A 

®o = 4 ' ^ 

Case (IV) ; r^ > 0^ r 2 > 0. Using the same notations as in 

ojqjrossion (2-6 #9), the likelihood function (2-6-8) can bo 
rewritten as 


-d* ( 1 ) 

(2-6-19) L(0,a 1 z(l) , a( 2)) = Const- o" [ l"exp{ ^ } ] 

<vj X. ^ "TX 


. [l-oxpl- 2 _ ^ uxpf- 


I - * 4jll> 


( 1 ) 

for z^ >■0/ > 0. 

r^+l - 


Consequently, 


( 1 ) 


log L(6,c7lz(l),z( 2)) = const.-d'*^' log a+r^^ log [l~e3q){~z^ 


+ rj log [l-o3!p{-y/o-4J'l.ii]-l’/a-£zi.^il 


.( 1 ) 
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The corresponding ML equations are given by 

and ^ ^ 

° a [ l-ej^C -z^^ a [ l-o3q){ ~ 2 ^^^ } ] ^ 

On simplification/ these oquations r^du-co to 

(2.6.20) r^/[exp{2^^^^}-l] + - f = 


p X*.! 4"1 

+ 4- ■'" =: 0 - 


( 2 .6 .21) 


r.z^2; 

^ ^2^^ .(1) P ^ n 

, ^ 4 ' '7 "' ■' '■■' f 2 ^ ^ 4 ^ 0 i 

GXp{ 2^^ +1 ? -1 GXp{ 2^2+1 ^ 1 


Eliminating [exp{ -1] in (2.6.21) by using (2«6.20), 


(2»6t22) 6 == log { 14r2y/{ P--d*a )} • 

Similarly/ eliminating [exp {zi^Lj-ll we obtain 

r2+l 


( 2.6 .23) e = ~ a 

ri+1 


log {l+rj^y/(d*a + fy~P)} • 


Equating (2*6,22) and (2*6.23)/ wo have 


(2.6.24) 


P-d*a 


]/[i +“r-^= — ]< 

d%+fy~P 


This equation can be rewritten in the follov/ing alternative 


forms : 
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(2.6.25) a = P/d*-r2y /[cl*{l+r^y/(d'^a+fy~P)}e^'^^~d*] 
or 

(2.6.26) a = (P--fy)/d* + r^^y/ [d*{ l+r^y/ ( P-d*a ) } ~ d'"'']- 

2 *6 *1 * Equation (2.6.24) has unique solution for o, and 
the solution lies in the interval ((P-fy)/d* P/d*), To this end, 
let 

f(a) = [n-r2y/(P'-d*a)]/[l4-rjy/(d*a+fy~P) ] . 

Note that, 

f(P/d*) = -oo and f ( (P-fy)/d'*^) =■- esq) { yd*/(P-.fy)} > 0. 

Hence, the equation f(a) =0 has atleast one root in this 
interval. Now, for concluding it has unique solution, it is 
sufficient to show that f(a) is monoton ically decreasing. 
Differentiation of f(a) v/.r. to a gives 

f'(o’) = ~exp (y/c7)y/a^ 

^ r^y J^oyd* r y -r.yd* 

r i — — — -n — — 1 r X *4* — 1- — — — 1 r ^ 1 

( d*a+fy-P ) ( P-d*a ( P-d*a ) "^ ^ ( d*cr+f y-P ) 

[l +J^]_y/ ( +f y-P ) ] ‘* 

< 0 for every o e ((P-fy)/d*, P/d*). 

Hence, f(a) is a monotonically decreasing function of cr in the 
above interval. 

Now, the ML estimator of 0 is given by one of the 
equations (2.6.22) or (2.6.23), and from expression (2*6.23), 



Q> 
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A (■ -, ) , . . 

it is easy to see that 8 < X since equgttion (2*6 *26; gives 

O “ 

^ > (P~fy)/d’‘'* The ML estimator of Q is the solution of the 
equation (2*6 *24) » It is difficult to show by second- derivative^ 
that these values really give the maximum of the likelihood 
function * But some numerical calculations/ show that the 

A 

likelihood function is really maximum at 0^ and 

Remark 2*6*2«» EJventhough r^ > 0/ r 2 0 in Case (W) , we can 
obtain some results for other cases frx)m this case by substituting 
and/or r 2 equal to zero* 

( i) If r^=:r2=0 then equation (2*6*25) gives = P/d*, as in 
Case ( l) • 

( ii) If rj^>0/r2=0 then aquation (2-6*25) gives cs^ =: P/d*, as i'- 
Case (ll)- 

( iii) If rj_=0/ r2>0 then equation (2*6*26) gives cr^ = (P~fy)/db 
as in equation (2*6 *16). 

( iv) If rj_=0, r2>0 and 0^ = then equation (2*6-24) reduces 

to aquation ( 2 *6 *18 ) • 


Re mark 2*6 *3 * A computer program for evaluating the ML 
estimators for all four cases is given in Appendix A * 


2 *6 *3 »• Comparison of the estimators ^ 

In general the explicit form of the I^IL estimators 
under the hypothesis is complicated* Hence the comparison becomes 
difficult* However, for ~ ^ the ML estimators are given 

in equation ( 2*1*12) » This can bo rewritten as 
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( 2.6 .27) 


e. 


min 




X. 


( 2 ). 


A 

O, 


= [< 30 ^ 


2 

- 2 
i=l 




O 


)J/d 


* 


Before deriving the MSE of the ML estimators, we first prove 
the following theorem. 

1 ) , „ ,^( 2) 

,-^,'^(1) ..(2) 


Theorem 2.6.1. Let U = ' H-n^X^ min( xj^ ^ ^ ^ 


and 


V 


min(X^-‘',X^‘'' ), 


where N = ’^^2/ 

( i) 2U/a has a distribution, 

( ii) 2N(V-0)/a has a X^ distribution, 

( iii) The random variables U and V are independently 
' distributed. 

( 1 ) (?) 

Proof . The jpdf of X^ and X^ ' is given by 


f( 


= nj^n2 exp n^x[^ ^+n2X^^ ^-NOi ]/cr 


for 4 ^ 74 ^^ > e, o > 


0 , 


T.,o-t- T, _ ^ +n 2 x :2 ^ ^“N min{ x^^ \ x^ ^ ^ } and v = min{ * 


This is not a one to one transformation •» Two inverse 
transformations are 

1) x^^ ^ = V, x^^*^ - V’^2 

and 

2) x^^'^ = v, = u/bj+vi 

whore u > 0 and v > 0 - The respective Jacobians of transforma- 
tions are 1/n^ and 1/n^ • The jpdf of U and V is then 
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1 

f(\3^v) =: exp [- ^ { n^v+u+n^v-N©} ] 


2 1 

+ —2 exp [~ g- £ n2V4-u+n^v-N0} ] for u > v > 0 

(2 -6 #28) r: { exp(->u/a)/cr} [n exp£ -N( v~0 )/a}/a ], u > O, v > e 
Then the marginal, pdf of U and V arc 
( 2 #6 #29 ) f(u)=:~e for u > 0 


and 

( 2 #6 -30) 


N( V--0 ) 


f ( v) 


N 

O 


for V > 6 


respectively- The equation (2-6 #29) shows that 2U/a has a 

distribution# Ebstein and Tsao (l953) also obtained the equation 

(2*6-29) by using a direct argument- Similarly^ 2N(V~0)/a has 
2 

a X 2 distribution- This can be derived by direct argument as 
well* Prom equations ( 2-6-28)/ ( 2-6 -29) and (2-6 -30)/ it follows 
that U and V are independently distributed* This completes the 
proof of the theorem* 

C oroll a.rv 2-6 -1 - '^ith and ©^ as defined in oqua.tion (2-6-27)/ 


A y 

( i) 2N(0^-0)/a has a distribution/ 


and 


( ii) 2d''a,VcT has a X^ distribution- 

2(d'''-l) 

Proof- The ML estimators given in equation (2*6*27) can be 
written as 

(2*6*31) Sq = V, 0^ = [2d0'Vc^ -i- 2U/0] 0/2d'*^ / 

where U and V are as defined in Theorem 2-6-1- Now/ part (i) 
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of the corollary follows from Theorem 2«6*1# By Corollary 2"3»1/ 
* 2 

2da /a has a X 2 ^ distribution • Since U is a function of 

( 1 ) ( 2 ) 

X. and X^ ■ , on applying a result similar to Theorem 
we see that a'^ and U are independently distributed# Consequently, 
26^0^0 has chi-square distribution vrith 2d+2 = 2(d'^-l)DF and 
the corollary follov/s# 

This corollary gives the following results ; 

E [2N(0^-0)/a] = 2 inplying -EC'Qq) = G+a/N, 

Var [2N(0^-0)/a] = 4 implying VarC©^) == 

E [2d*a jy/'a] = 2d*~2 implying = (d*-l)a/d* 

and Var [^d'h^/a] =s 4(d*-l) implying Var(a^) = ( d*-l )a^/d*^-. 

By these relations we have 

(2-6.32) MSE(0^) =20^/^'^ 

and 

(2-6.33) MSE(a^) = 

In equation (2-6-7) for r^^ = r 2 = 0, the variance of the LS 
estimators are 

( 2 -6 #34 ) Var( 0'J ) = d*a Vc d*( ^1+^2 ^ 
and 

(2-6-35) Var(0^^) = a Vf } <« 

A 

Eventhough is a biased estimator of d, from 
cou.ations (2-6-33) and (2-6-35) it follows that MSE(a )<MSE(a^)z 
whfirc a* is the unbiased estimator of a given by equation (2-6-6) 
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Further, as ojqpected MSE(a*) < MSE(a*) whore a * is given by 
equation (2»2#3). This can be seen from the equations (2«2»6) 
and (2*6 *3 5) and the fact that N^/Cn^+n^) < 2# 

Direct comparison of MSE(6^) and MSE(G ) seems to be 
difficult, but for - n 2 = n, the equations (2*6»32) and 
( 2 »6 »34) simplify to 

MSE( 0 ^) = a^y2n^ < MSECeJ) = d*(7V[ ('^*- 2 ) 2 n^ ] - 

Also for n^ =5, n2=15, s^ =0 and S 2 = 0 the equations 
( 2 *6 •32} and ( 2 *6 *34 ) give 

MSE(G^) = 0*005 > MSE(e*) = 0-00435 

This sh')ws that, nothing can bo said regarding the preference 
of the one estimator over the other arruong 0 and 0 • However, 
one can choose the estimator with the smailler MSE- 



table 2*5*1* The relative ef f icienqv " of a v;-r 






CHAPTER III 


TESTING OF HYPOTHESIS ABOUT LOCATION PaI^AiuETERS AGAINST 

ONE-SIDED ALTERNATIVES 

Introduction and test, statls'tics » 

In this chapter a test statistic from type II 
censored samplesis proposed to test the equality of location 
parameters of two exponential distributions against one~sided 
alternatives* The common scale parameter is assumed to be 
unknown* The null and the non~null distributions of the 
proposed test statistic are obtained* Some critical points 
and some values of power are tabulated* An approximation 
in terms of Student's t distribution for the null case is 
studied* 

Let be two independent samples 

from E(ej^/a) (i = 1,2)* Consider the problrim of testing 

Hq ; = 02 against on^.-sidod alternative hypothesis > ©2 

or H' : ^ 2 ' ^ unkno\«7n» 

For the right censored samplos (case r^^ = r 2 = 0) , 
Kumar and Patel (l97l) have proposed a test statistic (KP test) 
for testing against H 2 : 6^ ®2* Weinman (1973) 

oxtondod it for testing against one— sided alternatives* 
Against thair tost statistic is. 
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v7horG ( 7 ^ is the pooled ostlmator of a, given in equation (2*2-3) 
They had obtained the null and non— null distribution of as 


( 3 *1 « 1 ) 


and 

( 3*1 *2) 


P [Wo 5 


n^( l-n 2 c/d)~’ /n for c < 0 
l~n 2 ( l+nj^c/d)~ /n for c > 0 


P [Wq < dtp] 


~n^ e7?p(-n2 <¥*) (l'-n 2 c/d)'”^/n, c < 0 

n. (p 

Qd( 1 0)-n2G Q^£ | nj^c/d}/n 

+ e5q7(-n29)b^{ c^|n 2 c/d}/n/ c > C 


rcspe;ctivoly, where n = n^+n 2 / d = n^+n 2 ~r^"'r 2 -Sj-'S 2 ~ 2 ^ 


= d<^G^ cp = ^)/a, 

CO . 

Q_(xls) = 7 e ^ dy/(p— l)l, x > 0, p = 1/ 2/ -••/s > —1 

^ X ” 

IX) 

and L (xls) = / y^~^e"y'^®^Qy/(p-l ) I , x > 0^ p = 1, 2, - * * . 

^ o 

The expressions for the critical points obtained by 
Weinman (given in Section 1 *4) are actually interchanged, 

although his tabalated values are correct* By using the 
notations as in Section 1 *4, the correct critical points of W 



d 

n^ 


n, 


- 'TnV«,Tri“<n’ 





n. 


n^+n 




if n^(l-a) < n 2 d 
if nj^(l~(x) > ^2°'* 


are 
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>sJe. 




All the foregoing statistics are based on intuitive 
grounds# On the basis of LS and ML estimators/ v/c propose the 
following test statistics for testing against : 


'^LSE “ sECe^—e*) 


and 


'^MLE " (8 ^^-© 2 V®S'tii^3.tu of SE( 0^—0 2 )/ 

whi^'re 0* and denote the LS estimator and the ML estimator 
of (i = 1/2) respectively# Using oj^pressions for 0* and 6^ 

given in equations (2 #2 #3) and (2#4#3)/ and expressions for the 
variances of and given in equations (2#2#8) and 

(2 #5 #6)/ these statistics reduce tO/ 


O .1 .3) j, = [ C 

and 


i 

2/^i2 


X 


(3#1#4) T, _ rr^Xl) __v(2) 
where 0* is the LS estimator of O given in equation (2 #2*3)/ 


qi = bj_~b2, q2 = ^ ~ ^l’^'^2' 

rj^tl ^i"^^ 

b^ = S (n^~j+l)~^/ajL = S (n^-j+l)"^' and 
j =1 j =1 


m. 


log { n--/(n/-r^ )} for i = 1/2- 


1 ~ <• ' “i ^1' 

Both of those statistics against one-sided alternative arc 
equivalent to the statistic 

(3-1.5) T = {4^|i-4^|l}/U* - 

X « 
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Against very large values of T lead to the rejection of 
that is^ we reject against if T > where Cq, is 

determined so that p[t > c^lH^] = a»and a is the chosen level 
of significance of the test# Similarly/ against H', very small 
values of T lead to the rejection of We thus need the null 

distribution of T for finding the critical point c^^, which is 
derived in Section 3 •I# The non— null distribution is obtained 
in Section 3 #3# In Section 3*4/ the exact and approximated 
critical points are evaluated# Power function and its approxima- 
tion are studied in Section 3*5# 


3 *2 • Null distribution of the statistic T » 

We first prove the following lemma ; 

Lemma 3 #2 #1 # Suppose and ^2 are two independent random 
variates with pdf of as 


^1 


(3*2#l) f 2 _(z) ={l-exp(~z)} ^ e3(pC-(n^~r^) z}/B(rj^+l,nj^-r^) / z>0 

Then the pdf of Z = Z^-Z 2 is given by 


(3*2*2) g(z) 


2 j r„ (n 9 -r«+j)z 

H E . )b( r., +1/ f+j )e ‘ f z < 

j=0 J 


0 


1 ^ r. ~(n#”>r-, +j) z „ ^ n 

H E r-)B(r,+l/f+j)e ^ ^ / z > 0/ 

j=0 ^ ^ 


where H = { 1/B( rj^+l/n^-rji^)} / f s= 

B(p/q) = / u^^(l-u)^”*^ dU/ p > 0/ q > O# 
o 
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P roof >■ The jpdf of and is 




f(zj|^,Z 2 ) = H(l-e ) (l~e *) exp{~(n^-rj^)z^-(n 2 ~r 2 )z 2 ^ 


for 0 < 2 ^/ Z 2 

Making a transformation z = 2 j_~Z 2 and Z 2 = Z 2 ^ we see that the 
range of Z 2 and Z are Z 2 > inax(0, — z) and — <» < z < «» respectively* 
The inverse transformation is ^2 “ ^2* 

j acobian of the transformation is 1, and the jpdf of Z and Z 2 is 

X* 37 

(3«2»3) g( Z/ Z 2 ) = H [ l~exp{ -'( Z 2 + 2 )} ] ^ {l“exp(- 22 )} ^ 

- exp [~fz 2 ‘-( nj^~rj^) z ] , 22 ^ max(0,->z) /-«» < z < -o- 
The marginal pdf of Z is thus 


( 3 *2 *4) g( z) = 


/ g(z/22)<3.Z2 for z < O 

~2 


S g(z#22)<322 for 2 > 0. 


o 


For 2 <-0^ on putting y = Z+Z 2 / g( z) becomes 


~(n^-ri)z °° ~fz_ -(z-j+z) r- ~z„ r„ 

g(z) = He ^ ^ / e ^ [l~e ] ^ (l~e ^‘) “ dz. 


-z 




o 

r. 


JL r\ 

Now expanding by binomial expansion and simplifying, 

wo cot 


2 . r« (no-r_+j)z <» 


q(.)=HS (-Dihe = 2 

J=0 ^ 


; .,-Cf+J)Y (i_3y) 1 ay. 
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Substitxiting u = and simplifying, we obtain equation (?i2#2) 

--(z 2 +z) r- 

for 2 < 0» For z > G, we expand [ l~e ] of equation 

(3'»2»3) as a binomial sum and get the corresponding result of 
equation ( 3 .2 ;2) • 


Corollary 3«2»1 » If Y = Z + cp, then the pdf of Y is 


(3.2*5) 


g( v) = 


^2 _ (n2-r2+j )(y~ 

H E .^)B(r.+l,f+j)e for y < <P 

j =0 ^ " 

^1 j r. -(n^~r-+j)(y-<P) 

H E (-l)-'( /)B(r 2 +l/f+j )e ^ for y > <p , 

j=0 J ^ 


Proof . The corollary follows from Lemma 3.2*1-. 


Theorem 3*2*1 . The null distribution of T [where T is defined 
by equation (3*1*5)] is 


P [T<c|H^] 

' ^2 

H E (~1 )^'( .^)b( rj_+l, f+j )(l-h2( j )c)“’'V<i^2^ ^ ^ ° 

-> j=0 ^ 


"**1 X* 

l-H S (-•1)^'( .^)B(r 2 +l, f+j )(l+h^( j )c)“'^/dhj^( j ) for c > 0, 
j=0 ^ 


where h^ij) = (nj_-rj_+j )/d (i = 1,2), d = f--s^-S2~2, H and f are 
as in Lemma 3.2*1* 


Proof * Note that (i = 1/2) has the pdf given 

in equation (3.2*1) . Then under : 6^ =0 2 , we have T =: dZ/W, 
where Z = Z 3 _-Z 2 and W = da*/a* By Theorem 2.3*1, Z and W are 
independent with pdf of Z given in Lemma 3*2*1, and that of W 
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given in equation (2-^*6)- Hence, the jpdf of Z and W is 


g( Z/w) = 


1 ? 

^ 1 (n^-rotj)z j 1 

H ) B( +1 , f j ) G ' ' e \; ~ /( d-'l ) I , vj>0. 


j=0 


z<0 




rs J ^ 


j=0 


-w d-1 ) 1 , W>0, z>0 . 


Make a transformation t = dz/w, w = w« The inverse transformation 
is z = tw/d, w = w, and the Jacobian of transformation is 


1 9^ Zr w) 
a(t,w} 


The jpdf of T and W is 


w/d. 


g(t, w) 


2 , r -£l-h„( j)t}w , 

H E (-!)-'( /)E(r.+l,f+j)e ^ w /d | , w > 0, t < 0 

j=0 J 

,• t- -{l+h^(j)t}w ^ 

H s (-D-^C , )B(r,+l,f+j)e w /dl, w > 0, t > 0. 

j =0 


Integrating out w, wc got the marginal pdf of T as 


(3*2.7) 


f(t) 


" IC 

H S (-1)^( .^')B(r^-M,f+j){l-h2( t < 0 

j =0 ^ 

H E (~1)^'( /)B(r^+l,f+j){l+h.( j)t}~^^'‘'^\t > 0. 
j =0 J - J- 


The desired cdf of I now follows from equation (3*2*7)* 
Note that, for r^^ = r 2 = 0, the equation (3*2*6) reduces to 
equation (3*1*1) given by Weinman ^* (1973)- 
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^ 2 ( 0 ] <p) 


^ i "^1 

Q^(c 3 _| 0 )-K_ 2 ^(- 1 )'’( ^■^)B(r 2 +l,f+j) »exp{ ( j )d<p} 


^2 . r 

Q^{ CjL 1 ( j ) c}/dh^ ( j ) +K 2 (-1 ) ■^ ( ) B( r^+1/ £+j ), ' 


. ej<p[~h2( j)d<pj *L^{Cj_|h2( j )c}/dh2(j ) 


and the remaining notations are given in equation (3*1 *2) and 
in Lemma 3 -2 #1 # 


Proof . Let = 1/2). Then fol3-ov;s the 

distribution given in equation (3.2*1)/ and Y = Z^-Z 2 +<P 


^ ^r^+l''^r*^'+l^/^ has the pdf given in Corollary 3*2.1, Cl earl 


,(1) ,,(2) 

T can be written as dY/VJ® Using similar results as in the null 
case, we now have the jpdf of Y and W as 


where 
•*” ( Y f ^) 

and 

f 2(y/w) 


f ( Y/ w) 


fj_(y/w) for w > O/ y j 5 V 
f 2 (y/'w) for w > O, y > <p. 


XT 

2 -• r„ ( no-r, 5 +j ) (y~<p) ^ 1 

I-I S (->1)^( .^)B(r-+l/f+j)e ^ rw.d-1 

j=0 J 


-w^d-i/ ( ) j 


1 ■ r. -r-, +j ) (y- <p) „ ^ 3 ^-, 

H S (~1)^( .^)B(r^+l/f+j)e ^ ^ eV yi 

j=0 J ^ 


d-1) i. 


Making the transformation t = dy/w and w = W/ the corresponding 
inverse transformation is y = wt/d, w = w. The region w > 0, 
y < <p is transformed to w > 0, wt < d<p. Similarly/ the region 
w >■ 0/ y > (p is transformed to w > 0/ wt > d^. Note that /the 
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boundary y = is transformed to the hyperbola wt » d(P# (see. 
Figure '^-3*1)* The Jacobian of this transformation is w/d, and 
uhe jpdf of T and W is 

g^(t/w) for w > 0, wt < d<P 


g( t/ w) 


for w > 0, vt > d'<p f 


where 


and 


^2 . r h (j)(wt-d(p) , 

g^(t,w) = H 2 . ) B( rjL+1, f+j ) e “ " 

j -0 ^ 


e w /d ! 


^ r. -h. ( j ) ( wt~d «P) . 

g^Ct^w) = H E .^)B(r2+l/f+j )e 

j =£) ^ 


“*w a,, , 
e w /d I • 



FIGURE Joint pdf of (T, V7) for 0» 

The cdf of T upto the point c is the integral of joint 

density of (T^w) over the shaded region als shovm in Figure 3i«3*l 
This is given by 

oo c 

p[t<c1<p]=/ (/ gj_dt) dw = (say) for c < 0 


and 

00 0 “ dfp/\f 

P It < cl <P j = / (/ Pi dt)dw + / (/ 

““ o —00 d9>/ c o 


dt) dw 
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d<P/c C oo c 

+ -f (/ gj dt:)dw 4 - / (f d-t)dw 

o o ' d<p/c d(p/w 


o o 


(s 3 iy) for c > 0- 


Simpli float ion of I.'s (i = l/ 2 /’?, 4 / 5 ) Involves lengthy esq^ressions. 


altho\igh the method is straight forward. Simplification of one 


factor I- is illustrated below : 


1 c 


o o 


= H S ,^)B{r,+l,f+j) 


-h 2 (j>a<p Cj^ h,(j)wc -W,a -1 


/ ( e 


1 ) 


^2 r ~h2(j)d<}) 


^1 d~l ^ ^ 

r f w ■ e 
I s 71CTTI 


^ -w , - 

^ rdZTTT ^ 


^2 r -h2(j)d<J) 

3 ^ 


n [Ld(cilh 2 (j)c)-l+Q^(c^l 0 )] - 


Similarly/ 


L = H E (~l)^( . )B(r-,+l/f+j)e 
^ j =0 J 


. N { l~h„( j )c} "*‘V<ah.( j ) / 


I = H i (~ 1 )^( /')B(ri +1 / f +j ) G 
^ j =0 ^ 


'h 2 ^j)d^/dh^( j)/ 


I3 = H 


* ^ X* C 3 ) (p 

i (-l)^( /‘)B(ri+l/f+j) {l-e ^ }Q^(c^lO)/dh2( j) 


and 
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= H 


1 . r 

j=0 J 


B(r2+l,f+j) 

“h^Tjl [Q<i(c^|0) 


-- exp {h^(j)d<P} Q^{ Cj_ih^( j)c} ] . 


Combining these expressions and using Corollary 3 *2 #2/ the cdf 
of T as given in equation (3.3*3) is obtained. 

Bor calculation purposes/ we need some further 
simplifications. For this, note that 


QjjC x| s) 
and 


L^(xls) = 


2 exp{-x(l+s)}{x(l+s)} i 1 +s) }/ s > ~1 
j=o 


X /d I for s = 1 

r n 

[l~ 2 exp£-x(l--s)}{x(l~s)} j iy(l~s) for s /= 

j=0 


The proof of these results and a computer program for evaluating 
these functions ere given in Appendix B. 


Substituting these results in expression (3*3*1) we 
obtain the cdf of T as 


(3.3.2) P [t < c|(p] = Pj_(cifP) = 


Gj^(cicp) 
G^(cl <f)) 


where 


for c < 0 


for c > O, 


^2 . r„ B(r-+l,f+j) Gxp {-h,(j)d<P} 

G/cl(p)=H 2(~1)^(/) i 1 

j=0 J dh2(j){l-h2(j)c} 
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and 


d-l "'^1 . 

G^Cclcp) = 2 cJ--H E 

i =0 1 


j =0 


-TF^nr- 


d-~l exp(-G. ) ^2 

S -1-.- 1 -- - _ ^ H E 

i =0 i!{ l+hj^( j j =0 


^ dh^( j ){l~h«( j )c}'^ 


*^c 

. [ exp{-.h 2 (j)d(p} ~ e ^ 2 i4 { l~h.( j )c} ^ ] 

i=0 ‘ ^ 

for h 2 (j)c 1 (j = 0, 2/ • • •/ r 2 )' However/ for h 2 Cj)c == 1 , 
the ( j •fl) 'bh term of the last factor in 02 ( 0 ) <p) aixnply becomes 


H(-l)^' ( j^) B(r^+l/f+j) exp {~h 2 ( j)d(j)}c^/{ dh 2 ( j )dl } 


Note that, for cp = 0 , equations (3.-^*l) and (-^^3-2) reduce to 
the null distribution of 1., given by equation (3*2«6)« It is 
easy to show that, for r^ = r 2 = 0 , the equation (3*3 #2) reduces 
to equation (3*1 *2), the non~null cdf of the test statistic 
considered hy Weinman et al . (l973)» 


Usually, the null hypothesis is also taken as one-sided 
against ono-sided alternatives* Thus, if we take H* : 0 . < ©o 
and tost it against : 0 ^ > Q then we also require the power 
function for nt^gative values of <p • The following corollary 
gives the necossary distribution theory results* 


Corollary 3*3*1 * The non- null cdf of T for 

M^( c| cp) 

( 3 . 3 . 3 ) P [t < cicp] = P 2 (ci<ji) = 


M2( cliip ) 


9 < 0 is 
c < 0 

c > 0 , 



where 


6 2 


^2 . r 

M^Ccifp) = 1-Q^(c^|0)+H 2 (-1)^( ^?)B(r 3 ^+l^f+j) expC-hgC j)d(P} 

^ 

.Q^{cj_|-h2( j)c}/dh2( jO-H E (-l)^( .^) B( r2+l/ f+j ) 

j =0 ^ 

. exp{ ( 3 )d(P} L^£ I ( j ) c} /dh^^ ( j ) , 
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This is given by 


P [t < cicp] = 


oo 1 o 

!--[/ (f g2<3-t)dw + / (/ g 2 dt.) dw 
o o o c 

oo o 0 ° dcp/w 

+/ (J g2dt)dw + / (; gj^dt)dw] 


c 


C 


OO 00 

1 -X (/ g2dt:)dw for c > O. 

o c 


for c < 0 


On simplification. 


this gives the equation 


C 3 *3 *3 ) • 


Romarh 3 «3 »1 * Ai „'iOugh due to the complicated expressions 
involved, it is not possible to show that the tost based on T 
is unbiased, yet we feel that this test is unbiased# Limited 
calculations carried out in Section 3»5 strengthen this feeling* 


Remark 3 *'^'» 2 » Critical points and the power values for testing 
Ho : 05_ =02 against < ©2 can be obtained by relabelling 

the samples* 

3 * 4 # Exact and approximate critical points o£_.T « 

3,4*1. Exact critical points # The critical point is 
obtained by solving the equation 


(3.4*1) P [t > Cq^iH^] = a, 

whore a is the chosen level of significance of the test. From 
equation ( 3 . 2 . 6 ) it is cl oar that c^(< O) is the solution of 

^2 . r 

1~H S (-^ 1 )^( .^) B(r^+l,f+j){l-h2( j)cQ^}'' /dh2(j) = a, 
j=o ^ 

if > 3 ~a, where 
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2 , r 

P It < OiHq ] = h £ .^)B(r-+i,f+j)/dh„( j), 

j =0 J ■*■ & 

Similarly, if < i-a, then Cq^(> O) is the solution of 

Note that, if = 0,then 

Pq = l-BCr^+l/fVBCrp+l^np-r^) ♦ 

Now, if < l~a, then is given by 

Cq^ = d [{ (l-PQ)/a}^/'^ ~ l]/n^ > 0. 

Similarly, if r^ = 0 and 

Pq = B(r]_+1, f )/B( r^+1, ) > l-a, 

then is given by 

= d [l~{PQ/(l~a)}^/'^]/n2 

In particular, if r^ = r^ = 0, then P^ = n^/Cn^+n^) -and 
■ ^ ■ '■rn-+np7 l-a)-> U“ 2 

~ S HjCI-cc)* 

other than these cases, can be obtained by evaltiating P^ 
and then applying Newton ~ Raphson method to the pertinent 
equation* Note that if r^^ < r^ then it may be more convenient 
to calculate P [t > OiH^]* Some critical points are tabulated 
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in Tablo 3 »4 •! and Table 3 •4*2 for a = 0»05* These are tabulat-sd 
for following combinations of saitrole sizes and censoring patterns?, 

(i) for n^=n 2 =l 0 ,(r 3 _/r 2 )=( 0 , 0 )^( 0 ,l),( 0 ,, 2 ),(l, 0 ),(l,l),( 2 ,, 0 ) 
and d = 4(2)16 in Table 3*4.1^ and 
( ii) for ri^=10,n2=6(2)24^(r]_,r2)=(0,0),(0,l)^(0,2),(l^0),(l,l), 
(2,0) and d = 12 in Table 3 #4 *2. 


Note that and S 2 appear only through d* Due to many variables 
involved, it is not possible to give a large table of critical 
points* In those tables, we have provided mostly those values 
which were later used for studying the power function of the test- 

initial value for solving equation (3*4*1) can bo ta}ccn 
as an approximate critical point, given in Section 3-4*2* 


3*4*2. Student's t approximation* The asymptotic distrlbut J.on of 


(1) (?) 

sample quantiles and X^ is well known (for example, see 

■'rX XT 2 “tX 


David 1931, o* 255)* We however use the exact mean and variance 


of Z = X^^^, .-xi.^^, . under Prom the results given in Section 2*2, 


we have 


(3*4*2) E(Z) = Ba, Var(z) =Aa^, 

r,+l 

where B = bj^-b 2 , A = a^4a.2, ~ S (n^^-j+l) , 

J “1 

Ui+l 

b. = t (n.-j+l)”^ (i = 1^2)- 

«L. # A 

j=i 

Using the asymptotic normality we see that 
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= AN(0,i). 

following Tiku (l98l)/ we can approximate the null 
distribution of 

(3. 4 .3) " (T~B)/VA 

as Student's t distribution with d DF. Therefore « the exact 
critical value c^. is approximately given by 

( 3 *4 .4) c^ = VA + B = , 

where c^ is the upper o cth percentile point of Student s t 
distribution with d DF. 

3 .4 ,3 ♦ N ormal approximation » We now consider a normal 
approximation for the null distribution. For this we need the 
mean and variance of T* 

*■ 2 

From Corollary W = 2da /a has a ^2d* 

(3»4*5) E [l/VJ^j = (d-m~l)l/[2™(d-l) t m = 1, 2, , - (d-l) . 

Now/ T can be written as 

(3 •4*6) T = 2dZ/(aW). 

From the independence of Z and W, and using equation (3-4 *2)^ 
we get 

( 3 . 4 . 7 ) E(T) =d3/(d~l), Var(T) = d^{(d-.l)A+B2}/{(d-.l)2(d~2)} . 
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VJ-s now use the fact that 

Var(T)} V2 d 

where E(t) and Var(T) are given in equation (3»4*7)« Hencey 


(3.A^) Gq, = Gj^‘*{Var(T)} + E(T) = c*, 

where is the upper ath percentile point of noirmal 

d i st r ibut ion « 


The exact values c^ and the approximated values c* and c 
obtained from the equations (3*4 • 4 ) and (3 -4 *8) are tabulated 
in Table 3*4*3, for a = 0*05 and for some selected values of 
nij-n^/r^/r^ and d. On the basis of this and some other 
calculations, it is observed that, in general nomal approximation 
valuG C 9 is bettor than for r^^ > r 2 / otherwise c^ is bettor 
than cl^* 


3 ^5 • Power function and its approximation* 


The power of the test T for testing H* ; < ©2 

against : 0^^ > 02 is givean by 

1~P^( Cq^I (Jj) for (p > 0 
1-P2(cq,|^) for y < O, 


( 3 .5 .1 ) P [T > c^ 1 (p ] = 


whore c^ is the exact critical point, P^^CCq^I^) and P 2 (Cq.I <p) 
aro given in equations (3*3*2) and (3*3*3) respectively* For 
a = 0*05, the power of the test T is tabulated for various 
values of ^ in Table 3*5*1, Table 3*5*2 and Table 3*5*3 for 


SO^ 
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the following combinations of sample sizes and censoring 
patterns : 

(i) for n3_=n2=lO,d=l6 and ( r^, r 2 )=( 0 , 0) , (o, 1 ), (0^ 2), ( 1 , 0 ) , ( 1 , 1 ) , 
(2,0) in Table 3.5 .1, 

(ii) for n 3 _=n 2 = 10 ,d= 4 ( 2 )l 6 and r^=r 2 =:l in Table 3-5.2, and 
( iii) for nji^ =10, n 2 = 6 ( 2 ) 18 and rj^=r 2 =l in Table 3.5.3. 

The vax'iation in power due to different combination 
^1 ^2 showed in the Figure 3.5.1* From this and some 

other calculations the following points emerge. 

(a) Tablo 3*5.1 and Figure 3*5.1 show that the test T is more 
sensitive for variatioi in r^^ compared to variation in 
^2* Consequently, it is desirable to take greater care 
in handling first sample so that censoring on the left is 
reduced to a minimum for this sample. 

(b) It is clear from the Table 3.5-2 that fot fixed n^/n 2 , r^ 
and r 2 the power of the test is relatively loss affected 
by variation in right truncation which is represented by 
variations in d« 

( c) Table 3 . 5.3 shows that for fixed r^,r 2 and d, the power 

of test is not much affected by increasing the sample size* 


Since the power function given in equation (3.5*1) is 


very complicated, for large d and moderate critical points (c ) 


the following normal approximation is suggested. 

Let U = 

^I’^i r2+l u 
normal approximation of chi square distribution yields 


Note that, for large d, a 
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^ 2l 

a = AN(a,a /d) and for modorate value of c the effect of c„a’'^ 

vX (X 

is not negligible in the linear combination defining U-«, Then 
from Tn«’.oi:em and Corollary 2-3-1/ we see that under 

E(u) = (<?)+B-c^)a and Var(u) = ( A+c^J/d)cJ 
whore cp = (e^~e 2 )/a. Therefero, 

P [t > Cq,i(p] = P > c^Kp] 

!L 2 

= P [u > OlC,o] 

(3-5.2) ^ 1 ~ §(h^), 

ixi 2 

where b* = -(<p+B--CQ.)/(A+c^/d)^'^^ and §(x) = S dy/'^ 1 2^ ) , 

— OO 

The exact power [left hand side of equation (3-5-2) J and the 
approximated power [right hand side of equation (3-5-2)} are 
tabulated in Table 3-5-4 for a = 0-05- As can be seen from this 
table, the normal auproximatioin is quite satisfactory in this 
case* However, some calculations performed for small values of 
c^ show that this approximation is not that good for such 


situations 
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TABLE -^,4 

•1 -« Exact upper 

critical 

points 

of the test 


statistic T 

for a = d 

-05 and n^^ 

=nj=10. 









\( , r2) 

d 




(1,0) 

(1/1) 

(2,0) 

(0,0) 

(0,1)^ 

(0.2) 

4 

0-3113 

0.1901 

0.0818 

0.6401 

0.4730 

1 *0100 

6 

0-2807 

0.1776 

0 .079 3 

0.5575 

0.4280 

0.8614 

8 

0.2668 

0.1717 

0.0780 

0.5209 

0-4054 

0 -7960 

10 

0 -2589 

0.168 3 

0 .077 3 

0-5002 

0 . 39 2 5 

0.759 3 

12 

0 .2538 

0.1 661 

0 .0768 

0.4370 

0 .3842 

0.7358 

14 

0-2503 

0-1645 

0.0765 

0.4778 

0.3784 

0.7195 

16 

0.2477 

0.1633 

0 .076 2 

0.4711 

0.3742 

0 .7075 

TABLiE 3 ^4 

•2 • Exact upper 

critical 

points Cq^ 

of the test 


Sto 

.tistic T 

for (X = 0m05f n^=l0 and d =: 

12. 








^2 

(0,0) 

(0,1) 

(0,2) 

(1,0) 

(1,1) 


6 

0-2194 

0.095 2 

-1 -4357 

0.4436 

0 .29 28 

0.6863 

8 

0.2396 

0*1^70 

0.0321 

0.469 3 

0 .3474 

0.7154 

10 

0 -2538 

0.1661 

0 .0763 

0 .43 70 

0 .3842 

0.7358 

12 

0-2644 

0.1876 

0.1097 

0 .5000 

0 .4110 

0.7506 

14 

0.2726 

0.204^ 

0 .1 351 

0.5099 

0.4315 

0.7618 

16 

0 .2792 

0.2176 

0 #1553 

0.5178 

0.4476 

0.7707 

13 

0.2846 

0.2285 

0.1718 

0*5 242 

0 .4606 

0*777? 

20 

0.2891 

0.2^75 

0 .1856 

0.5295 

0.4714 

0 .78 36 

22 

0 *29 29 

0.245 3 

0*1972 

0.5339 

0.4805 

0.783 5 

24 

0.296 2 

0.2519 

0.2072 

0.5377 

0 .488 2 

0.79 27 
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TABLE 3 •4.3* Comparison of exact (%), Student's t approxlmatior i 

( c^') and normal approximation critical 

points of the test T for cc =:0<»05» 


^1 


0 

0 

0 

1 

1 

2 


0 

1 

2 

0 

1 

0 


n^^n 

2-10, d = 

13 

n^=:n2=15 

, d = 22' 


% 

c* 

1 

^2 


c* 

c^ 

2 

,2519 

.2505 

.26^2 

• 1618 

-1620 

.1665 

.165 2 

.2074 

.2197 

.1038 

.1318 

.1359 

.0766 

.1518 

.1712 

1 .0543 

.0941 

.1002 

.48 20 

.4296 

.4604 ! 

-2939 

.2746 

*2855 

.3811 

*3744 

.3935 

.2393 

.2374 

*2440 

.7 270 

.6 240 

.6827 

.4358 

.3908 

.4111 


TABLE 3«5*1- Power of the test T for testing against 
for cc. s 0 «c 557 n^=n 2==10 2 in<i d _ 



(0,0) 

(0, l ) 

(0,2) 

U ,0) 

(1-1) 

(2,0) 

~0 -50 

.0003 

.0003 

-0003 

-0009 

.0009 

.0017 

-0-40 

.0009 

,0009 

.0009 

,00 20 

-0020 

-0035 

-0 -30 

.0025 

-0025 

.0025 

.0046 

.0046 

.0071 

~0 -20 

.0063 

.0068 

.0068 

-01 04 

.0105 

.0139 

-0 .10 

.0134 

.0184 

.0184 

.0231 

,0232 

.0267 

-0 -05 

.0303 

-0303 

.0^03 

-0341 

.0342 

-0367 

0-00 

.0500 

.0500 

-0500 

-0500 

*0500 

.0500 

0.05 

.08 24 

.08 24 

•08 24 

-0726 

.0725 

-0675 

0-10 

-1358 

-1358 

.1355 

-1042 

.1038 

-09 03 

0-15 

-2 231 

-2218 

.2146 

.1476 

.1467 

«1 19 4 

0.20 

•3540 

*3416 

-31 36 

.2054 

*2035 

-1559 

0.26 

-5137 

.4758 

-4205 

.279 2 

.2754 

.2006 

0-30 

-666 2 

-6016 

.5249 

-3683 

.3613 

.2539 

0-35 

-7851 

-7070 

-6198 

-4684 

-4566 

-3155 

0 -40 

-8665 

.7895 

.1020 

-57 21 

.5543 

.3843 

0 -45 

-918 3 

.8514 

.7703 

^6709 

.6472 

-4581 

0.50 

.9504 

•3966 

.8 256 

.7577 

•129 A 

.5341 

0 .60 

J 9817 

.9514 

.90^ 

-3329 

-3523 

.6803 

0.70 

.9933 

.9778 

.9479 

.9498 

.9 259 

.3016 

0-30 

-9975 

.9901 

-9723 

.9301 

.9646 

.3331 

0.90 

.9991 

.9956 

-9362 

.9924 

-9836 

.94 21 

1 .00 

.9997 

.9931 

-9931 

c 9972 

-9926 

.9 722 
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TABLE 3 *5 *2 # Power* of 

the test T 

for testing 

again st 



for 

a = 0 

.05, n^=n2 

=10 and 

^1=^2 = 

1 ."■ 










<P 


6 

3 

10 

12 

14 

16 

o 

in 

• 

0 

1 

-0010 

.0009 

.0009 

.0009 

.0009 

-0009 

■.'00<^ 

-0 . 'iO 

.002 2 

.0021 

.00 21 

.00 21 

.00 21 

.0020 

.0020 

-o,so 

.0050 

.0043 

.0043 

.0047 

.0047 

<^047 

.0046 j 

-0*20 

.0110 

.0103 

.0107 

.0106 

.0105 

.0105 

j 

.0105 ! 

-0 .10 

.0239 

.0236 

.0234 

.0233 

.0233 

.02 32 

.02^2 

~0 .05 

.0^47 

.0<t45 

.0344 

.034^ 

.0343 

.0342 

.0342 j 

0.00 

.0500 

.0500 

.0500 

.0500 

.0500 

.0500 

.0500 ' 

0 .05 

J311Q 

.0716 

.0720 

.0722 

.0723 

*0724 

.0725 

0 .10 

-099 2 

.101 2 

.10 22 

.1029 

.1033 

4037 

-1033 

0 .15 

.1357 

.1403 

.1423 

.1444 

.1454 

.1462 

-1467 1 

0.20 

-1303 

.1901 

.1952 

.1935 

.2007 

•2023 

.2035 ! 

1 

i 

0 ^25 

.2341 

.2504 

.259 7 

.2653 

.2700 

.2732 

.2754 1 

0.30 

.2942 

.3193 

.3343 

.3449 

-3520 

• 3573 

-3613 i 

0.35 

.359 2 

.3955 

.4173 

.4320 

.4425 

.4505 

*4566 

o 

o 

-4266 

.4740 

.5026 

.5219 

.5353 

.5463 

-5543 

0 .45 

.4942 

.5519 

.5363 

.609 3 

.6256 

.6370 

.6472 

0.50 

.5600 

.6260 

.6644 

.6394 

.7069 

-7197 

.729^ 

0 .60 

.6795 

.7539 

.79 37 

.3173 

.0 337 

.3440 

-0528 

0 .70 

.7772 

.343 7 

.30 23 

.9016 

-9131 

.9207 

.9259 

0.30 

.3512 

.9123 

.9376 

.9501 

.9572 

.9616 

.9646 

0.90 

.9041 

.9517 

.9635 

.9759 

.9793 

-93 2i 

.9336 

1 .00 

.9401 

.9745 

.9347 

.93 37 

-9907 

.9919 

.99 26 
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TABLE 3*5 

• 3* Power 

of the 

test T for testing 

against 

Hj, 

for a 

= 0.05 

y = 10, 

d = 12 

and r^=r ‘2 f 

1 . 








6 

8 

10 

12 

14 

16 

■nhhbi 







-0 -50 

.0009 

,0009 

.0009 

,0009 

-0009 

.0009 

O 

0 

1 

.00 21 

.00 21 

.0021 

.0020 

.0020 

.0020 

-0 .30 

.0048 

.0047 

.0047 

.0047 

.0047 

.0046 

~0 .20 

.0107 

.0106 

.0105 

.0105 

.0105 

.0105 

O 

• 

o 

! 

.0235 

.0233 

.0233 

-0232. 

.0232 

.0232 

-0 .05 

.0344 

.0343 

.0343 

.0342 

.0342 

.0342 

0 .00 

.0500 

.0500 

.0500 

.0500 

.0500 

.0500 

0 .05 

.0718 

.07 21 

.0723 

.0724 

.0724 

.0725 

0.10 

.1018 

.1028 

•1033 

.1036 

.1037 

.1038 

0.15 

.1418 

.1442 

.1454 

-1461 

•1464 

.1466 

0 .20 

.19 29 

-1981 

.2007 

.2021 

.2028 

.2033 

0 .25 

.2550 

.2650 

.2700 

.2727 

.2742 

.2751 

0.30 

•3 256 

.3430 

-3520 

.3568 

.3595 

.3612 

0.35 

•4010 

.428 3 

.4425 

-4504 

-4548 

.4576 

0 .40 

.4771 

.5153 

-5358 

.5472 

.5538 

.5578 

0.45 

.5502 

.5990 

.6256 

.6405 

.649 2 

.6544 

0 .50 

.6179 

.6754 

.7069 

.7247 

.7349 

.7411 

0 .60 

.7323 

.798 3 

-8337 

.8535 

.8647 

.8713 

O 

o 

.8179 

.8811 

.9131 

-9 301 

.9394 

.9447 

0.00 

.8787 

.9325 

.9572 

-969 3 

.9756 

.9789 

0 .90 

.9 205 

.9627 

.9798 

•9874 

.9910 

.99 28 

1 .00 

.9485 

.9797 

.9907 

-9950 

.9968 

.9977 
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TABLE 3.5 *4 « Exact and approximated power of the -test Tj. 
for a = 0.05,n^=n2=l0 and d“ = 16- 

0-076 2 0*^742 0-7075 

(r^/r 2 ) ( 0 / 2 ) (l,l) ( 2 , 0 ) 


Exact 

i 

Approx- 

Exact 

ApDrox - 

Exact 

Approx - 

-000 ^ 

-0001 

-0009 

•0001 

-0017 

.0003 

-0009 

.0006 

-0020 

-0004 

.0035 

.0010 

•0025 

-0027 

-0046 

.0018 

•0071 

.0031 

•0068 

.0099 

.0105 

.0065 

.0139 

.0085 

.0184 

•0304 

•0232 

.0201 

•0261 

.0212 

-0303 

•0497 

-0342 

•0333 

-0367 

.0320 

A3500 

•0777 

•0500 

•0528 

.0500 

.0470 

AD8 24 

•1164 

.0725 

#08 04 

.0675 

.0672 

•1355 

.1671 

-1038 

.1178 

•0903 

.09 35 

• 2146 

•2302 

.1467 

•1661 

-1194 

.1 268 

•31 36 

-3047 

•2035 

•2256 

*1559 

.1675 

•4205 

•3884 

.2754 

•2955 

.2006 

.2158 

-b 249 

•4777 

» 36 1 3 

• 3741 

.2539 

.2713 

-6198 

.5681 

-4566 

•458 3 

*3155 

•3332 

-7020 

•6550 

•5543 

.5444 

-3843 

.3999 

.7703 

.7344 

.6472 

.6285 

.4581 

•469 7 

•8 256 

•8034 

•7294 

.7068 

•5341 

.5405 

-9029 

.9047 

•8528 

-8357 

.6803 

.6761 

•9479 

.9611 

.9259 

.9 206 

^016 

.7916 

.9728 

•9867 

.9646 

.9673 

.8881 

.8784 

•986 2 

•9962 

.9836 

.9885 

.9421 

.9 361 

•99 31 

.9991 

.99 26 

.9966 

.9722 

.9698 j 






CHAPTER IV 


TESTS OF HYPOTHESIS FOR LOCATION PARAMETERS 
AGAINST TWO-SIDED ALTERNATIVE 


4 . 1 * IntrorJuction and test st.atist.ics « 


In this chap cor, test statistics based on the LS and 
ML estimators are proposed for testing H^ : 0^ = © 2 / against 
the two~sided alternative H 2 ! ^ ® 2 * null and non-null 

distributions of the proposed test statistics are derived in 
Section 4.2, In Section 4.3, the LR test statistic for testing 
H against H is discussed. An approximation to the critical 
point is investigated in Section 4*4. In Section 4.5, the 
performance of tests is studied and comparisons with Tiku's 
(l98l) test and LR test are made* 


Epstein and Tsao (l953) discussed the LR test based 
on right censored samples. For same problem, Kumar and Patel 
(1971), also proposed a test (KP test) based on V cr I, 

where 0 * is the pooled estimator of a, given by equation (2.2-3). 
They obtained the null distribution and have tabulated some 
critical points of • Duboy (1973) and Weinman et (1973) 

derived the power function of the KP test. Weinman al. 


compared the KP tost with LR test, 
expression for the LR test as given 


They used the power function 
by Paulson (l94l). 


Recently, Tilca (l98l) generalised the KP test for 
left censored samples by proposing a tes. 
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U 


(X 


( 1 ) 


r^+1 




Ho obtained the nu.1.1 distribution of U • I'u 3 rther# he studied a 
Student s t upproximation of the null distribution to obtain 
approximated critical points, Khatri (1981) derived its non- 
nu] 1 distribution. There are some practical limitations for 
ovaluatinci the power function of the test by using the non— null 
pdf given by Khatri, As a particular case of our statistic/ the 
power function of Tiku^’s statistic (u) is given in Section 4- ,2, 


For the one- sided case wo have obtained the statistics 
*'lsE 5Dy equationsC 3 *1 *3) and (3,1*4) respectively - 

Botyi o.re ecfuivalcnt to statistic T given in equation (3*1,5), 

I'or testing : 6^ =02 against H 2 : / ©2 we may use eithcir 

(4*1.1) \r.E " ' ^^rJ^+l~^r2+l^/^'*^"’%’ " '^-"^il = ^1 

or 

(4.1.2) = lCx(,J^^-X^^|p/cT*-q 2 l=,T-q 2 l = (say), 

whore = bjL-k'2' ^2 “ 

ri+1 

b^ = E (n.-j+l)“^ and m^ = log { n ./(n .-r^^)} (i = 1/2)- 
"^1=1^ 

Hoto th-it/ q 2 = 0 if r^/n^ = r 2 /n 2 - Ho hav^ not boon able to 
identify cas<js whor(3 q-j = 0 except for the trivial case of 
ri sr ^ rin^ Tliuso st'^.tistics ctr^ of cnc foirm 

(4.1*3) v = 

1 ^ 
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where q is ci suits’s.blc constant which is either equal to or 
equal to q_ 2 * Without loss of generality, we can take q ^ 0, 
s i neci 




whore q 




r2+l r^+1 

and the distribution is obtained by interchanging 


the samples* 


iV.ir q =s 0, V reduces to the Tiku^s statistic U* Purthor/ 
f'.ir r^ z= r 2 = 0 and n^^ = n 2 the test statistics ’^lsE' ^MLE' 

.and the LR test given by Epstein and Tsao (1953) are equivalent. 

Thu test procedure is to reject against if 

V > Cq^/ where is determined so that P [v > Cq^I^q] = cc./ 
is the; chosen level of signif icancer. 

4.2* Distribution theory . 

This section is devoted to deriving the null and non~null 
distributions of the statistic V and its special cases. 


T he o rem 4 . 2*1 . For q > 0, the cdf of V under is 

^G^(c) , 0 < C < q 


(4*2.1) P [V < C|H^] 


lG 2 (c) , q < c < <»/ 


where 


G,(c) = HS (~1)^( .^)B(r 2 +a/£+j^ [{ l+hj( j Hq-c)} 


-d 


j=0 


{ i+h^ ( j )(q+cH 
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P_(c) x: 1-H .^)B(r2+l,f+j){l+h5_( j)(q+c)}~‘Vdh^( j) 


- 11 


JC 

^ j ^ ^ (q~c)} '’'Vdh 2 ( j ) 

aijcl the remaining notations are same as in Theorem 3»2»1» 


E.^iSSl!i* Note that V = lT~qi^ where T has pdf given in equation 
(3-2 *6)* I''or th.e random variable Y = T~q, the cdf is given by 

Fj^(y) for y < -q 
F 2 (y) for y > "q. 


p [Y < yiig 


where yj[_(y) = 

and p 2 (y) = 
Since V - iY|, 


H E (-”1 )^ ( .^)b( r^ +1, f+j ){ l-h-jC j )(q+y)}~’'V<^ho( j ) 
j=0 J ^ ^ 

. r 

l-H 2 (~l)^( /)B(r,+l,f+j)fl+h- ( j)Cq+y)}"'V<S5^i (j) 
j =0 J ■*• 

P [V < ciHq] = P [Y < cIHq] - P [Y < -ctH^ ] . 


# 


Consequently, for q > 0, 


*' [V < C|H„] 


0 for c < 0 

P 2 (c)~P 2 (~c) for 0 < c < q 
F 2 ( c)-Fj_(-c) for q < c < =°. 


On simplification, equation (4 •2*1) follows* 


The above theorem immediately gives Corollary 4»2*1 
for q = 0, which agrees with the null distribution derived by 
Tiku (1981) - 

C orollary 4*2.1 * For q = 0, the cdf of Tiku^s tost statistic U 
under is 



8.1 



for 0 < c < <=o<. 

Pnr the KP tost, wo have q = 0 and =0. The 

cdf is now given in Corollary 4*2.2. This agrees with the cdf 
obtained by Kumar and Patel (19 71) . 

Corollary 4.2.2. For q = 0, r^ = r 2 = 0, we have the cdf of KP 
test statistic as 

(4.2.?) P [u^ < ciH^] = l-{n2(l+n^c/d)~’^+n^(l+n2c/d)‘''^/(n^-:-n . 

for 0 < c < «>. 

The distribution of the statistic V under H 2 : 
can Ise obtained directly from the distribution of T given in 
Theorem '='.'^.1. The derivation of the non-null distribution of 

V is similar to that of the null distribution of V obtained in 
Theorem 4 *2 .1 . This distribution is given in the following 
theorem : 

Thc ox~om 4*2*2 * For (p = {6^-^6 2 ^/^ ^ ^ *3 ~ 

V under H 2 is given by 

(4.2.4) P [v < c|<p] = Pi(cl (p) 

''P 2 (q+cl <p)-F2^‘5”'Cl<p), O < c < q 

P 2 (q+cl <p)“Fj_(q"Ci <p), q < c < <=°/ 
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where F^(.|<p) and are given in equation O .3 •! ) • 

The equation (4t2.4) reduces to equation (4.2#1)^ the 
null cdf of V for (p = 0« For q = 0^ Theorem 4i»2»2 gives the 
non-null cdf of the Tiku's statistic. 

Corollary 4.2.'^ » The non— null cdf of the Tiku^ s statistic is 


(4 .2*5) 


P [^ < 



-exp{hj^( j )dcp} Q^{ 1 hj^( j )c}/dh^( j )+H 


^2 

E 

j=0 




) 


.B(r^+l,f+j) exp{~h 2 ( j )d(p} 1 h 2 ( j )c} 

~{ l+h 2 ( j )c}''*^]/dh 2 ( j) for 0 < c < <», 
where the notations are same as in Theorem 3*3.1. 


The following lemma can be used to obtain the non— null 
pdf of V as well as U* 


Lemma 


be as 
Qd( '1 

( i) 

( ii) 


4 .2 *1 . Lot Cj = d<p/c, cp > O and O^^Cxls) and L^(xis) 
dc-fincid in equation ( 3.1.2). The first derivative of 
. ) and L^(.i.) functions w-r- to c are given by 


Q^ic^\0) cf o Vc Kd)^ 

-ca,(H-h.c) ^ j 

i ^ ^ {c^(l+hj_c)} 


d e 

Q^Cc^lhj^c} = "•dhj^ 

1=0 


(l+hj^c)^"^^ il 


a -Cifl+hici 


■4* O 


/C JU) 


and 
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( i±i) ~ 


dh. 




d -c (l-h^c) {G. (l~h^c)} ^ 

[1 - E e ^ 2 (- A . -I ^ )] 

i«0 






for G l/h„* 


G TTd) 

££22.1: * •* i'!G proof follows at once by writing these functions as 

a sunt and then difforontiating term by term w»r» to c or by 
writing them as integrals and applying the Leibnitz rule for 
partial derivatives* 

From Corollary 4*2*3 and Lemma 4*2*1^ we got the non- 
null pdf of the Tihu's statistic U for (p > O as 


(4*2*6) frr(xx) = ri E .^)B(r,+l,f+j) 

j=0 ^ ^ 


-■H 


d * , ' ^2 . r 

£ Ui^Cl+h. ( j)u}^’~‘^''Vit + H S (-1)^'( ?)B(r.+l,f+j ) 

X X J _Q J 

Tr 


i^O 




2 » X' 

.xp {-h2(j)dcp}{l+h2|( + H _ E^(*-l)^' ( j^) 


j =0 


P f r< +1 , f +j ) -^2 C j ) -u., (3 


{l-h2(j)u} 


m c° 


r [ui(l-h,(j)u}] Vil] 
±=0 


where 


for u > 0, u l/h 2 (j) (j = 0 ,ly*..,r 2 )/ 

^ d ^(/u ^ at me singularity pomtS/ u « l/h^Cg ) 
the pdf fyCu) is taken as zero* 

It may be noted that, Khatri^s (l98l) expression for 
the pdf of U as given in Section 1 *5 » differs slightly from 
the above expression due to some minor integration errors* Also 
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the derivation of the cdf from the pdf given in equation (4*2 *6) 
i5> extremly difficult* Consequently, the evaluation of power 
function usincf Khatri/s approach is not easy, although one may 
resort to numerical integration techniques* But equation (4*2*5) 
yields the power function immediately* 

Corollary 4*2*'^ reduces to the following corollary for 

r-j^ ™ ” 0* 

Corollary 4*2*4* The non— null cdf of the KP test is given by 


(4*2*7) P [ < ci<p] = Q^( Cj_ 1 0)-n2 e3qc)(nj^<J>)Q^{ c^^l n^c/d}/(n^+n2) 


+ n^ e3<p(-n2<p) [L^{cj_^l n2c/d} -f l+njc/d}*" ]/(n^-!-n2) 

for 0 £ 

The cdf agrees with the cdf of the KP test derived by 
Weinman ^ (1973) and Dubey (1973) and given in Section 1*5» 

4 *3 • The LR test statistic * 

Lot Q ^ parameter space and 

^ Q ^ Q ^ 0 > 0] be the null hypothesis subset of Q » 

sot of tTO 

indepesndent type II doubly censored samples from E(ej_,a) (i = 1,2) 
with likelihood f sanction 


-dj 

A r. r 1 ft 
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exp [- ~ { 


!=r.+l J ^ ^ 


forei<4i)l < ... < (i =1-2)- 

VJe c 3 .ir>cuss the case > 0 , r2 ^ 0 in detail* Considerably 

simpleir expressions hold for other cases* From equations (2 • 4 *’^) 
and ( 2 » 4 '» 4 ), the ML estimators of 0^^8 2 and o in the parameter 
space Q are 

= ^r^li "■ ^ {nj^/(nj_~r^)} (i = 1,2) 

and 0 ^ X'iispoctively # Hence^ the maximum value of the 

likelihood function under q is 


(4 •3 .1) L(q) = 


i*) ^ ”i' yiyt 


In the' null hypothesis subset , the ML estimators of 9 and o 
for [soG, Section 2 . 6 ] are given by 


(4.3.2) == log {1 + ~rir~rrT:^ =^r?ii~^o + 


S - -r, 41-^0 tx . -o 


P-d C 7 


whore is the solution of the equation 


and trie notations 


p„a* g^ ^ cl* g^ + fY-P 


are same as in Soction 2 .6 


Prom equation ( 4*3 » 2 ) wo immediately get 


^ = nvea* a„+Cf+ri)Y-Pl 

(J 

o 


l-*oxp { 
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1-exp {- ^ = r2Y/{P-d* + r2Y} 


Further, 


^ il • **^3 • 

S [ r, (X^.^^*-0 ) 4- S.(X^^^ ~0 )]= -—{P+f \ -0 )} 


X 1 


Then the rndxlinum value of the likelihood function under 63 is 


LiH) =[ n 


i=l '"r d*a^+(rj_-ff)Y-] 


r- ^ r_Y 

■1 "■ [— rl- 


P-d^c7„+roY 


cfS +fY~P 


'd*0^+(rj_4-f)Y~P 


]' e^(-P/ 5 „) for Y = > o. 


12 ) ^( 1 ) 


I'low/ the likelihood ratio is given by 

A ^ A ^ 

X = L( 63 )/L(q) 


( 4 *3 *3 ) 


whore Const 


r "i“T^ 

- , .a d* Y ^ ^(d*a +fY-P)^ exp (-P/a ) 

Const - { ~} o ^ ' o 

{d^"aQ+(r 3 L+f)Y-P}^l'^^ {P-d^a^+r 2 Y }"^2 

for Y > 0 , 

2 n, -(n.-r.) 

. = cxp(d*)- 71 n^^Cn^-r^) "■ "" - 


For Y = ^ ^ obtained by replacing n^,n2/r3_,r2 

by n2#n^# ^2*^1 t’ospectively in ecfiBtion ( 4 » 3 » 3 )- The LR test 
then rejC'Cts if X < ^(x' where x^^ is chosen so that 
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4.4- Critical points of the tests U and X . 

4-4--1- ii^act cx'itical points of the test statistics- The 
ciii-ical point c^^ of V is obtained by solving the equation 


(4-4-1) P [V > c^lHj = a, 

v;hero a is the chosen level o£ significance- From equanion 
( 4-2-1 ) i^ 'is clear that, is either the solution of 

^1 1 r ^ 

(4.4-2) H I^(-'l)-^( j )B(r2+l/ f+j) [{l+h^( j)(q*-CQ^)) 

~ {l+hj_( j)(q+CQ^)}’''^]/dh^( j) = 1-a 
(case in which < q) or the solution of 

(4-4.3) H [ Z (-1)^( MB(r2+l.f+j) { l+h^( j ) (q+%)r /dh^( j ) 


+ I (-l)i''hB(r,+l,f+j){l-h2(j)Cq-c„))“'V®2(j)] =a 
j=0 J 

(case in which c^ > q) according as > l-cc, or < l~a 
respectively/ whore 

^1 ’ ^ 

= P [V < qlH^] = ^ i;^(-l)^'( j^)B(r2+l/f+i) 


- [l~{ l+2h^(j)ql“^]/<^l^l 


The value of can thus be obtained by first evaluating P^ and 
then applying Newton- Raphson method to the relevant qu 
unlike the one-sided test statistic discussed in Chapter III/ 
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huro it i .s not possible to get a compact expression for Cq. even 
for = r 2 = 0 case* 

i!ow/ critical points c^^\c^^^ and c^^^ of the tests 

cx cx cx 

V^/V,^ >ind TJ arv:; nothing but the solution of the equation (4*4*l) 
vrith q = ^ respectively* It is clear from the 

orfial'.ior. (4*2*1)/ that thu critical point c^^ of the test U a*? 
Gim:|:>lY thv: solution of the eq[>.iation (4*4*3) with q = O* Some 
of th'..Ge critical points arc tabulated in Table 4*4*1 for 


studying the; nerformance of tlrieso satisfies. 

.'Jjru’'. hlie Lfi tust statistic X given in equation (4*3*3) 
is v.iry coinplientod, wo have not attempted to derive the 
distribution of' X * Howovor, we have obtained simulated critical 
points Xq. based on 10/000 iterations/ for studying the 
p>,rformanc... of thi. statistic. These are tabulated in Table 4.4*1 . 

i‘'or r, - r. = 0, the tabulated critical points are exact* 

1 2 

Jn .'ioctlon 4*5, we have conpared the performance of 
thC'iO etatietios. It Is very difficult to compare the power 
functions of the test statistics. But Tables 4.5.1, 4.5.2 and 
sonn other numerical calculations show that, in between Vj, and 

there is very little difference of the power values. However, 
is nliahtly blotter than V^- H®oe we are^suggesting a method 

for obtaining an approximate value eg for . ihis 

value for solving equation (4-4*l) 
also bo t.iXon as an initial value lox u y ^ 

for c‘lb An idontioal approach can be used for c^ and 

cc 



In Section 3 *4 / 


4.4 -2 • Approximated critical point c^'’ of V, • 
we studied an approximate null distribution of T. In the 
present case^ let so that ' Then from 

equation (3.4»7) wo have 

(4.4-4) F,(Y^) = q^/(d~l),Var(Y^)=d^' C (d~l)A+qJ/[(d-l)^(d-2)] 

with = '3. Wow, it is clear fi'om Section 3-4 that Y^/'/A has a 
Student^ s t distribution v/ith d DF» We improve it slightly by 
considering tw constants a(> o) and b such that mean and varianc-e 
of the random variable a(Y^+b) is equal to the exact mean and 
variunci:- nf a Student^ s t random variable with d DF. Tha^ is^ 


E [a(Y^+b) ] = 0 and Var [a(Y^+b) ] = (V(d~2) . 
Then on u.<e no equation (4-4*4) we got 


( 4 .4 -5 ) 


S incx- 


a = (d-l)/[d{ b = -q^/(d~l) - 

p [Vj > = 1 - P[iqi < 

= l-P [a(-c(^hb)<a(Y-+b)<a(o^^hb) ] 

Ca/ «** 


( 4 •‘i ) 


= 2P [Ta > ' 

b = 


l) 4- /a 

% = 


■a^ 


Whor.:. Is th.! (a/2)th parccntilo point of Studont's t 
dlstrilxition with d DP. Soms of tho exact critical points 
[solution of thi oration (4.4.1)] and its approxi^aatod value 
cj [given in the negation (4.4.6)] arc- tabulated in Tabic 4.4.2 
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foe different combinations of sample sizes and censoring patterns. 
Mote that, for very small values of d and | ( r/n^-r/nj) , > o.i, 
the approximation is not satisf aotor-.r, 

oxm.ij.:ir to the approach discussed in Section 3.4, 

w-: also studied a normal approximation for the critical points 
of V^. But those are not as good is the Student's t approximated 
vo.hios for most of the casos- H..nce, these values are not 
t-ibulatod. 


4 -5 . Power of the tests ^ ^ 


The power of the tost V for testing s = 82 against 
H 2 : ^2 given by 

l-Pl(c^l<P) for 9 > 0 


( 4 .5 .1) p [y > c |(p] = 


|^ 1 *'P 2 ( for (p < O, 


where c^ is tho exact critical point, given by 

equation (4 •2.4) and obtained by interchanging 

n,, , n^/ r,j , r 2 by ecni^tion (4*2*4) and evaluating it 

for t (PI . 


Tli(; power functions of ^ obtained from 

equation (4-5*l) by replacing (q,c^) by 

and respectively* Some of these power values are 

tojoulatod in Tables 4.5*1 and 4*5*2* 


Th*e power 'if the LR test statistic X is given by 
P [X < Cjj^Kp]. As wo had mf-mtioned in Section 4-4, the distribution 
theory of X is very complicated* Henco we obtained the povrer 
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by applying Monte-Carlo technique, using 1000 iterations. 

However^ for r^ = = 0, the exact power was obtained by using 

tne power function expression as given by Paulson (l 94 l) [see, 
also Weinman et 1973]* 

T<or the following sets of sample sizes, censoring patterns, 
and various values of g wc studied the comparative performance 
of all thu four test statistics for a = 0.05. 

(i) n^=lO, n 2 = 6 , d=l 4 ,(r^/r 2 )=( 0 , 0 ) and (r^/r 2 )=(l,l) in 
Ta}>le 4.5*1. 

(ii) ’.' 1 ^= 112 = 15 , rj^=l,r 2 =» 3 , d = 18 and d = 24 in Table 4.5*2. 

To hioh light the shape of the power functions, three such cun^'-js 
for a n 2 - 15, r^ = 1 , r 2 = 3 and d = 24 arc drawn in 
Figuri.„ 4 *5 *1 . 

Since the difference in power values of and "^2 is very 
small, w*.; have not drawn the power function of V 2 in Figure 4*5 ■•■i. 
In all casos studied, it was observed that the "average'' power 
of was always greater than the "average" power of ^ 2 * where 
the average was taken of the two values corresponding to and 

-cp . 

A normal approximation similar to the one studied in 
3ocrt:ion 3*5, for power function is also considered. Since 
= T-q^, the power of the test is given by 

P [V^ > > 4^'] 

= f [' 5^1 & + p -4^*3 



( 4 - 13 - 2 ) 




y 2 


P [VL > o] + P [ w, < o] 


where W. 


4jir42lr^‘5i+4^^)a^ and 


.( 1 ) 






.(1) 


a 


)a 


* 


Similcir to the results obtained in Section 3-5, v/e see that 
W^. (i = 1./2) has an asyrrptotic normal distribution with 


= (tp - = p^a, 

Var(W^) =-- 

Edi^) = (<P + c^^ber 5 

and VarCw^) = { A+(q^«c^^^ ) s g^a^. 

Hence from er^fuation (4-5*2)/ we immediately get 


(4*5-3) P ^ + ^C~P2/'^g2^* 

The exact power values obtained from equation (4-5-1) and 
its approximated values given in equation (4-5-3) are tabulated 
in T.able 4-5-3 for a = 0-05 and selected values of n^/n 2 /r^/r 2 

and d- 

Using the po«r function as a tese, the following concluaions 

may b>,- drawn- 

Ca) All the throe teats V,,V2 and U aro biased. Ho^.ever, the 

extent of the bias is different for different tests. 

(b) There is a considerable loss of power for all tests even 
with mild censoring on the left as is evident from 
Table t.5.1. Bit the loss of power is negligible with 

variations in censoring on the right- 



The test statistic U proposed by Tilcu (l9Sl) is more 
biased than and In general the statistic U shows 

poor performance as iCr^^/n^-r 
tne bias of U seems to increase with increase in d. 

Table 4 «S <*1 shows that the test statistic V 2 is rcilativcly 
nx 3 ro, biased than . 

V/o do not have sutficient evidence to conclude that, the 
I;R tost statistic X is unbiased for > 0 and/or r 2 ^ Q/ 
since some simulated values are less than a as in 
'j'ablo 4 #5 •! - However, for r^^ = r 2 = 0/ IXibey (1973) 
established that, the LR test statistic is unbiased [see. 


2 ^n 2 ) I increases. Further, 


also Kh.itri 1974]. 


T.*ble ^.b.l. Tabic 4 .S *2 and Figure 4.5*1 show that there, 
is very little difference in the oower values of X and 
- Since tho statistic X given in egaation (4.3*4) is 
voiy complicated while the statistic given in equatio*^ 
(4.1*'=') is considerably simple, we strongly recommend the 
use of t<2Ct statistic in such situations. 

The normal approximation for the power function as fairly 
good if =D (say) is small, even for small 

vJuos of nj, and n^. Howwar for large values of D, th.. 
approximation is not that good even with moderately large 


values of n, and n,. This is well illustrated in Table 
4 . 5 . 3 , for the caso nj^ = 15, n^ = 25 and d = 28 . lor 


Cr ,r ) = (0,4), the maximum difference (among all tahulet 
vo,Le!) between approximate and exact values is 0.0274 for 
» = 0 . 30 , whereas for (ri,r2)=C4,o), tho maximum dlffercuc 


is 0.0566 for <P = 


.40* 
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table 4.5.]^, Exact 


tests and exact or 

siraulated poweiTBt \ for a = oIOB, n,--=10, 

n2=*5 and d = 14. ““ — ~= 


CP 


=r^=0 


v^~" 

V2=u 

X 

-1 *00 

*9977 

-998 2 

•9971 

-0*80 

-OB 33 

•9871 

.9795 

-0 -GO 

.B9 27 

-913 2 

.8739 

-0 *50 

*7619 

-798 2 

•7308 

-0 •4Ci 

.'55 24 

-5987 

•5158 

-0 *30 

-3X94 

-357 2 

• 2915 

“0 • 2b 

-2249 

-25 3 6 

• 204 2 

~U • 20 

-1545 

.1743 

-1403 

-0 *18 

•1064 

•119 2 

•097 2 

— C *1/! 

.0755 

-083 2 

.0700 

-0 *OG 

-0575 

•061 2 

•0549 

0 .Qf'i 

.0500 

•0500 

.0500 

C' *05 

-05 23 

•048 3 

.0553 

0 -10 

.0661 

-0567 

-07 29 

0 -IS 

.0055 

-078 2 

•1079 

0 -20 

-147 2 

-1130 

•1679 

0 • 2S 

: .2276 

-13 27 

• 258 6 

0 .30 

-3358 

-2746 

•3761 

/» f** 

. * j 

-•rf j' ') 2 

-3879 

.505 2 

0 .40 

-5341 

•5*990 

.6 273 

c 

0 

.73 20 

-7 243 

.3121 

Cj .60 

-=39(54 

•8631 

-9125 

0 *80 

-9707 

-9710 

•93 21 

1 -0C1 

*9957 

-9941 

-9964 


r^ ^r 


^1 

^2 

U 


.9 388 

♦9490 

•9561 

.9 23 

-8089 

•8323 

•8498 

.779 

.5580 

•5907 

•6169 

.535 

•4099 

•4404 

*4657 

•403 

•2759 

.3001 

.3 205 

• 254 

-17 29 

•1894 

• 2036 

•155 

•1346 

.1476 

•1588 

•124 

•1045 

•1143 

•1229 

.083 

•08 18 

•0888 

.0950 

-078 

•0655 

•0701 

.0741 

.057 

-0550 

.0573 

.0594 

.05 2 

.0500 

.0500 

.0500 

.048 

-0504 

.0479 

.0457 

.041 

.0566 

•0512 

•0464 

.049 

.0696 

.0605 

.05 25 

•014 

.0905 

.0769 

.0648 

•076 

•1209 

•1017 

•0843 

•106 

-16 20 

*136 2 

•1125 

-173 

-2147 

•1315 

•1505 

-221 

-278 2 

-2376 

•1938 

.^8 

•4Z35 

.3765 

.3 239 

.455 

•5356 

•5 312 

•47 23 

.613 

-3 258 

•789 2 

-7443 

•835 

.9403 

-9 24 2 

.90 23 

.941 


^Simiilatrrid power based on 1000 samples for each sample size* 
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TABLE 4 *5 *2* 

Exact 

power ( 

of the 

tests Vj 


and simulated 

power 

of the 

test X for CL 

= 6.05 


=15, 


T 

and 

3^ 





(p 

d=18 


d= 

24 




^2 

U 

A 

^1 

^2 

U 

“TID — 

A 

-1 *00 

*9960 

*9964 

•99 70 

*997 

*9973 

.9970 

-9967 

-994 

-0 *90 

•^9 3 S 7 

*9871 

-99 21 

-987 

•9919 

•99 24 

.9949 

-939 

0 #3 0 

•963 2 

•9650 

•978 2 

-950 

*975 2 

•9765 

•9861 

.973 

-0 *70 

♦9 1 4-jB 

*9163 

•9435 

•390 

.9 336 

-9 367 

.9596 

««96 

-0 *60 

•8177 

-8 234 

•870 2 

•768 

•8461 

•8517 

•3966 

*779 


•6700 

•6773 

-74 21 

*57 3 

*6978 

.7058 

.7741 

*588 

— 0 *40 

♦<1-S 6 3 

-4933 

•5641 

.386 

•5045 

*5130 

•5910 

.396 

~0 .3f) 

•3073 

*3134 

•37 24 

• 250 

-3143 

•3 212 

•3877 

.219 

-0*20 : 

•1703 

•1741 

*2130 

•112 

•1710 

*1754 

• 2189 

.113 

“0 *15 

•1 221 

•1249 

•15 36 

*103 

•1216 

*1248 

•1567 

.105 

-0 *10 

*0364 

*088 3 

• 1077 

.07 2 

•0864 

.08 84 

•109 2 

.064 

n *00 

•0500 

•0500 

•0500 

.05 3 

.0500 

.0500 

.0500 

*056 

0 •05 

*0474 

.0463 

.0338 

•058 

.0485 

.0472 

.0336 

.066 i 

0 *10 

•056J 

•0541 

.0237 

*069 

•0596 

.0566 

.0233 

.093 

0 *15 

•03 14 

-0764 

•0135 

•131 

.0869 

.0315 

.0180 

-173 

0*20 

*1 26 2 

•1180 

-0179 

,*191 

•1358 

•1269 

-017 2 

.223 

0 -30 

• 2032 

-27 25 

•0331 

.374 

•3097 

*2931 

.0319 

.418 

0 ..l-O 

*5143 

•4955 

-0843 

*618 

.5 430 

•5 244 

.0337 

.669 

0 .BO 

•7 206 

.7054 

*19 27 

.795 ' 

•7451 

.7310 

*1973 

.3 25 

0 *60 

-3597 

*0503 

• 359 3 

*904 ' 

•875 3 

.8670 

*3763 

.9 26 

0 -70 

•9 36 3 

.9314 

*5507 

.953 

*9445 

•9404 

•5793 

.967 

0 *1C 

•9732 

*9709 

*7 210 

.907 

*9771 

.975 2 

•75 24 

.933 

0 *90 

•9 :9 3 

*9384 

-046 2 

•995 ■ 

: *9910 

.990 2 

*3707 

.99 3 

1 *00 

*9960 

.9959 

.9 229 

*995 

*9966 

*9963 

.9338 

.999 


*.SijtTu] ated power based on 1000 samples for each sanple size* 



TAB]jE 4*S*3* Exact and approximate power valu es of the test 
St at istic for a = 0,05 « 


(nj_=n2=10,del6) 


(p 


r^=0, r2=l 


Exact Approx- 

Exact Approx. 

~1 .00 
-0 .90 
~0 .80 
~0 -70 
-0 .60 

~0 .50 
.40 
-0 .30 
-0-25 
-C .2t) 

-0 -15 
-0 -10 
-0 .05 ! 

0 .00 

0 .05 

C -10 
' 0-15 

0 -20 

0 -25 

0 -30 

0 *40 

0 -SO 

0 .60 

0 -70 

0 .80 

0 .9 0 
t -00 

-9786 .9832 

.9551 .9568 

-9093 -9038 

.8 293 .8 140 

•7045 .6351 

.5 381 .5 284 

.3593 .367 2 

-2110 .2271 

.1560 -1709 

•1143 -1253 

.0344 -0905 

.0647 .0663 

.0536 .05 21 

-0500 .0474 

.05 36 .05 21 1 

.0647 .0663 i 

.0844 .0905 

.3143 -125 3 

.1,560 .1709 

.2110 .2271 

.3598 .367 2 

.5 381 .5 284 

.7045 .6851 

.8 293 .3 140 

.9098 .9033 

-9551 .9568 

.9736 .933 2 

.99 38 .9956 

-9841 .9851 

-9606 .9577 

.909 2 .6990 

.3 1 09 .7960 

.6555 .6475 

.46 24 .47 21 

.2313 .3021 

.2094 - 2291 

.15 26 .1675 

.110 2 .1186 

.0301 .03 25 

.0605 .0590 

.0500 .0473 

.0481 .0490 

.0558 .06 34 

.076 2 .09 25 

.1145 .1377 

.1774 .1996 

.2675 *2776 

*4967 .4683 

-7066 .667 2 

.8456 .8 275 

.9 250 .9 273 

.9649 .9755 

.9840 .99 84 

.99 29 .9936 
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TABLE 4 *5 *3 • Contd . 


^ 0.60 
>0-55 
-0-50 
-0 -45 
-0 -40 


( n^=15,n2=25, cl=:23 ) 


r^=4,r2=0 

r^=s0, rj=4 

Exact Approx - 

Exact Approx - 


•3719 
-3139 
-75 22 
-6673 
-5657 


-3717 

•3041 

•7185 

•6131 

•5091 


-996 3 
.9739 
•9546 
-9141 
-8443 


.9913 
•9735 
•95 27 
•9065 
•8332 


-0 .35 
-0 -30 
-0 -25 
-0 -20 
-0 -15 


•4511 
•3314 
• 2185 
•1264 
•066 2 


• 39 96 

• 2979 

• 2108 
•14 24 
•0941 


-7412 
-6075 
•4538 
-3179 
-203 2 


•7306 
-6034 
•4641 
-3291 
• 2136 


■0 -10 
■0 .05 
0 .00 
0t05 
0 *10 


-0399 

-0334 

-0500 

.0701 

•0984 


•0648 
.05 26 
•055 2 
-0714 
•1003 


•1222 

.0729 

-0500 

.0507 

-0796 


-1 269 
-0715 
•0464 
.05 21 
•0931 


0 -15 
C -20 
0 -25 
0 *30 
0 *35 

0 -40 
0*45 
C -50 
0 -55 
0 -60 


•136 2 
-1B47 
-2448 
•3166 
-3985 

•4875 
•579 2 
.6634 
•7499 
•8199 


•1418 
•1957 
• 2614 
•3371 
•420 2 

•5070 
♦59 35 
.6757 
•7500 
•3142 


•15 33 
• 2913 
.4722 
.6490 
.7379 

-8814 

.9377 

-9639 

-9851 

-9931 


.1758 
• 3012 
.4577 
•6 216 
-7659 

•3732 

-9403 

•9753 

•9916 

.9975 
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FIGURE 4-5*1. Power functions of V-i^U and X for n^ rn2r15 =1,r2=3 and d=24 



chapter V 


GEI'IEIlAJj 11-3 ED STATISTICS FOR K RIGHT CENSORED SAMPLES 


5 • 1 * Introduction • 

oo f-ar wG considered the case of testing the equality 

of location parameters of two exponential distributions. In 

this chaptei:, we px'opose two test statistics for testing 

^^o * “ ‘■'^2 = = e against the alternative hypotheses 

IL : 0. > max (0.) and : at least one 0^ (i = 1/2 /..-,k) 

2<j<K -J i 

is different from 0, based on K(> 3) independent right censored 
sain}:)lc.'S . The necessary distribution theory of the proposed test 
statistics is discussed. Some critical points and power values 
art.; tabulated. Finally# we compare the performance of these 
statistics with the test statistics proposed by Khatri (1974) 
and .Cingh (1983). 

Lot ^ 1/2# ---#K) be K independoit 

sarapl-s fran E(G^#a)# v;hcro > l. For simplicity of notations 

lot x (i = I/2/.../K) and ^ ) = min(X^#X 2 / .* ./^j^;)/ that 

is# X^ is tiu; minimum of the i th sample and is the minimum 

of all th.„ obsor’/ations. 

Khatri (1974) derived the LR test for testing against 
H 2 . It is given by 

K 




i=l 


(5 -5 . 1 ) 


(5-1-2) d0* 


arci same as in 
of tlr:i LR test 


K 


(i) - .. 41 ) K 


^ 3.^X 

earlier chapters* He obtained the power function 
and showed that 


(5 *1-3) P [U^_ > c|H^] = {B(d,K~l)}“l dy» 

e 

Fror.-i ocp lation (5*1-3), it is easy to shov/ that dUj^/(K--i) has 
an F~distrlbution with 2(K~l) and 2d DF, which is denoted by 
^'2(K~i),/d* <lcrivod two union intersection test 

Gtatj sties from two different view points- These are given by 

(5*1*4) = [max{n^(x£-x^j^ P /n2(x2-XQ) ), » ..,nj^(3{^~x^^) )}]/da"' 

and 

(5*1*5) = [m‘ix{n2(x2”X^)# -*-/njr(xj^x^),n^(xj^-X2)/ **• 

n^(xj^~xj^)}]/da^* 

In all three cases, the test procedure is to reject if 
Ui > otherwise accept where the constant c^^ (i=l,2, 3) 

is dotorminc'd by solving the equation P [u^ > = a* 

Khatri (1974) also obtained the null distributions of U 2 and U^, 
and th..-ir power functions under the assunptions 6^ > 
v/hich can bc' <;chievcd by renaming the populations. He has 
provided s;orne critical points of ^2 and U^* However, no powc:r 
comparison studios for U^, U 2 and been done* 
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Singh (1983) also discussed the LR test procedure for 
testing against The test procedure is equivalent to 

relectlng if where 

K 

^4 (K-i)a*} and P [u^ > c^'^^ ] = a. 

fie has nhown that, the distribution of is 2d ' 

did not study the power function of U^# Note that^U^ = dU^/(K-i) 
and tno power function of has been already given by 
Khatri (1974)- 

Although there are number of tests for testing 
agaiur/t Ji 2 based on some theoretical considerations as describ'd 
above., we propose another test based on 

T = { max X. ~ min XA/a*o 

This is mainly for comparing the performance of various tests - 
Iharti’tT, a generalization of T 2 for left censoring is easier 
than that of the other tests- Note that, if all samples are of 
equal sizes, then T 2 is equivalent to U 2 » 

5 * 2 - To st statistics and their null distributions - 

Similar to the one-sided test statistics of the two- 
c ample case, v/e propose 

(5-2-1) = {X^ - min(X2,X3, ---,Xj^)}/0^ 

for testing against 
H 2 # we propose 


Similarly, for testing against 
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(5 '"2 *2) T„ = { max x.- ~ min x.}/a*« 

l<i<K ^ i<i^ 

The test procedure is to reject if T^ > c^., where 

^ ^ ^ ~ i/2)/ and a is the chosen level of 

s igni f j.curace <• For obtaining the critical ooints c^ and C 2 / v^re 
nc>cd tno distri 'out ions of and T 2 under These are 

obtainu<l ffrom the follov/ing lemmas s 

Lemma 5»2»1 * arc K independent random variables 

with pdf oI given by 


(5 *2 *3) (w) = n^ exp (~nj^w)/ w > 0 (i = 1/2, ••(•/K), 

ths;n the* pdf of Z = VJ^~rain(W 2 / • « VJj^) is 

n^(H--n^) exp {(N~n^)z}/N/ z <0 


(5-2 -4) f(z) 


n-(N-n. ) exp {-n. z}/N/ 


z > 0/ 


I A4,^ \ XM C 

wheri'i II n^+n2+ * • *'+njj.* 

P roof - ii'.t Y rr rnin (W 2 /^ 3 / “’■-/Wj,)- Then the caf of Y is given 
by 


K 


K 


= ^ Gxp(-n^y) = l~exp(“N 2 y^ /Where N 2 - S n^. 

~ i=2 ' ^-2 


Consct.jucntly, the pdf of Y is 

f(y) = N 2 exp ('-N 2 y), y > 0- 
The result now follows/ on applying Lemma 3-2-1- 

Lcntca S.2.2. Let ( 1 = 1. 2 be K indepondbnt randonp, 

v, 4 rlati..s with pdf given by aquation (5.2-3). Then tho pdf of 
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Z = max Wj ~ min W. is 
l<i<K l<i<K ^ 


(5 *2-5 ) 


g( z) 



K -n . z 

j=l,j^i ^ 


r\ { l-exp(-n, z)} , z >0 

h=l^}vii, j 


wru-'r«j IJ ■■ n^4-n2+ * . •+nj^, 

P roo L ■ '-■'in CO V/^ (i ts 1^2/ »"’/K) arc K independent random 

V- ri'-itof! v;ith pdt (w), honco^tho cdf of Z is [soe^ David 1981/ 

1 

p • 26 J givi-n by 


K o® K 

(5»2«6) G(z) = r. / -f^ (w) T\ (w+z)~F^ (w)}dw/ 

i=l -CO i j=l/j/^i j j 

where- for xj > 0, (w) = .w) is the cdf of Wj • 


Subs'ti tutiriff f-.-r (w) and F (x^), v7o havo 

i j 

K oo -n.w K -n .w -n . (w+z) 

G(z) = E ; n.c n {c ^ -G ^ }dw 


1=1 o 


3=1/ j/i 


K oo 


-Nw 


S J n.e‘'-' n 
i=l o j =l/j/^i 


K_ -n.z 

{ l-'G } dw 


^ ’^i ^ 

= E ^ n 


-n.z 
{1-e -J } 


i=l ' 


(5 .2.7) 


E G.(z)/ 
i=:l 


where G.|(z) = -s" FT {1 


n * K 

i TT — exp ( —n J z )} (1 = 1/2/»'*/K) 


K 


NOW/ log G/(z) = Const. + ^ 


log { l-exp(-n^z)} 
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Diffcrv ntiating w-r- to z, on both the sides of the obovo 
cquativ'^n, wc havtj 

= 5 r / ^ s 

Cj j ( z) L • z//{ 1- exp (-n. z)}*] , where g.(z)= — rg — 

-> J j -J “"x 9z 


Hence, 


K K 

■tf 2 n. expC-n.z) 7"! 

^ ^ h:=l,h;6i,j 


From cKiuv.tion (5.2*7), the pdf of Z is given by 


{ l-exp(-n|^z)} - 


K 

g(z) = S gHz) 
i=l 

x^rhich givers the required equation (5.2*6). 


T heorg g 5*2 *1 » Under the null hypothesis the statistic 
doflned in equation (5 *2*1), has tbe following cdf : 


(5-2-3) P[T^<clH^] 


n^{l-(N~n^)c/d}~V'^, c < 0 
l-(N-n^)(l+nj^c/d) 'Vn/ c > 0, 


K K 

whore IT S n,. and d = S (n-~s.~l)» 
i=l ^ i=l ^ 

Proof - The random variable = (Xj_~6j)/o (i = 1,2, ---/K) has 

the pdf civen in equation (5-2* ■^)» Consequently, Z = W^-' min 

2<i<K 

follows the distribution as given in equation (5-2-4)- Then 
under = d2/VJ, where W = dec^/q* Now, applying Theorem 3-2-1 

xfo got tbo required equation (5»2-8). 

For general K and the distribution theory of the 
statistic T 2 is very complicated, although it is possible to 
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follow rhe same approach. We therefore consider the simplifying 
assumption K = ^ in Theorem 5.2*2. 


T hi'jorcm 5 > 2*2. 


For K = 3 ^ the null pdf of is given by 


C5.2.9) f(tj) = ^ l: nj^(N-nj^)f(l+r, tVdr‘^-l_(l+ I i^tVdr* 

h=l/fv^i ^ 

for t 2 > 0 , 

where W 

£ ££ ££* I'^ote that/ = {X£~6^)/a (i = 1 / 2 , 3 ) has the pdf given 
in orfu.-ition (5.2.3). From Lemma 5.2*2, for K = 3, the pdf of 
Z = m'C-:( W^/ W 2 / Wj) ~ min ^ rex^ 7 ritten as 

1 ^ ^ 

= ■^ S nj(N~nj ){e3{p(~n .z)-e 3 {p(~ E n, z)}, z > 0 . 

■ i=l ■ h=l/}v^i ^ 

Since T 2 - dZ/W, where W = dcr^/a, wo obtain the desired pdf of 
T 2 by proccodinq on linos similar to that of the proof of 
Thcor^jtn 3 . 2 . 1 . 


If 

in cqU'ition 


Of = n 2 = ... = = n, then the statistic T 

( 5 . 2 . 2 ) is equivalent to T^, v;here 


2 


given 


(5 .2.10) T 3 = nT 2 = dU 2 . 

For ecjual soiAjOle sizes, V 7 e therefore use T^. The null 
distribution of T 3 can be obtained as before. It can also be 
derived from the cdf of U 2 as given by Khatri (1974). Note 
that/Khatri has used v/s for U 2 and p for d. Then by using 
Khatri' s expression we have 
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p[T 3 < CIH^] = P[U 2 < 

(5.2.11) = 1- i^-i)^(^"J'){i+Cj+i)c/a}“‘^, o>o 

3 =0 ^ 

For studying the porformance of tiiese statistics, the 
non-null distributions are derived in the next section* For 
K > 3 / the exprcssi>’ns for pot^/or functions are very complicated* 
Hence, We hav.j only considered the case K = 3 in detail- 
Without loss of generality, we have assumed 9 ^ ^ ^2 ~ 
this can be achieved by relabeling the sample* 


5 »3 • Hon- null distributions of the statistics * 

The non-null distribution of T^ is derived by using 
th'.- following lemma : 

L. .imma 5 * 3 * 1 * Lot W^,W2 and W3 be independent random variates 
with Pdf of Wj_ given by 


(5-3-1) (w) = C!xp {-nj_(w-(pj_)} / w > <Pi (1 - 1/2/3^/ 

i 

where -Pj. > <P2 > <P3- Then tba pdf of z = vl^-mlnCw^^Wj) Is 

(n2+n3)b^ oxpi (n2+U3) 2} / -<=° < z < <^2 

n^b2 exp(-n^z)+n3b3 exp(n 32 ), a2 < z < dg 

n3_(b2+b^) exp (-n^^s)/ a3 < 2 < 

where aj = (<Pi-<?j) (3 = h = "l V(hl+'> 2 +" 3 >' 

bj = a"2 {nj^(l2-n3®3+n3'‘2'/an2+n3)(ni+“2‘‘’"3'’ ' 

bj = n, o,.p C-n3a3)/(npn3) and b^ = n3 e^(n,a3)/(n,+n3). 


(5-3.2) fCz) 
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cdf of Y = min 

P[Y<y]=i^j3_^p < y]3 {1 ~ P [W3 < y]}, 

1 ^ rom Lain, Lha pdf of Y given by 

^3 oxin{-.n3(y~(p3)}, CP 3 < y < 

( n2+n3) wxp{-(n2+n3)y-(-(n2^2+n3^3^a (^2 < y < oo„ 


f (y) 


Ij'-jViT, 1. trf^’iTi th(! jpdf of Y and and mvaking a transformation 
y. rr. — y and vr = V 7 , wo get the marginal pdf of Z as 


00 


(n.3.3) f( 7 . ) 




Pj^( z,y)dy, 


*0= < z < a. 


<P. 


S p^(z,y)dy +/ Po^Z/yay, CC2 < z < C&3 
<^2 <p^~z * 

^2 

/ paz/y)dy + ; p^(z/y)dy^ a- < z < «>, 
(^2 ^3 ^ 


3 3 

\/hcrc,- p^(z,y) ^ nj^(n2+n3) exp (-n^z + ^i^i 


and ppCz^y) = Gxp + n 3 <P 3 -(nj^+n 3 )yl • 

Jiq’aatl'>n (5*3-3) now gives the required equation (5»3«2) on 
s Impl 1 1 ication * 

Th.-xjr^^a Por 0^ > ©2 ~ ®3' non-~null distribution of 

is given by 


p < C 10 ] = p(cie) 



Ill 


(5 «3*4) 


{l-*(n2+n3)c/d}*"^ , c < O 

Q^Ccsio) + L^{ C2I (n2+n3)c/d} 

~^3 [b^(c2l n^c/d) ~ L^(c 3 I n^c/d)] 

~Cb2Qd{c2ln^c/d}+b4Qd{c3inic/d}] , c > 0, 


where = da^/c, a. = {G^^e^/a (j = 2 ,^), Q^( M *). V-l.) 
■are qiven in equation (3*1 *2) and (i = 1^2/ 3^4) are given in 
LorrBnri 5* 3*1# 


£E2£1£* bet = ^i/*^ = 1/2/3)* Then are independent 

random variates with pdf given as in equation (5*3*1) with 
^i ~ of Z = W^-min( W 2 / v*J 3 ) is given in Lenuna 5*3*1 * 

Civ;arly, can bo v/ritton as dZ/V7, where W = d£7*/cT • Now/ us ana 
simiT az' arguments as in Theorem 3*3*1/ the cdf of T^ upto the 
point G is the integral of the joint density of(T^,w) over the 
shaded region shovm in Figure 5 * 3 » 1 * 



FIGURE 5*3*1* Region showing the cdf of T^^ upto the 

point c* 
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Consociuently, 

00 c 

^ ^ '*^1 — ~ ^ ') dw = J^ (say) for c < 0 

O ~oo A _ i - 


and 


00 O ~2 c 

['■*-"1 < ^•= -T (s p^ra-t-)dw+; u p.d-L. )dv 

o —00 -• -^ ^ ^ - 1 . J. 


2 c 

is 

o o 


00 '^CCy'w 


'3 c 


1- / (/ p^dt^_)dw -f / is P2dt^)dw 


^ 2 ° 


'2 ^^2/ 


w 


day'v/ 

+ -f is p-jdu. )dw + S is PoClt. )dw 
day'w ^ /._ - x 


ctog/w 


'”^ 2 '*^ 3 ’*^ 4 "*^ 5 *^ 6 "^7 (say) for c > 0 , 


v/h>.ro Pj = (n 2 +n 2 )b^ oxp [ -w{ I“>(n 2 +n 2 )t^/d} ] w^/d! ^ 

cl cl ' 

P 2 - ^±^2 { -wC l+n^t^/d)} w /d! + n^b^ exp{ -wC l+n^b^/d)} w /'d I .. 

and 

Pj = nj_(b2+b^) exp {-w(l+nj^tj_/d)} v/Vdl« 

Simplifications for J^/s (j = 1/2, .■>»,7) aro similar to that 
oi Ij^^s of Ti.oorv^m Tlieir sjinplifxv.a forms aro as 

follows : 

b^ {I~(n 2 +n 3 )c/d}'"^" 

^2 ~ ^1 

J 3 := [ljp{ C 2 I (n 2 +n 3 )c/d} ~ 1 + 

=; exp { (n 2 -tn 3 )a 2 "l} Qd^^ 2 *^^ 

J 5 =b^ [oxp(-n 3 _a 2 ){Q,^(c 2 l 0 )-Q^(c 3 l 0 )} ~ {Qa(c 2 l n^c/d) 

Q^(c 3 in^c/d)}]+b 3 [ L^(c 3 !n 3 c/d) 

-L^C c 2 1 n 3 c/d )-exp ( n 3 a 2 ){ Q^( c 2 1 0 3 !°^^ 
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” ^2^ ®^(“'n^cX2)-ejp(~n^a2 )} Q^( Cji0)+b2{ exp(n2Ct2^“®xp(’^C(;2} 

J2=b4 [exp(~nj_a3)Q^(c3iO)-.Q^{c^!n5^c/d}] , 

By combining these expressions, we get the required cdf of 
given in equation (5»’^*4).« 


For unequal sample sizes, the non-null distribution 
of is very complicated* It can be obtained either by proceeding 
as in Theorem 5*3*1 or by the method given by Khatri (1974)* 

For equal sample size case n^ = = 03 = n and 0j_>62 — 

the non—n\all distribution of T 3 obtained in a similar manner or 
deriv, <,l from Khatri'^s result is given by 


(5.3*5) P [T3 < c|0) ]= Q^(c 3 t 0 ) - 2 g^£ (l+G/d)‘'‘^-(l+ 2 c/d)’''^} 

~ ( 96+^7 ^^d^ ^ 2 ’ 2 c/d) 

- ( ) Q^( C 3 1 c/d) + 2 g 3 Q^( C 3 1 2 c/d) 
-(g^+g 2 ^^d^'=^ 2 *'^/^^'^^% ^^(c^l 2c/d) 

+( g 2-g;j_ ) L^( C 3 1 o/d ) . G > 0,. 


where Cj 


dTj/c, rj = n(8^-6j)/h (j - 2 / 3 )/ 


= oxp(-r 2 ~r 3 )/ 6 / 92 = ejq) {'~r^)/2, 93 - exp( 273 - r 2 ) 76 , 

g^ = exp(7 3 )/ 2 / 95 = oocp 2, gg = exp 2 , 

97 = exp (2r2-.r3)/6/ and d, Qa( -i •) and L^( *1 *) are given b. 
oqua-tion ( 3 * 1 * 2 )* 

One of the objectives of this chapter is to study the 
performnee of the statistic q for different combinations of 
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and d along with a comparative study of the statistics 
" 1 ^ 3 / U 3 and for n^ = ^2 “ ^^ 3 " Methods of obtaining the 
rcquirvid critical points for all Laoc procedures are 

given in Ihe next section* 


5*4* The critical points of the best statistics « 

From equatjon (5*2*8)/ the upper 100a percent critical 
point c. ^ of the test statistic is qiven by 

^ X / Ou J. *' 

h - !t 3^' mrlp ' " 1 /'" i 1 -“ 


( 5 *^- -1 ) c 


l/OC 


aN ^ n. 


/ n-x/N < 


whore a is the chosen level of significance# Some critical points 
of are tabulated in Table 5 *4. 'I* 

The critical point c^ ^ of the test statistic is 
given by P [T 3 > 03 ^q^ 1 Hq] = a* Prom equation (5*2^11)/ it is 
clear that, 03 ^^ is the solution of the equation 


(5 *4 *2) 


(?Txn i+( 3+1^03 

3=0 


a => 


As WO had monti/.'.ned in earlier soctiors/ -che statiota.es T.;^ and 

arc ver-/ closely related bo each o-tl.or* Although IChatri (1974) 
has provided a table for the criticea points of U 2 / he has not 
mentioned the procedure for solving oquat-ion (5*4»2)* Hero wc 
suggest Nuvhxon-Raphson method for obtaining 03 ^^ by solving the ^ ^ 
equation (5*4*2) with an approximate critical point as the a.nieiaj 
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valuo# Note that^ for d not too small^ {l+(j+l)c 

3 /CX 

decrca.jeo rapidly as j increases • Hence, the approximate critical 
point may bo taken as the solution of the equation (5*4*2) 

cr>rrerjpondincj to the term j = 0. Gonsequently, 


(5 *4 *3) = df ^11 

3,a - 1 } ^ 

It is cl i;ar from the equation (5#4<»2), that for K = 2 ^ c? given 

3 /CC 

in oqiiatien (5*4*3) is the exact critical point. 

Some exact and approximate critical points are tabulated 
in 'I'able 5*4*2 for a = 0*05 and for some selected values of K 
one! d* Frc'jm the extensive study made in this direction/ we 
conclu<lr* that, except for large K and very small d, the approxi”- 
mition is reasonably good* Note that, is equal to the 

critical point of U 2 tabulated by Khatri (1974), 

For the statistic U^, Khatri (1974) has tabulated the 
critic a no inns for a = 0,05, 0,01 and for some selected values 
oL nj_/h 2 /n^ and d* For equal sample size case, the critical 
point in the solution of 


(5 *4 *4 ) 


K~2 
E {■ 
j=0 


l)'^( 


K -2 
i- S (■ 
j =0 


• 1 ) 


^ (^TJ ^ ^ l+u+ 2 ju) 

J +1 


~d 


= Kx* 


As tor Tj, the solution of this equation for terms corresponding 
t^) j = 0 / namely 


(K-i)(K+2)(1+u)*"'^ - 2Ka 


{ 5 *4 #S ) 
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gives a fairly good approxima-tion to IhLi exact critical point 
for large values of d* This could bj used as an initj.al valuc- 
'/lug uguaLion (5^4 4 /iJ* TabJL 6 i yiv’es exacL and 

appr'.ximicc critical points of U 3 for 0 : ■-■= 0:^05 and selected 
valuv-s of K and d» 

As Khntri (l974) and Singh (l96S) have shown that 
h.as an ®' 2 (K“’ 1 ) critical points of can 

bo • .btained easily* 

Gome critical points of these test statistics T^^T^^ 
U 3 and arc tabulated in Table 5 *4- *1 for K •= 3/ n^=n 2 =n 3 =n/ 
a = 0 *05 and a = 0 *10 • These are used for studying the power 
functi' n those tests* 


5 *5 • Power of the tests* 


The pov/er of the test Tj_ is P ^ 

C-, ,, is the exxet criticod point given in equation (5^4^l)>^ 
By making use of the non-null distribution of T^ given by 
equation (5 some power values of 6^ > 02 t 

evaluainxl o.nd tabulated in Tables 5 "5 d,- 5 -5 -<2 and 5 -5 -^3 xor 


a = 0*05 and different coiabinations of 

Table 5*5*1 shows that^ for fixed 
d, the p -wer :■£ the tost increases very rap 
But this is not thu case for changes in n 2 r 
as is soon in Tables 5*5*2 and S'B-'S* 


_ . n, and d 
2' o 

values of ^ 2 -^ ^3 
r'dly as iXj increases* 
n^ and d for fixed 
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pOWGJT of 'T' *4*-. i*j_ • 

tne test is obtained by using 
^c^ation (5 .,. 5 ). si.co the derivation of the power functiona 
of and T 3 Involved lengthy calculations, the slinulated 

power v.-auoo of those tests are also tabulated along with the 

exact values in Table .c; .a - 

lauxe t.b.4. This serves as a check for 

t h eo Ti !t ic* 1 } express ions <» 


The power functions expressions of the tests and 
provided by Khutri (l974) are extremely coirplicated even for 
”1 = ^2 “ ”3 Hence, only the simulated power values are 

ol)talntx3 l)y using 1000 iterations# Some of these values are 
tabul.'iteci in Table 5 #5 #4# The nature of Tj^ is laritirely different 
from thv. remaining statistics* It can be used for testing the 
imfu.ility of (i = 1, 2/*»ii,K) against a specified alternative 

^ (0 2 JtQ -^0 ••gSy^) (with suitable relabelling if necessary)* 

It can bo seen from Table 5*5*4, that the power of T^ is 
con side r.ibly higher than that of other three statistics* 


Since the power values for all the tests are provided 
only for 0 ^ ©2 — difficult to compare the 

perform jict., of those statistics* However, some conclusions 
can bo drawn from Table 5*5*4* Both tests T^ and perform 
cquiliy well* Thu* test II 3 performs slightly better than 
when large, while T 3 performs better than U 3 if 

is small* If small, then is better than T^ 

while for largo, the reverse is the case* Similar 

conclusions are expected for other values of K* The LR test 
statistic is recomraaaded for testing against a general 
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alternative hypothesis, since its critical points are easy to 
ov-iluatfi from the F~distribution even for unequal sairple sizes 
Aaaiin;t a spciciFiad alternative like 6^^ > max (© 2/^3/ " ^ 
statistic is racoinm ended • 


the 



TABLE 5 '»4 •! * Exact critical points of the tests ^3 

and for K '= 3 and rij^ ^ 'n 2 = ' n^' '= 


— 

Tests 

n = 11, 

d = 30 

n = 21, 

d = 60 

a = 0-05 

a = o»io 

CC = 0 #05 

a, = o.ic 


0 -2459 

0 .1780 

0-1261 

0.0913 j 

T 3 

3.9036 

3*1125 

1 3.7878 

3.0402 i 

1 

^3 

t 

0.1236 

0.0977 

0 .0600 

0 .0477 1 

i 

^4 

2.52 

2.04 

2.45 

1 *99 

{ 

1 

i 
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table S-4.2. 


Exact critical point c^ q, of in top row 

and its approximated value c* as given in 
equation (5»4»3) in bottom row for a = 0«05 



6 .9443 
C .9443 

5 .1433 
IS. 1433 

4.4590 

l.*1590 

4 *1028 
4.10 28 


9 *6651 
10 .6491 

6 .8972 
7.2599 

5 .8589 

6 *0595 

5 .3 222 
5 *4564 


11.4509 

13.4919 

8.0316 
8 -7446 

6.75 33 
7.13 26 

6 .0937 
6 .3397 


12.7893 
15 .8885 

8 <^786 
9 .9 266 

7.4181 
7*96 28 

6 .6648 
7.0112 


13.8618 
18 .0000 

9 .5566 
10 .9 248 

7*9494 
8 .649 1 

7.1203 

7.5594 


15 .5 267 
21.6643 

10 *6096 
12.5775 

8 .7736 
9 .759 2 

7.8 258 
8 .4337 


6 

3*885 3 

3 .885 3 

4 .9959 

5 .0959 

5.69 31 
5.8716 

6 .2073 

6 .4547 

6 .6166 

6 *9 266 

7 . 2497 
7*67 22 

7 

3. .7 389 

3 .7 389 

4 .7770 
4*8567 

5 .4 245 
5.5638 

5 .9006 

6 *0909 

6 .2789 

6 .5 149 

6 -*8631 

7 .1804 

8 

3.6337 

3*6337 

4-6 201 
4.6867 

5 .2321 

5 .346 2 

5 .6809 
5-8349 

6 .0369 
6.226 2 

6.5860 

6 .8373 

9 

3*5546 

3*5546 

4.5022 

4.5597 

5 .0876 

5 .1845 

5 .5159 

5 .6453 

5 .855 2 

6 .0129 

6 .3777 
6.5348 

10 

3 .49 28 

3 .49 28 

4 .4104 

4 -4613 

4.9751 ' 
5 .0597 

5 .3875 

5 .499 2 

5.7137 

5 V3439 

6.2156 

6.3913 

15 

3.3153 

3.3158 

4 .1479 

4 ♦IS 21 

4-6537 

4 .7076 

5 .0205 

5 .0393 

5 .•='094 
5.^90^ 

5 .753^8 
5 ^5 29 

2fJ 

3.2317 

3*2317 

4.0 235 
4-0510 

4 .5016 
4.5436 

4 *8469 

4 .3991 

5 .1131 
5-1735 

5 -6322 
5 .6057 

50 

3.08 7 3 
3. OB 7 3 

3.3106 
3-8 234 

4.2415 

4 .2667 

4.5501 

4.5793 

4.7909 
4.8 239 

5 -1565 
5 .1941 

100 

3.0411 

3 .0411 

3.74 27 

3 .7573 

4 .1585 
4.1793 

4.4554 

4 .4795 

4 .6367 
4*7129 

5 .0363 
5 .0653 


16 *8016 
24 .8 328 

11.4168 
13.9 386 

9 .4053 
10 .65 14 

8 .366 2 

9 .1 26 2 

7 .7340 

8 .2 57 

7 . 3095 ^ 
7 .6987 

7 .0050 
7 .3108 

6. 7761 j 
7.0 261 

6 .5978 
6 .8084 

6.0874 
6 . 205 2 


S>^455 . 
5 .9 295 1 

I 

5 .4314 
5 .47 22 

5.299 3 1 
5 *3302 
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TAiUih; Exact critical point of in top row and. its 




approximated value 

aiven by ecru at ion 

"('5 . 4 - 5T 




in bottom 

row for a 

= 0.05. 












x: ■" 1 

N'' 

Cl 

2 

3 

4 

5 

6 

8 

10 

_ ... 2s - 

2 

3-^7 21 

3 *47 21 

4*5099 

4.7736 

5 .1599 

5 *708 2 

5*1317 
6*48 33 

5 * 26 22 
7*1650 

6 .0087 

8 -3541 

6 *2515 

9 .39 23 

3 

1 *7144 

1 *7144 

2*1559 
2*218 3 

2*4305 

2*5569 

2*4610 
2*8 259 

2*5454 

3.0548 

2-8057 

3.4395 

2-99 38 

3 #76 22 . 

4 

1*1147 

1*1147 

1.3776 

1.4028 

1*5399 

1*5900 

1*5687 

1*7356 

1 .6 260 
1*8574 

1 .7717 
2-0584 

1.89 33 
2-2237 

5 

0 *8 206 
0-B205 

1*0031 

1*0164 

1*115 2 
1-1411 

1*1390 

1*2368 

1*1812 
1*316 2 

1-2803 

1 -4457 

1*3661 

1-5508 

6 ^ 

0-6475 

0 -6475 

0-7858 

0*7940 

0*8702 

0.8860 

0 *8900 
0.9560 

0 .9 230 
1*0136 

0.9973 

1 *1070 

1 .06 25 

1 -18 22 

7 

0*5 341 
0-5 341 

0 .6447 
0*6503 

0*7120 
0*7 226 

0*7 287 

0 .777 2 

0*7556 
0*8 220 

0^147 

0-6942 

0 #3668 i 
0 .95 20 1 

B 

0 *454 2 
0.454 2 

0.5461 

0.5501 

0 *6013 
0.6093 

0*6161 

0 -65 ^9 

0 6339 

0 *6904 

0 -6377 

0 .7433 

0.7303 

0.7955 

9 

0.3949 
0 #39 49 

0.4733 

0.4764 

0.5 207 
0.5 265 

0#5333 

0-6640 

0.55 23 

0 -5946 

0 #59 44 
0.6435 

0.6311 

0 .63 24 


10 

15 

20 

50 

100 


Op349 3 
O #349 3 

0*2211 

0#2211 

0*1616 

0*1616 

0*0617 

0*0617 

0 .0304 
0 *0 304 


0 .4175 
0*4 200 

0*2623 

0.2633 

0 *19 10 
0 .19 16 

0*0725 

0.0726 

0 .0356 
0*0357 


0*4537 
0 *4633 

0*369 

0.2039 

0 - 2005 
0*2096 

0 .0730 
0 .079 1 

0.0337 
0 *0333 


0 *4699 
0.4956 


0 .4369 
0 .5 219 


0 .5 231 
0*5639 


0.5549 

0-5971 


0.2939 

0 .3073 

0 .3044 

0 .3 231 

0-3260 

0-3473 

0 345 0 1 
0-3663 ‘ 

0*2136 

0.2229 

0*2211 

0*2337 

0 -2364 
0-2505 

0 .2499 
0-2633 

0-0307 

0-0838 

0 *03 35 
0.0376 

0-0390 

0-0935 

0 .09 33 
0.09 3 2 

0-0396 

0*0411 

0-0409 

0.0429 

0*0436 

0.0457 

0.0459 

0*0479 



TABLK Exact power of the test for a = 0.05/ 0 ^= 20 / 

n ^=15 / d Iz'To'', K=:3~' and 0'^=O , ' 


25 

30 

.0500 

.0500 

.1745 

.2240 

.5561 

.7 207 

.9013 

.9497 

.98 27 

.9913 

.9970 

.9985 

,9995 

.9997 

.1206 

.1477 

.3949 

.5180 

.7715 

.8663 1 

\ 

.9533 

.976 2 

.9918 

.9959 

.9986 

.999 3 , 

.9998 

.0999 

.3694 

.48 20 

.6954 

.7706 

.8919 

.9368 

.9779 

.9888 

.9961 

.9980 

,999 3 

.9997 

.9999 

,9999 

.8503 

.83 36 

-9320 

.9488 ! 

.9759 

.9859 i 

.9951 

.9975 

,9991 

.9996 

.9999 

.9999 

.9666 

.9740 

.9843 

.9886 

.9946 

.9969 

.9989 

.9994 

.99 25 

.9942 

.9983 

.9987 

.9996 

.9997 

.9999 

.9999 


0 



10 


15 


20 


. 0 (> 

. or . 

.0500 

.05 

.00 

.0641 

• 10 

.( X ) 

.08 24 

• IS 

•oo 

#1058 

.20 

.00 

•1359 

.25 

.00 

.1745 

.30 

.00 

.2241 

.05 

.05 

.0584 

-10 

.05 

.0749 

1*15 

.05 

.096 3 

• 20 


.1236 

.25 

.05 

• 1588 

|.30 

.05 

• 2038 

.35 


.2616 

flO 

*10 

.07 23 

• 15 

.10 

.09 28 

♦ 20 

# to 

.1191 

-25 

.10 

.1530 

• 30 

.10 

.1964 

.35 

.10 

.25 21 

.40 

-10 

•3 228 

• 20 

.20 

.1169 

.25 

.2 CJ 

.1500 

.30 

.20 

♦ 19 27 

.^5 

.20 

.247^ 

.40 

.20 

.3166 

.45 

.20 

.402^ 

.30 

.30 

-19 21 

.35 

*30 

.2466 

.40 

♦30 

.3153 

.45 

.30 

.4012 

-40 

.40 

.3157 

•SO 

•50 

.5006 

•60 

.60 

.7312 

.70 

.70 

.87 23 


.0500 
.03 24 
.1359 
.2241 
.3674 

.57 » 
.78 23 
.0690 
*1137 
.1875 

.3081 
.4891 
.699 2 
.8660 
.1074 

.1771 
.2909 
.4611 
.6596 
.8 268 

.9301 

.S45 

.4506 

.6425 

.8018 

.90 27 
.957 2 
,.6411 
.7995 
i 69 S 9 

.9517 

.3986 

.9768 

.9948 

.9988 


.0500 
.1059 
.2241 
.46 24 
.7696 

.9409 
.9886 
.08 24 
.1743 
• 3630 

.6474 
.876 3 
.9712 
.9946 
.1632 

.3395 

.6004 

*8136 

.9407 

.9864 

.9974 
.5868 
.7912 
.90 29 
.9591 

.9867 
.9970 
,8999 
.95 29 
.9733 

.9909 
.9777 
.99 50 
.9989 
.9998 


,0500 

.1359 

.3663 

. 754 ^ 

.9490 

.9910 
.9984 
.099 2 
.268 3 
.5989 

.8769 
.9758 
.9958 
.999 3 
.2509 

.55 26 
.8035 
.9417 
-9886 
,9980 

.9997 

.7778 

.8968 

.9561 

.9870 

.9975 

.9996 

.9504 

.9770 

.9902 

.9971 
.9889 
.9975 
♦999 4 
.9999 


c / 
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rAVibB b«S»2* 


Exact power of the test for a = 0,05, 
^ ~ 15 f d ss 22/ K =: 3 and 0^ =0* 


«1 "2 \ 
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!3*5»2 Confcdi 



00 

.0500 

.0500 

.0500 

.0500 

.0500 

.0500 

OCl 

• 1053 

•1058 

•1058 

.1059 

.1059 

.1059 

00 

. 22'’0 

• 2240 

.2240 

.2240 

.2240 

.2240 

( X ) 

• :550 

*4558 

.4563 

.4567 

.4570 

.457 2 

flf ) 

•7434 

•7474 

•7502 

*75 22 

.7536 

•7547 

Of) 

*9 219 

•9 269 

.9303 

.93 26 

.9343 

.9354 


•98 26 

•9854 

•9871 

•9883 

.9890 

-.9896 

(/) 

•9967 

.9977 

*998 2 

.9985 

.9987 

.9989 

(.(: 

•9994 

.9997 

.9998 

.9999 

.9999 

.9999 

10 

.0395 

•1315 

•1531 

.168 2 

•1790 

*1371 

10 

.2070 

• 2122 

.316 2 

.3469 

.3690 

.3853 

10 

.^873 

-4914 

.5603 

.6074 

.6406 

•6645 

10 

.6190 

.7209 

.7858 

.8 283 

.8570 

.8768 

lO 

.834 2 

.8903 

•9 250 

» 9 d 6 3 

•9593 

.9637 

10 

*9514 

•9711 

.98 21 

.988 3 

.99 21 

.9944 

ao 

*9394 

•9945 

.9969 

.993 2 

.9989 

.999 3 

10 1 

.9930 

.9991 

.9996 

•9993 

.9999 

.9999 

20 ' 

•3118 

•457 3 

•5444 

.5998 

.6 369 

•66 27 

'20 

•4685 

.6534 

.7545 

•8135 

•8498 

.8732 

■20 

.618 2 

•7966 

♦88 26 

.9 267 

.9505 

.964 2 

■ 20 

.7677 

.8954 

.9498 

.974 2 

.9853 

•9915 

■20 

.6994 

.9597 

.9831 

.99 25 

.9965 

.998 2 

■20 

•9705 

.9894 

.9960 

.9984 

.9993 

.9997 

■30 

-57 24 

.7841 

.8791 

.9 257 

.9502 

.9641 

•30 

*6765 

.8705 

.9428 

.97 22 

.985 2 

•9913 

■30 

•768 3 

.9 250 

.9736 

.9899 

.9957 

.9980 

■ 30 

•8591 

.9615 

.9888 

.9965 

.9988 

•9996 

w 

fiJyO 

-7406 

.9 204 

.97 3 

.9897 

.9957 

•9980 



TA!’,LK 5.5-3 


i^act power o f the test for a = 0-05, 

= 30, = 15, K = 3 and 9^ = 0» 


1 






O2 \, 1 

5 

10 

J )0 

•fX) 

-0500 

•0500 

p()S 


• 2129 

• 2213 

»l(i 

-00 

- 540 ^ 

•6 347 

*15 


-618 2 

•9074 

*20 


: .9438 

•98 23 

• 25 

- 0 (.i 

-9854 

.9969 

*30 


•9966 

•9995 

•10 

•10 

: -3773 

.4309 

-15 

•10 

•6486 

•7 250 

•20 

-if.» 

•8465 

•9064 

• 25 

•10 

•9475 

•9191 

.30 

•10 

•935 2 

•9960 

• 35 

•lO 

•9964 

•999 3 

•40 

•10 

i .999 2 

-9999 

• 20 

• 20 

•8088 

•8609 

• 25 

• 20 

.0007 

.9380 

• M) 

i 

- 2.0 1 

.96 35 

•9791 

• 35 

•20 

•9879 

.996 2 

•40 

• 2 fi 

•9966 

•9991 

•45 

• 20 

•999 2 

•9998 

.30 

•30 

•9551 

•9689 

.36 

•30 

S 19 S 

•986 2 

.40 

*30 

• 99 18 

• 095 ^ 

•45 

.30 

- 997 ^ 

•9989 

•40 

•40 

•9899 

.9931 

.50 

• 5 C) 

•9977 

•9985 

•60 

•60 

•9995 

.9997 

•70 

•70 

•9999 

.9999 


15 

20 

25 

30 

.0500 

.0500 

.0500 

.0500 

.22'='1 

• 22^1 

.2239 

.2240 

•6736 

.6959 

.7104 

•1201 

-9317 

.9416 

.9467 

.9497 

-9880 

.9898 

.9907 

-9913 

.9979 

.998 2 

.9984 

•9985 

-9996 

•9991 

.9997 

.9997 

.4541 

.4673 

.4759 

•48 20 

•749 2 

•7603 

.7666 

•7706 

.9 224 

.9 297 

.9340 

•9 368 

.9847 

.9870 

-9881 

•9888 

.9973 

.9977 

.9979 

.9980 

*9995 

.9996 

.9996 

•9991 

.9999 

.9999 

.9999 

.9999 

.8734 

.8787 

.8317 

.88 36 

.9440 

.9465 

.9479 

.9488 

.98 27 

.9843 

.9853 

.98 59 

.9966 

.9971 

.997 3 

.9975 

.9994 

.9995 

.9995 

.9996 

.9999 

.9999 

-9999 

.9999 

.9717 

SI 29 

.9736 

• .9740 

.9875 

.9381 

.9884 

•9886 

.9961 

.9965 

.9967 

.9969 

.999 2 

.9994 

.9994 

•9994 

.99 ^7 

.9940 

.9941 

.994 2 

.9986 

♦9987 

.9987 

•9987 

•9997 

.9997 

•9991 

-9997 

.9999 

.9999 

.9999 

.9999 


TABI/K 5 *5 *4 - Exact and simulated powers of the tests and 
and simulated pov/ers of the tests U_ and fot 
a "= "o’.05, = 0 and n^=n 2 =n 2 =n» 




I m()C) 

-05 

i-lO 
; -15 
1 - 20 


® 2 

T 

1 

.00 

.050 

• 

o 

o 

.037 

)»0C 

•ISO 

.00 

j .260 

6 

c 

.441 


1 *25 AX) 

I *•^0 *00 

U^S A3C> 

1*40 -00 

*50 A )0 


.669 
.850 
*9 44 
.981 
-997 


10 

.10 

15 

-10 

'20 

.10 

-25 

.10 

-30 

-10 


•117 

.202 

-541 

.731 


.35 

.10 

.40 

•ic 

.50 

•10 

.20 

^20 

! .25 



-^70 
*947 
.99 3 
-334 
.5 21 


.30 

-35 

•40 

.50 

.30 


.20 I -699 


*20 

.20 

.30 


-8 27 
.008 
-98 2 
.696 


(n = 11/ d = 30) 


t: 


T. 


T 


* 



.066 

.050 

.099 

.056 

.147 

.076 

.271 

-111 

.469 

.173 

*674 

.279 

•846 

.439 

.940 

.6 29 

.980 

.795 

.998 

.961 

-112 

•08 2 

.206 

•107 

.342 

•153 

-518 

.2^2 

.734 

•358 


.05 3 

.049 

.054 

.058 

.078 

.075 

.128 

.129 

.200 

-211 

.286 

.319 

.439 

.49 2 

.6 29 

-65 3 

.733 

.8 29 

.960 

.97 2 

*07 2 

.063 

•109 

.100 

.149 

*141 

• 247 

• 254 

-371 

.411 


.05 3 
.064 
•03 2 
.156 
*226 

• 266 
.39 2 
.601 
.654 
.890 

.115 

*165 

.236 

.319 

.461 


.869 

.947 

.993 

.333 

-516 

.688 
.8 25 
.906 
.934 
.673 


-5 22 
.691 
.917 
.207 
.271 


.5 2S 

.667 

^17 

.225 

.277 


-5 29 
-731 

.940 

.169 

.236 


.5 30 
.660 
.358 
.4 29 
.5 26 



^Simulated powers based on 


1000 sanvples for each sanple size 
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TADL.K 5 *5 »4 • Cont.d* 





(n = 

It 

60) 



®1 

0 

T 

1 


T 

3 

3 

u* ■ 


-oo 

.00 

•050 

.042 

.050 

.048 

.051 

.043 

-05 

.00 

•143 

•135 

.073 

*064 

#066 

•07 2 

-10 

.00 

•407 

•416 

-160 

•174 

.175 

•18 3 

*15 

♦00 

.849 

•853 

-403 

*39 2 

•433 

• 360 

-20 

-00 

•981 

.980 

-804 

.810 

.866 

.660 

-25 

-00 

.998 

.999 

-969 

.976 

•989 

•906 

-05 

-05 

•112 

•111 

-079 

.07 3 

•075 

•09 Q 

•10 

-05 

•318 

.319 

-143 

.154 

.157 

•180 

• 15 

.05 

•111 

.716 

.3 28 

.3 29 

•328 

•335 

-20 

•05 

•946 

-947 

.68 2 

.675 

.679 

.670 

• 25 

.05 

•99 3 

#99 3 

.9 24 

.9 27 

.9 32 

.9 23 

•10 

•10 

•307 

.313 

.19 3 

•199 

.146 

*400 

-15 

•10 ! 

•636 

•686 

.351 

•349 

•347 

.6 28 i 

• 20 

•10 

•900 

.913 

•661 

.661 

.699 

-783 j 

-25 

.10 

•981 

.986 

•88 2 

•890 

.908 

^98 i 

• ^0 

•10 

-998 

.997 

.974 

.98 2 

.987 

.958 

• 20 

-20 

-888 

.864 

•756 

•749 

•701 

•898 

• 25 

• 20 

•961 

.960 

•88^ 

.88 2 

•839 

♦944 ! 

-■='0 

• 20 

•983 

•935 

.957 

.950 

.956 

♦966 

• 30 

•30 

j 

! 

*986 

.939 

*969 

•975 

•974 

♦99 3 


■^Simulated powers based on 1000 samples for each sample size- 


CHAPTER VI 


GHijEUAIilZaJ STATISTICS FOR THE EQUAL SAMPLE 
CASE U'iiEN OM2 OBSERVATION IS MISSING ON THE LEFT 


r. . t « Introduction and test statistics * 

In Chapter V, tests for the equality of location 
parameters of K (> 3) populations are discussed, when the 
smallest observation is available in each sample. Here the 
prob.1 eni is extended to the case when the smallest observation 
in rnissinq tut the second smallest observation is available 
in ouch sample of size n, that is. 


(i = 1,2, 

arc the available observations from the i^ population with 


( n—s^) ^ 2 “ 


Similar to 

Jl) 


of Chapter V, the test V^ defined by 


V. = {X;-"' ~ 

1 2 2<i<K 




= 6 against 


( 6 .1 -l) 

ju 

is proposed for testing s - ^2 " * 

H : 6<, > max where 

^ ^ 2<i<K- 

‘i (n-s-2). 


K 


(6.1.2) do* = 

i=l 3=2 


Tht- gene 


ralization 


of statistics T 2 


Chapter V are ^2 "^3 


.spectively. These for ogual sample 
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i'ir'A.. caS'j are given by 


( *1 -3 ) 


'lur’ 


V 2 = { max - min 

l<i<K ^ l<i<K 


( 5 ♦ 1 »^* ) 


V3 = [ max { (X^^J^-X^„^b}]/(7-^. 
2<j<K ^ ^ 2 2 


^ o i:h ol t,hu30 are proposed for testing against s atleast 
')n<, Cj j.s different from 0. Compared to two-sample case (K = 2), 
the iiR tf'St is much more complicated even for K = 3* Even the 
') riv. ifion of ML estimates under H is far from simple and 
a full study in itsolf» Consequently, we have not 
otudjod tho LR test in this case and have loft it as an open 


p ro bl (i-m * 


The test procedure is to reject if > c^f whore 
fn*- constants c^^ are obtained from solving the equations 


(6.1,5) P [Vj_ > CjLlH^] = a, 

'•/bore (X is the chosen lavol of significance* The required null 
distribation of those statistics are derived in the next section. 


d . 2 * Distribution theory . 

Lumma 6.2.1. Lot Y- (i = 1 , 2 , ...,K) bo K i-i.d. random variates 
'v/ith pdf 

(6.2.1) f(y) = n(n-l) [exp{--(n~l)y}-oxq)(-ny)], y > 0 , 

then Z = {Y^-min(Y 2 /Y 3 , ...,Yj^)} has the pdf given by 



130 


( 6 -2-2) 


wj lOre 


f ( z) 


f;J_(z), Z < 0 

f 2^ 2 ) , z > 0^ 


K~2 


= (K-i) s {® 


34=2 35Z 


j =0 3 '•3, -Jo Jo'J 


r-1/ 


K-2 


'1 -'2 -^2 -"3 


^-'(n-l)z 


io(z) = (K-^l) s' (~i)3(^7^')n^"’J'(n~l)J'^2 ^ ^ , 

j=0 3 


^1-^2 


^2-^3 


■jl rr (n*-l)K+j, = j^+ 1 , 33 = j^+ 2 , = (n~l)(K-.l)+j 

'iod + 1, 

Pr oof s From eqi.iation (6.2*1), the cdf of is 

(6.2*3) F(y) = [ !+( o-l)exp(~ny)-n e 3 <p {~(n-l)y}]< y ^ 0 
The pdf of Z2 = min ( > » .,Yj^) is then given by 

£(^2^ = (K-l) [ne 

~ ( n— 1 ) ( K— 1 ) ! 7 . ^ _ „ „ T . , 

'(l~e ^) [l~(^^)e 


(n-l)z„ -nz^ „ „ --(n~l)z, ~nz. 

2_(n-.i)e 2]K-.2 n(n^i){-e 2„^ 

. ~(n-l)(K-l)z„ ~Z„ r.~1 

(K-l)n^-l(n-l)- * 


e 


Expanding { ( l~(n--l)exp(-z„)/n}^"*^ as a binomial sum , we get 


2 J=0 J 

* exp [ --{ (n~l ) CK -1 )+j} Z2 ] * 
Now, tho jpdf of Z2 and Z^ = Y^ is 
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K-2 


(6.2.4) fiz^.z ) = (K~l) S '(~l)J(^:2)nK-J(n-i)J‘+2(i^e 2) 

j =0 3 


.exp [~{(n~l)(K-i)4-j}z ] 


* [G5<p{~(n~l)z^}~e3{p(-nz^)] 




r-U.ko a transformation z = ~ ^ 2 ' range of the 

hransformod variables arc z^ > max (0,--z), -«. < z < oo, 
'rhun.tho marginal pdf of Z is 


( 6 . 2 .*5) f(z) 


OO 

/ £iz 2 , Z+Z 2 ) 6 Z 2 / z < 0 


CO 

J* f ( ^ 2 ^ ^"^^ 2 ^ dz 2 / z 6 y 

.. o 

where f( ./ .) is given as in equation (6.2»4)» On simplification/ 
the equation (6-2.5) givc-s the required density function of Z. 


Lemma 6.2*2. Let (i = 1/2/.../K) be K i.i.d. random variates 
v;it}i pdf given as in equation ( 6 . 2 * 1 )/ then the pdf of 


Z = max (Y. ) - min (Y.) is 
1<1^ ^ l<i:^ 


\K~i~1/K“-l\ i+ 1 / >. < *'( n~l ) z, j—l 

(6.2.6) f(z) = K S (~l) ^ ^( . )n^^^(n-l) { 1-e } 

j=0 


- ( nK-n- j ) ^ / z > 


0 . 


Proof. The cdf of Z is (David/ 1981, p- 26) 


G(z) = K X f(y) {F(y+z) - P(y)}^ ^ 


( 6 .2.7) 
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= (K-1) S ^ _-(a+l 

^ «siO ^ •!» 




i-jCvj) = (K-l) s’ .^,-(a+l)^^ ) 


- (l+nVj^/ai-<«+lV(3233>] - 

and (i = l,..,^ 5 ) are given in equation (6'.2»2)« 
iIiL'.ool- Note that, - e^)/a (i = i, 2 , ,..,k) has the 

d ifjLrijDut ion given as in equation (6*2*1) • Hence, under 
= dZ/W, where 2W = 2da*/a has a distribution and 

-'I ;z - min (Y^, Y^, • ■••,Yj^) has the pdf given in equation (6*2*2) - 
Vy writing the jpdf of Z and W, and making a transformation 

™ da/w, w = w, and integrating w*r» to w, we get the required 
piii of given in equation (6-2*11) ■» 

For K = 2, equation (6*2*11) simplifies to 

: [n{ l-(n-l)v^/d}"'‘^“^~(n~l){ l-nvj^/d} ^ ^], < O 

f{v^) = 

: [n{l+(n-l)v^/ar‘^"^-(n-l){l+nv/d}-<^~h, > o 

v/hich agrees xirith the null distribution of T, given in equation 
( 3 * 2 - 7 ) for = r^ = 1 and n^^ = n 2 = n case* 

For K = the null distribution of is given by 
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I" 2(n~l )v^/d} l-( 2n~l )v^/d} 

(f-:i-12) t(v^)^ ^ B^{l-2nv^/d}“‘^“'l] V < 0 

[®5^ 1-+^n~l)vj_/d}“*^'*^-Bg{ n*nv^/d}"''^''^],Vj_ > 0, 

whon; = 2n(n~l)/{3(3n-l)(3n~2)}, = n^On-^l), 

^^3 = 3n(n-.l)(2n--i), B. = (n-l) 3n~2) , 




= n(5n~3) and Bg = (n~l) (5n-2) * 
it tr. easy to establish the following relations among 
^2 ~ ^4 ~ ®5 "* ®6 “ 4 n - 2 , 


B, 


B- 


B. 


•. r I ^ "^4, 1 

*"1 ^ 2n--2 2n-l 2n^ “ 3 ' 


ir'f.i 


£■ 


B, 


3, 




2 

= I‘ 


1 ^ li-i n 

relations can be used for showing that f(v^) given in 
^i.ion (6-2-12) is continuous at =0 and it is indeed a pdi 

Since the null distribution of statistics V 2 and are 

coifiTvi iciit ed/ only the case K = 3 is considered and discussed in 
th:; iol lowing theorems : 

Thu^ rem 6«2»2» Under the distribution of V 2 for K = 3 is 


(6.2*13) f(v2) = 3B3_ [~B2£l+2(n-l)v/d}~^"*^-i-B3{l+(2n-l)v2/d} 

-B^£ l+2nvy'<i} ^ ^+Bg{ l+(n~l)v2/d} 

- BgCl+nv^/d} / V2 > O, 
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where are given in equation (6#2*12)» 

ProoC- hote that, j_)/cr (i = 1,2,3) are i»i-d* 

raiKlo-:* Vf,iriates with cortmon pdf given in equation (6*2*l)» Now, 
f txw) betcOTa 6*2-2, for K = 3, the pdf of 


Z = rncix (Y^,Y2,Y2) - min (Y^,Y2,Y2) 


e 


an he written as 


= ayw, where w = **/»• the trensformetdon 

V,, =■ dz/w, w = w, and integrating w-r-to w, we obtain tne 
* . 
rewutn.-d wdf of given in equation (6.2*13}. 

■rh,.nr...m 6 . 2 .^. The null cdf of the statistic V, for K = s Is 

-d - ^ .---d 


(6.2-14) P [v, < olH^] = l-lAi{l+nc/a}"'^-A2{l+(n-l)c/ai 

+ A,t l+2nc/d) "‘^-A4t l+( 2 n-l ) c/di “ l+2( n-1 ) ydl 
+ A,£ l+anc/dr'^-A,! l+(3n-l)c/«"*-^8f 3h-2)t/« 


-d 


-a 


+ Ag{ i+ 3 (n-l)c/d} c > 


wlicr^ A^tTCri' 


,)[i.(n-l)/( 2 h-l)-(n-l)V 3 -nV( 3 n- 2 )t 2 n^n-l)/C 3 n-l) 3 , 




[l+y c 2 n-l )-( n-l ) V ( 3n-l )-h V 3t2h( h-l > V C 3n- 2) ] 


4 
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N = 2n(n~l)^ [l/( 2n~l)~rv'(3n-2)+(n-l)/(3n-l)] / 
A^.r2n^(n-l) [l/( 2n-l)-n/(3n-2)+(n~l)/(3n-l) ] , 
ami = [l-2r]/3 + 2(n-l) 3n~2) ] • 

= 1/2,3). Note that, s sre 

;! .j.d* ivinilom variates with pdf fly) and cdf F(y) as given in 
<'qu<ii.io;u5 (6.2*3) and (6.2*3) respectively. 


Let Z = max {(Y^~Y 2 ),(Y^-Y 3 ),(Y 2 -Y^),(Y 3 -y^)} * Thus, 

p [21 < y.\n^] = P [Yj-Y^ < Z/Y^-Y 3 < z, Y^-Y^ < Z/ Y 3 ~y^ ;< z] 

= P [Y 2 > Y^~z, Y 3 > Y^-z, Y 2 < Yj^+z, Y 3 < Y^+z] 

= X P [y^~z< Y2<Y^+z,Y^~z < Y3<Y^+zlY^=y] f(y) dy 

.—00 

CO 

(6 *2*15 ) - S Q( z,y) dy, 

nu.y) = £(y) N [ F (y«)-q(y-a ]. since f(y) = 0 for 
j=2 

V < hence. 


((i, 2 . 16 ) p [Z < zlHol = ' Qa,y)dy + / Q(2,y) dy. 

o ^ 

Not<. that, the first term of the equation (6.2.16) reduces to 
/ f(y){PCy+z))^ dy. Now substituting for f(y) and FCy), we 

O "x 

c..:t the staplifiod form of the equation (6.2*16) as 


(6. 2*17 ) P[Z<z1Hq] 


~nz , -'(n-*l)z 

1+A^o ~A2e +«-3^ 4 

, -2(n-l)z^ Q-*3nz_j^ ~(3n'"l)z 

+ A^g 7 

_j^^^-(3n-2)z^9Q*-3(n-l)z^ z > 0 . 
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Under = d2/W, where 2'^ = 2d0'^y{7 has a X^^-distribution » 

Thus^ 

r oo 

(r..2.aB) P [V3 < CIH^] =; P [z < cW/d!W=w] I » 

o "" 

Sul'jvStituting from equation (6»2»17) and simplifying the resulting 
expression, we get the required null cdf of 

6 *3 * Momr^ts of the statistics under for K =e 3* 

\ -*d-* 1 

Note that /the pdf of contains factors like (.l~av; 
and ( l+av) This allows the evaluation of moments of hy 

using tho following leirma t 

L errina 6 » 3 «i > For a > 0 and 0 < h < d 

(i) dv = (-1)’'^ B(h+l^d-h)/a^‘’'^ 

^C 50 

and 

(IJ) v^d+av)”"^"^ dv = B(h+l,d-h)/a^‘^^. 

o 

I’roof# Making a substitution t = (l~av) / we have 


/ v^( l-av)"*^”^ dv 


mQO 


1 f , ^ ^(a-h~i) 

I (£=i) 

a. 

O 


t „ 


(- 1 ) 




(-1)*’ B(h+l,a~h)/a’^’''^. 


Part (ii) of the lemma can 
V k: -u in part ( i) • 


be nroved by making the substitution 





Now, from Leimma 6»^*1 and equation (6-2*12), the 
■‘ill! trornent of about zero for K s 3 and h < d is given by 

= r3^.B(h+l,d-h)d^‘''^- [(- 1 )^ B 2 /( 2 n- 2 )^‘^^-(-'l)^ B^/ ( 2 n-l 




ivx|jZ'essions for the moments of V ’2 3 -nd can be written down 
i,n a similar manner* 

Mc aairk 6 * 3 * 1 * It is easy to see from the pdf of end ^2 
q ivon in equations (6*2*12) and (6*2*13) respectivelyj that 

£(¥ 3 ) = 3 S(V^). 

This can also be justified from the fact that 


V2 = + VI2 + 

.(i) 


mlB <xa’)}/0* (1 = 1.2.3) ^"'3 «1 = h' 


v/hnro W, « {X^ mxii 

G »4 * Critical points o f the tests for _ K * 3 j; 

The UDper lOCO, percent critical point c^ of the test 

statistic Vi is the solution of ecruation (6,1.5). Prom equation 

(6.2.12), we have (O) = ^ and hence c^ is either the solution 

1 

o£ 

( 6 .4 .1 ) Bi [ ; i+( n-1 ) Cj/di ■'* - -I { l+nCi/di ] - a 


or of 
( 6 *4 *2) 




~d 


B 

_ ,^-l^l-( 2 n~l)c-^/d}' 

X2n-1) ^ 

B/1 «d n 

+ tl-2nci/d} ]=a 
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according as a < 7 / 3 1 or a > 2/3 respectively. From equation 
(6,2,13), C 2 is the solution of equation 


(6,4.3) 3B^ [. ^^i^2(n~l)c/dr'' t C )c/d} 


~d 


^ 4 - ® 

2 ^ {l+2n c/d}’'^ + { l+(n~l)c 2 /d}' 

B. 


— §■ { l+nc 2 /d}~'^] = a. 

‘Similarly, the critical point c^ is obtained from the 
Ofjuation (6-2 -14 ) . 


Somvo critical points of ^3 tabulated in 

Tables 6*4*1, 6*4*2 and 6*4*3 respectively for a = 0*05 and K = 3* 

Those are calculated by using Newton Raphson method with 
approximate critical points as the initial values* Approximate 
values arc obtained by using the fact that the null distribution 
of {V,-E(V,)}/{Var(V. )}^/^' is approximately normal* This gives 

mjm mpIm «pik 

the approximate critical points as 

= cf.{Var(V,.)}^/^ + E(V.) (i = 1,2,3), 

Upp t a* J- X 

wiicro is the upper 100 oc percent of nCO^I)* Some numerical 
calculations show that upper lOOa/2 percent point of N(0,i) 
distribution gives closer approximation for V 2 and than the 

100a percent point* This may be due to the fact that the 
distribution of is in the interval (- 00 , 00 ), whereas the 
distribution of V 2 and Vg are confined to the interval • 
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Some exact and approximate critical points for all the 
three tests are tabulated in Tables 6«4«4 for a = 0»05 and K = 3* 
It is clear from the Table 6*4*4, that the approximate critical 
points are reasonably good for large values of d* 

6 »5 • Perforrnan.ce of the tests* 


It docs not appear simple to evaluate the non— null 
f’iistributions of these statistics# Consequently, we use Monto— 
Carlo techniques for the calculation of power of these tests# 
Sori.j of these vahaes based on 1000 iterations arc tabulated in 
Tabic 6*5*1 for a = 0*05 and K = 3# 


On the basis of these calculations, the statistic is 
rocotntiondod for testing against a specified alternative lilco 
> niiix ( 02 - 63 ) • c°n=lusion is similar to the use of 

01 Chapter V against For testing against the alternative 

Hj, it is otaserv^ad that Vj performs better if is very small, 

otherwise the performance of is better* 

In Table 6*5-2, the power values of the test T^ (of 
Chapter V) and are tabulated for a - 0.05, K = 3, n^=n 2 =n 3 =n =16 
and a = 6(12)42. These calculations show that there is a 
considerable loss of power duo to censoring of the smallest 
observation. This highlights the imoortanco of srmllost 
Observations for tosti^ag the ovality of location parameters. 

,^iT-aTm "hv Greenberg and Sarhan (196 2) 
Similar conclusions wore drawn by Greenne p 

for the estimation of location parameters. 
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TADJjL 6 • 4 , 1 , Eicact critical points of the test 



for a =0.05 and K = 

h. 



nXi 

6 

9 

15 

18 

27 

3n— 6 

4 

1.39 38 






5 

1 .07 21 

0.9866 





6 

0 .07 22 

0,8026 





7 

0.7357 

0-6770 

0-6 340 




8 

0 .6 36 3 

0 .5855 

0.5483 

0.5 394 



9 

0.5607 

0.5159 

0.4831 

0.4753 



10 

0.5012 

0.4612 

0.4319 

0.4 249 



11 

C .45 ^ 2 

0.4170 

0.3905 

0 .3841 

0.3739 


12 1 

0 .4135 

0 .3805 

0 . 35 6 3 

0.3506 

0*M12 

0.3394 

13 

0. 38(^3 

0 .3499 

0*3277 

0 .3 224 

0.3138 

0*3107 

14 

0 .35 20 

0 .3 239 

0 .3033 

0 .2984 

0 -2904 

0.2866 

IS 

0 . 3 277 

0-3015 

0-28 23 

0 .2778 

0.2703 

0-2659 

16 

0 *3065 

0.2B20 

0 .2640 

0.2598 

0.25 28 

0-2480 

17 

0-378 

0 .2648 

0 - 2480 

0 -2440 

0.2375 

0.23 24 

18 

0-2714 

0.2497 

0.2338 

0 .2300 

0 .2239 

0^2187 

19 

0*2567 

0 .2361 

0 .2211 

0.2176 

0.2118 

0.2065 

20 

0 . 2435 

0 .2240 

0 - 2098 

0-2064 

0 .2009 

0 *1955 
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TABLE 6#r»2- Exact critical points of the test V 2 
for (X = 0 *05 and. K = • 


n 

6 

9 

15 

18 

27 

3n~6 

4 

2.0681 






5 

1 

1 .5888 

1.4193 





1 

c 

1.2919 

1 .1540 





7 

1.089 3 

0.97 29 

0.889 2 




8 

0.9419 

0 .8413 

0 .7689 

0.7518 



: 9 

0 43 299 

0.7412 

0.6774 

0.66 23 



10 ! 

0 .7417 

0 -66 24 

0 <»6054 

0.59 20 



11 

0.6706 

0 .5989 

0.5473 

0 .535 2 

0.5155 


12 ’ 

0.6119 

0.5465 

0.4994 

0 .488 3 

0 -.4704 

0 .4669 

1 

13 

0.56 27 

0 .5023 

0 .459 3 

0.4491 

0.4326 

0-4 26B 

14 

(J .5 209 

0.4651 

0-4 251 

0.4157 

0 .4004 

0.3930 

15 

1 0 .4848 

0.4'=»29 

0.3957 

0 .3869 

0 -37 27 

0.3642 

16 

0.4534 

0 .4049 

0 .^701 

0 .3618 

0 .3486 

0.3394 

17 

0 .4 259 

0.3803 

0 -3476 

0.3398 

0.3274 

0.3177 

18 

0.4015 

0 .3585 

0-^276 

0.3204 

0.3086 

0.3987 

19 

0.3797 

0.3391 

0-3099 

0 .3030 

0-2919 

0.2818 

20 

0.360 2 

0.3217 

0-2940 

0-2374 

0.2769 

0 .2667 
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TABLE 

6 -i • 3 • 

Exact critical points 

of the 

test 




for a = 0 

-05 and K 

h 


m 

6 

9 

15 

18 

27 


4 

1 .9 39 0 






6 

1 -4895 

1.3349 





6 

1-2111 

1 -085 2 





7 

1-0 211 

0.9149 

0.8387 




8 

0-88 30 

0-7911 

0-7 25 2 

0*7097 



9 

0-7779 

0.6970 

0-6389 

0 .6 25 2 



10 

0.695 3 

0.6 229 

0.5710 

0 .5588 


! 

I 

11 

0-6 286 

0.563 2 

0.516 2 

0.505 2 

0.48 7 3 

i 

! 

1 

12 

0-5736 

0-5139 

0-4710 

0-4610 

0 .4447 

0 .4415 j 

13 

0 -5 275 

0 .47 26 

0.433 2 

0 -4 239 

0.4039 

0 .4036 

14 

0 .488 2 

0.4374 

0 .4009 

0 *39 23 

0 .3735 

0 .3713 

15 ■ 

0.4544 

0-4071 

0.3732 

0 .365 2 

0.35 23 

0-3446 

16 

0-4250 

0 .3308 

0 .3490 

0 .3415 

1 

0.3 295 

0.3 211 

17 

0 -399 2 

0 .3576 

0 .3 273 

0-3208 

0 .3094 

0-3007 

IS 

0-376 3 

0-3371 

0 .3090 

0 .30 24 

0-2917 

0»QB21 

19 

0.3559 

0-3139 

0.2923 

0.2350 

0.2759 

0 . 2563 

20 

0.3376 

0-30 25 

0 - 277 2 

0.2713 

0*2617 

0 -25 
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TABLE 6 *4 *4 * Cojnparlson of exact and approximate critical 
points of ^^,^2 ^3 a=0^05 and K = 3* 


n 

A 

'^1 

^ i 

f 

^3 

\.X 

Exact 

Approx* 

Exact 

Approx -j 

Exact 

Approx- 

10 

6 

-5012 

*485 3 

•7417 

•7740 

.6953 

•7258 

10 

12 

*44 26 

•4156 

•6 26 2 

•6 371 

•5399 

•5991 

10 

IS 

*4 249 

• 3968 

•59 20 

•6004 

•5588 

•565 2 

10 

24 

*416 3 

• 3881 

-5756 

•58 33 

•5439 

•5494 

20 

6 

• 2435 

-2358 

•3602 

•3758 

•3376 

•35 24 

20 

12 

•2150 

• 2019 

• 3041 

•3094 

• 28 64 

.2909 

20 

18 

• 2064 

*19 28 

•2874 

.2915 

-2713 

• 2744 

20 

24 

-20 22 

•1885 

• 2795 

-28 3 2 

• 2641 

• 2667 

20 

30 

1 

*1998 

•1861 

.2748 

• 2784 

-2599 

• 26 23 

20 ! 

36 

•198 2 

-1845 

• 2718 

• 2753 

-2571 

• 2594 

20 

4 2 

• 1970 

-18 34 

-2696 

*2731 

• 2551 

• 2574 

20 

48 

•196 2 

•18 25 

• 268 0 

• 2715 

• 25 36 

• 2559 

20 

L_ 


*1955 

J 

•1819 

• 2667 

• 2702 

• 25 25 

• 2547 
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TABLE 6 #5 ■ Simulated power values of the tests 

and Vg for a * 0 *05^ K = 3' and 0^ =* • 


(n = ll) 



o 

d = 15 

d = 27 


®2 


^2 

^3 

^1 

^2 

^3 


•0 

-0 

*061 

t05 2 

.054 

•061 

-045 

.045 


♦ 1 

►0 

•114 

.058 

.056 

•113 

-066 

.065 


-2 

•0 

*245 

*101 

•105 

.248 

•106 

•101 


• 3 

.0 

•468 

• 233 

.249 

•504 

.235 

• 257 


♦4 

.0 

•674 

.361 

.413 

*730 

.407 

.450 


*S 

.0 

*844 

.558 

.617 

.897 

.638 

.691 

! 

-6 

.0 

•951 

•7 39 

•776 

*966 

.818 

-864 


•7 

.0 

•98 2 

•864 

•89 2 

*990 

.9 28 

-941 


.8 

♦0 

.99 9 

.949 

.970 

.998 

.978 

.986 


-9 

.0 

.999 

.978 

.989 

.999 

-995 

.996 


• 2 

• 2 

*147 

.106 

*076 

•149 

•119 

.08 3 


•A 

• 2 

•5 37 

• 285 

.309 

*59 2 

• 326 

*34 2 


•6 

-2 

.380 

.60S 

.65 2 

.913 

.679 

.734 


•8 

.2 

.99 2 

•83 3 

*9 20 

*99 3 

.951 

*964 


•4 

.4 

.5 29 

.369 

*304 

.557 

*409 

#345 

1 

1 

1 

1 

• 6 

^4 

.844 

.611 

*656 

i868 

.677 

.711 


•8 

*4 

.97 2 

•871 

.893 

.980 

.911 

*9 35 


.6 

#6 

.850 

•710 

.655 

.871 

.770 

.716 


•8 

•6 

.960 

.881 

*897 

-968 

.9 26 

.9 34 


»8 

•8 

.965 

.913 

*904 

.978 

*940 

#9 25 


1 

J — — 
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APP ENDEC A 


SOLUTION OF THE HL EQUATIONS UNDER 
THE HYPOTHESIS 0^=0 2 =® 

For a given set of ordered sample observations 
5 ; ,,, < x_ _ and v.. < v_ .0 < •*• < v, 


< -,^^2 < -■ < and 


0 0 « ^ y ^ M 

-^ 2-^2 


an algorithm is presented to compute the ML estimates of 

location parameter 0 and scale parameter a under the hypothesis 

6^=6o=Q<’ Id Section 2»6, for x 11 < Yr- >1 ^dd for different 
X X‘2|"rX 

values of r^ and r 2 / the f'lL estimates 0 and a of 0 and a are 
given as follows s 

(a) Case s rj^= 0 /r 2=0 

§ r x^ and a = P/d*, 

(b) Case ; r^ > 0 , d 2=0 


A 

0 


x^ 4 .]^"'^ (l+d^/f) '-dd ff = P/d*, 


(c) Case : r^= 0 ,r 2 > 0 


( i) i^ a < b, then +'l“^ (l+r^f) and a = (P-fO)/d‘''* 

( ii) jf a > b, then 0 = x and a is the solution of equation 

r^+1 

(A«l) exp (Q/o) = 1 + r 2 Q/(P*"^'''^ 

(d) Case : r^ > 0, ^2 > 0 


A 

0 = X, 


ti+l 


'0 log { l+r^Q/(<3.*a+fQ“P)} 



15 2 


A 

ana 0 is the solution of the equation 


r,Q 


(A. 2) exp (0/0) = {1 + — ~}/{l 4 1 


P~d''0 


d*0+fQ~P 




where P = 






X.+S^ X 


i=r^+l 


1 -n,-s.-..i:„ yi+^a^n.-s 


j=r.+i 


2 2 


r^+1' 


Q = Yr^+r^-hl' ^ = f~s^~S2, 

a = Q/log( l-i-i'2/f ) and b = (P~fQ)/d*- 

As remarked in Section 2-6, equations (A.i) and (A«2) 
are solved by Newton- Raphson method with (2P-fQ)/2d* as an 
initial value* 

The subroutine ESTMAT computes the ML estimates THETA 
and sigma of 9 and 0 respectively foi~ given values of n. ,n-,r^,r„, 
and Sjy and vectors X and Y. E’or this subprogram, the 
r^iquirod accuracy AC, and txn^ expected num.ber of iterations NI 
aro supplied from the main calling program* The failure 
indicator PI takes the value 1, if the itv^ration procesdure does 
not converge in NI steps* 

LANGUAGE 


Port ran 10 



IbJ 


STRUCTURE 


SU Bi WT INE ESTMAT ( N1 , N 2 , 

Pqraal parameters 

^^1/^^2/S1,S2 

X, Y, NI; AC, THETA, SIGMA , PI ) 

N1 

Integer 

input 

m 

size of the first sample 

N2 

Integer 

input 


size of the second sample 

R1 

Integer 

inpu.t 

• 

number of smallest 
observations missing in the 
fii'st sample 

R2 

Integer 

input 


number of smallest 
observations missing in the 
second sample 

SI 

Integer 

input 

o 

number of largest observations 
missing in the first sample 

S2 

Integer 

input 

o 

number of largest observations 
missing in the second sample 

X 

Real 

input vector of : 
length Nl-Ri~si 

available ordered observatio is 
in the first sample 

Y 

Real 

input Vector of : 
length N2-R2-S2 

available ordered observations 
in the aecond sample 

WI 

Integv^r 

input 


upper bound for the number of 
iterations in i^iiich the 
iteration pirocoss is expected 
to converge 

AC 

Real 

input 


desired accuracy 

THETA 

Real 

output 

0 

ML estimate of 0 

SIGMA 

Real 

output 


ML estimate of 0 

FI 

Integt-'r 

output 

7 

failure indicator 


1 if the iteration does not 
converge 


0 otherwise* 
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C THIS SUBROUTINE COMPUTES THE ML ESTIMATES OE THETA 

C AND SIGMA BASED ON TYPE II DOUBLY CENSORED 

C samples prom TWO EXPONENTIAL DISTRIBUTIONS 

C 

SUBROUTINE ESTMAT ( Nl^ N2, Rl/ R2, SI, S2#X, Y, NI, AC,THETA^ 

1 SIGMA, PI) 

INTEGER R1,R2,S1,S2,FI 

DOUBLE PRECISION F, D, P, Q, A, B,X(999 ) , Y(999 ) , Z(999 ) , S( 999 ) , 
ITHETA, S IGMA, SUM, AC, U1 , U 2, U"^, U4, FN, FP 
IP(X(Rl+l) .LE,Y(R2+l) ) GO TO 40 
C INTERCHANGING THE SAMPLES 

C 

N11=Ni;NR1=Ri;NS 1=S1 
N1=N2/R1=R2/S1=S2 
N 2=N1 1 ; R 2=NR1 ; S 2=NS 1 
DO 10 I=R1+1,N1-S1 
10 Z(l)=Y(l) 

DO 20 I=R2+1,N2-S2 
20 Y(I)=X(I) 

DO SO I=R1+1/N1~S1 
30 X(I)=Z(I) 

C CALCULATION OF F,D Q 

C 

40 F=Nl-fN2~Rl-R2 

D=:F-S1~S2 
Q =Y( R2+1 )-X( Rl+1 ) 



SUM OF Ar.L OBSERVATIONS 


SUM=0 

DO 50 I=R1+1,N1~S1 
SUM=SUM+X( l) 

DO 60 I=R2+1,N2-S2 
SUM=SUM4-Y( I) 

CALCULATION OP P 

P=SUM+Sl-rt-X(Nl-Sl )+S2^Y(N2~S2)->F*X( RI+1 ) 
calculation of a and b 

A=Q/DL0G( 1 .+RV5’) 

B=-.( P-Pi^Q )/D 

SEP.-iRATI0N OF DIFFERENT CASES 

IF(R2.ED.0) go TO 75 

IF( (R1.EQ.0) .^^D.(A.LE„3)) GO TO 80 
£STIMi\TING SIGMV BY NEWTON-RiiPHSON METHOD 
WITH INITD^L V.iLUE S(l) 

S ( l) =( 2 )/ ( 2 .'U3 ) 

1=1 

Ul=l .+R2^«-Q/ ( P-Da-S(T) ) 

U2=l »+Rl*Q/(D*S( I) +Ri-Q“'P) 

U 3 =R 2'>'rQ'^D/ ( ( P~D ( I ) 2 » ) 

U4=Ri^Kl4tD/( (D->tS( l)+BS-Q~P) ■)«J-2* ) 
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C THE FUNCTION OP SIGM^i IS DENOTED BY^PN" 

C ^-.ND ITS DERIVxiTIVE IS BY^'PP" 

C 

FN=DEXP(Q/S( I))-U1/U2 

PP=^DEXP(Q/S( I) )4J-Q/( S( l)^H«-2*)-(U3^c-U2+Ul*U4)/(U2 w2*) 
S(l+l)=S(l)~FN/FP 

IP(DABS(S(l+l)-S(l)) .LE.AC) GO TO 70 
IP(I-GE.NI) GO TO 85 
1 = 1+1 
GO TO 65 

C FIN/iL V.JLUE OF SIGMTi. AND THETA 

C 

70 SIGKTi=S( I+l) 

THBT.i=X( Rl+1 )~SIGM/#-DL0G( 1 ,+Rl*Q/ ( D*SIGMA+F*Q~P ) ) 

GO TO 90 
75 SIGm=p/D 

THErA=:X( Rl+1 )-SIGI'iS^'tDLOG( 1 .+R1/P ) 

GO TO 90 

8 0 SIGM^'i=(P~B5-Q)/D 

THETA=Y( R2+1 )~SIGHi*DLOG( 1 #+RVF) 

GO TO 90 

C /ASSIGNING THE VI^UE TO THE FAILURE INDICATOR 

C 

85 FI=1 

90 RETURN 


END 
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C 

Q **■}(•* -;t- * * # 

CALLING PROGRAM 

C -Vr ^ » -,!r -/• ■JS' -Yf * 

c 

C SOLUTION OF THE ML EQUATIONS 

INTEGER R1,R2,S1,S2,FI 

DOUBLE PRECISION X( 3.00 ) , Y( 100 ) ^ AC, THETA^ SIG^ 

C INPUT VALUES, THE DATA X(l) AND Y(j) ARE SIICTLATED 

C VALUES FROM E(2, l) DISTRIBUTION ARRANGED IN AN ASCENDING 

C ORDER OF l/t^GNITUDE 

C 

Data N1,N2,R1,R2,S1,S2,NI,AC, (X( I), 1=4, 14), (Y(J), J=3,12)/ 

1 17, 13, 3, 2, 3, 1, 20, *00001, 2.16140, 2.21918 , 2.30073 , 2.34808, 
1 2.91879, 3,14132, 3 .21995, 3 .347 28 , 3 ,35 224, 3 #45513, 3 .6 2251, 
1 2.15 214, 2-.18 279 , 2*29610, 2*30496 , 2 *41608 , 2 .53415 , 2 *5 7 010, 
1 2*95 275 , 2*968 00,4 *3 2659/ 

CALL ESTM.iiT( Nl, N2, Rl, R2, Sl, S2, X, Y, NI,AC, THETA, SIGMZ^, FI ) 
IF(FI~l) 10,30,10 
10 P RINT 20 , THET/i, S IGM^\ 

20 FORMAT dOX,' THE ESTm^TE OP THETA IS^ F9 *5/l0X, ' THE 

1 ESTIMATE OF SIGM/\ IS',P9*5) 

♦ 

GO TO 50 

30 PRINT 40, NI 

40 FORMAT ( 5 X, ^ THE TTER^iTION DOES NOT CONVERGE IN', l4,^ STEPS'' ) 

50 STOP 

END 

OUTPUT 


THE ESTIM'VTE OP THETA IS 1*979 21 
the ESTm-*TE OP SIGM/i IS 0.97834 



appendix b 


EVALUATION OF Q^(xts)AND L^(xls) 

The procedui'es for evaluating the functions Q^(xls) and 
Ii^(xls), where d is a positive Integer are presented. Function 
subprograms QD(D^X, s) and LD(D,XyS) are given for calculating 
these functions# 

From Section 3*1, we have 

Q^(xis) =/ dy/(d~l)l , x > 0. 

X 

Note that, this integral is convergent only for (l+s) > 0, that 
is for s > -1- By substituting z = y(l+s), we get 

Qjj(xis) = / e ^ dz/{ (l+s)^(d~l)i} . 

x( 1-i-s ) 

(B.l) = g-xCl+s) , s > -1. 

j=0 

Similarly, for finite positive values of x, we have 

L^(x|s) =/ y'^-l dy/(d-l)!. 

o , 

This converges for all s, but has to be treated separately for 
s = 1 and s /: 1 cases- Thus 

:ic^/dl , s = 1 

(B.2) L^(xis) = ^ 

f 1- S e~^^ x( l~s)} V-3 !]/*^ ^ ' s/: 1- 

I'- j=0 



Tho FUNCTION QD(D, X, s) computes the summation given 
in equation (B-l) for D = 1,2,.../ X > 0 and S > ~i. For 

= 1,2, »../ X > 0 and < S < oo, FUNCTION LD(D,x, s) computes 
the value of L^(x|s) given in equation (B.2). 

LANGUAGE 

Fortran 10 

C FUNCTION SUBPROGRAM FOR COMPUTING QD(X|S) 

C 

DOUBLE PRECISION FUNCTION QD(D,X,S) 

DOUBLE PRECISION SUM, X, S, XI, FAC 

INTEGER D 

SUM=0 .0 

DO 10 J=0,I>.1 

XI =1 .0;I=J 

IP(J-EQ.0) GO TO 90 

Xl^X*! 1 .+S) 

IF(J.EQ.l) GO TO 90 
DO 80 K=2,I 
pac=dfloat(k) 

8 0 XI =Xl^t.X* ( 1 .+S )/FAC 

9 0 SUM=SUM+X1 

10 CONTINUE 

QD=DEXP(-.X.jf( 1 .+S) )* SUM/ ( ( 1 ,+S )^H-D) 

RETURN 


END 



c FUITCTION SUBPROGRAM FOR CCMPUTING LD(XtS) 


DOUBLB PRECISION FUNCTOT LD(D,X,S) 
DOUBLE PRECISION X^ S, SUM, XI, FAC, SI 
INTEGER D 

IF(S.EQ.l) GO TO 21 

SUM=0 »0 

3X) 15 J=0,D-1 

XI =1 -o;i=J 

IP(J*EQ.O) GO TO 90 
Xl=3Si- ( 1 .~S ) 

IF(J.EQ.I) GO TO 90 
DO 80 K=2 , I 
fac=dfloat(k) 

3 0 XI =Xl-;i-X.ii-( 1 *“3 )/FAC 

90 SUM=SUM+X1 

^5 CONT INUE 

1,D=( 1 .-de:.cp(-x*( 1 -s))*SUM)/( ( 1 

GO TO 28 

21 Sl=l- 

DO 25 J=l/D 

FAC=D FLOAT ( J ) 

25 S1=S15'<-’V^’^^ 

ld=si 

o-a return 


END 



